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Exchange Forces and Electrostatic Forces Between Ions in Solution 


RonaLp W. GuRNEY 
Wills Physical Laboratory, University of Bristol, England 
(April 5, 1938) 


Though the treatment of ionic dissociation by thermo- 
dynamic methods has been fruitful, the mechanism can be 
better understood by the use of statistical mechanics. For 
this purpose one needs to know the potential energy of the 
ions as a function of temperature. The work D required 
to dissociate a molecule or molecular ion consists of two 
parts, one part Dg due to electrostatic forces, and the 
other Dyon due to nonelectrostatic forces. As the former 
varies rapidly with temperature, for each species the value 
of dD/dT is a guide to the value of D,). It is shown that 
for acid molecules and molecular ions in aqueous solution 


the relative values of dD/dT for different types are in the 
order predicted. Further, the characteristic temperature @ 
at which the dissociation constant passes through a 
maximum should depend, not on the value of D, but on 
the relative magnitudes of Dnon and Dg). It is found that 
the known values of 6 fall into line in accordance with this 
rule. Between two Hg* ions, or two Te™ ions, there will be 
exchange forces of attraction similar to those which give 
rise to the molecules Na2 and Br2; thus in aqueous solution 
the ions (Hgs)** and (Tez)~~ are stable, although for 
them Dg) is negative. 


$1 


LTHOUGH the treatment of gaseous mole- 
cules by the methods of quantum mechanics 

has now become habitual, little has been done to 
use quantum-mechanical forces to describe the 
stability and the dissociation of molecules in 
solution. The author has described! the mecha- 
nism which determines the heat of ionic dissocia- 


_ tion, and suggested that the details of dissocia- 


tion at different temperatures may be under- 


_ stood by recognizing the magnitude of the ex- 
_ change forces between pairs of simple ions. 


_ Properties of benzene derivatives, 


Similar ideas have been used, particularly by 
L. P. Hammett,? in accounting for the effect of 
various substituents on the thermodynamic 
etc., in 


aqueous solution. 


, Gurney, Ions in Solution (1936), Chapter 13. 
Hammett, J. Chem. Phys. 4, 613 (1936); J. Am. Chem. 
Soc. 59, 96 (1937). 


A greater insight into the mechanism of such 
processes than is possible by a thermodynamic 
treatment, can be obtained by the use of statis- 
tical mechanics. For this purpose we require in 
the first place a knowledge of the mutual po- 
tential energy of the ions. Take then a positive 
and a negative ion first in a vacuum. Leaving 
aside the Coulomb forces between them, we 
know that the nonelectrostatic forces will 
contribute either (a) a repulsion only, or (4) an 
attraction as well as a repulsion; that is to say, 
this contribution, plotted against the separation 
r, will have the form either of curve a in Fig. 1 
or of curve b. To obtain the total potential 
energy, this contribution must be combined with 
the electrostatic energy. If the latter is given by 
curve c, on adding the ordinates to those of 
curve a, curve d is obtained; and similarly from 
curve b curve e is obtained. 

Between two positive ions or two negative ions 
there will be an electrostatic repulsion, repre- 
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sented by curve f in Fig. 1. But in some cases 
there are in addition strong exchange forces of 
attraction similar to the familiar forces between 
neutral atoms. In the diatomic molecules Naz and 
Kg it is well recognized that the exchange forces 
of attraction are almost as simple as in the 


iz P 


Fic. 1. 


hydrogen molecule He, the single valence elec- 
tron in each atom taking part. Consider now a 
singly charged calcium ion Ca*; it has the same 
electronic configuration as the neutral atom of 
potassium, namely, with one 4s electron outside 
the core ; only the charge on the atomic nucleus is 
one unit higher. If then we bring together two 
Ca* ions, exchange forces of attraction will be 
present. But, although two K atoms combine 
to form Ks, two Cat ions will not in a vacuum 
combine to form (Ca)2** owing to the electro- 
static repulsion between the nuclear charges. 
If we had some means by which we could cut 
down the electrostatic repulsion without at the 
same time cutting down the strength of the ex- 
change forces of attraction, there would be the 
possibility of obtaining such diatomic ions as 
(Caz)*+* and (Hge)*+* etc. Now a polar solvent 
has just the necessary properties. The orienta- 
tion of the polar solvent molecules causes a 
large effective reduction in the Coulomb forces. 
Exchange forces, on the other hand, arise from 
the rapid circulation of electrons—motion which 
is so rapid that the heavy solvent molecules can 
make no response to it. If curve f in Fig. 1 repre- 
sents the potential energy of a weak electro- 
static repulsion, then on combining this with 
curve b we obtain curve g. In aqueous solution at 
room temperature the familiar (Hg2)** ion must 
possess a potential minimum of this kind suffi- 
ciently deep for stability (confirmed by the 
presence of a Raman spectrum). 
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Similar exchange forces between halogen atoms 
give rise to the molecules Fe, Cle, Bre and I,. 
Now the oxygen atom has one electron less than 
fluorine, and the atoms of sulphur, selenium and 
tellurium have one electron less than Cl, Br and 
I, respectively. The singly charged negative ions 
O-, S-, Se~ and Te~ have the same electronic 
configurations as the halogen atoms, and _ be- 
tween any pair of these ions exchange forces of 
attraction will be present. In any solvent of high 
dielectric constant the electrostatic repulsion 
may be so small that diatomic ions with two 
negative charges are stable. As long ago as 1923 
it was, in fact, discovered* that in aqueous solu- 
tion at room temperature tellurium ions exist in 
the diatomic form (Tez)~~. 


§2 


In the remainder of this paper we are concerned 
with the potential energy between ions bearing 
charges of opposite sign, for which curves d and 
e in Fig. 1 give the appropriate forms. Let the 
depth of the minimum of curve e be denoted by 
D. In curve d the depth of the potential minimum 
is due entirely to electrostatic forces of at- 
traction ; but in curve e these electrostatic forces 
are superimposed on curve b, which already has 
a potential minimum due to the nonelectrostatic 
forces. We may therefore write 


D= Dponeiectrostatic t+ Detectrostatie: (1) 


Species of ions for which the potential energy has 
the form of curve d in Fig. 1 may form short 
range ion pairs. It is only for species for which 
the potential energy has the form of curve ¢ 
that we can draw a distinction between stable 
molecules in solution and ion pairs; the depth D 
of the minimum is the dissociation energy of the 
molecule in solution. The importance of separat- 
ing this energy into the two parts D=Dyont+Da 
lies in the fact that exchange forces between 4 
pair of atoms or ions are characteristic of the 
species concerned and are independent of 
temperature; while D.; depends on the nature 
and temperature of the solvent. With rise of 
temperature the curve c of Fig. 1 becomes steeper, 
and the depth of the minimum becomes greater, 


3 Karsanowsky, Zeits. f. anorg. allgem. Chemie 128, 
17 (1923). 
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the increase in D being equal to the increase in 
D... In a solvent of dielectric constant K we 
shall not expect the increase in D,, to be strictly 
proportional to 1/K owing to the familiar elec- 
trical saturation of the dielectric in the neigh- 
borhood of each ion. 

Although the methods of this paper are ap- 
plicable to the dissociation of all species of 
molecules and molecular ions, the only types for 
which recent accurate data are available are acid 
molecules and ions, where a proton (or deuteron) 
is transferred from the solute to a solvent 
molecule; we shall accordingly discuss only 
these. Some years ago Wynne Jones‘ examined 
the data on the strength of acids in different 
solvents, and found that they grouped them- 
selves, as we should expect, into three classes, 
according as the proton is transferred from (1) 
a neutral molecule, (2) a negative ion, such as 
(H2PO,)-, or (3) a positive ion such as (NH,)*. 

Since then successful attempts have been 
made to interpret the effects of polar substituents 
at room temperature.’ Hammett® has also dis- 
cussed the change of entropy with temperature, 
recognizing clearly the separate contributions 
from the electrostatic and nonelectrostatic 
forces. 

To use the methods of statistical mechanics, we 
need, in addition to the potential energy, a 
knowledge of the energy levels of the protons. 
It is well known that in molecules such as HCl 
the spacing between the vibrational levels is so 
large that even at temperatures above room 
temperature practically all the molecules are in 
their lowest vibrational level. This is due to the 
small mass of the proton and will be true of the 
protons with which we shall be concerned here. 
Suppose first that in a vacuum we have a mole- 
cule of an organic acid ROOH and a water 
molecule side by side; in the ROOH molecule 
the proton is in a vibrational level inside the 
electron cloud of one of the oxygen atoms. 
The water molecule is capable of accommodating 
an additional proton, to form (H;O)+; ice., in 
the water molecule is a set of vacant vibrational 
levels for a proton; we are concerned only with 


* Proc. Roy. Soc. A140, 448 (1933). 
Schwarzenbach and Egli, Helv. Chim. Acta 17, 1183 
(1934) ; Schwarzenbach and Epprecht. ibid 19, 493 (1935). 
Hammett, reference 2, p. 101. 
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the lowest vibrational level of the set. If the 
vacant level in the HO molecule is higher than 
the initial occupied level, in order to transfer the 
proton from one to the other, a certain amount of 
work will be required. When the proton has been 
transferred, leaving behind a negative charge on 
the oxygen atom, an electrostatic field has been 
created in the surrounding vacuum, having a 
density of energy equal at every point to E*/87, 
where E is the intensity of the field; for a mole- 
cule in solution similar considerations will 
apply. The work required to complete the 
separation of the ions is entirely electrostatic, 
and D., finally becomes equal to the integral 
of the energy in the field surrounding the nega- 
tive ion added to the energy in the field sur- 
rounding the (H;O)* ion. The latter part is of 
course the same for all acid molecules. From the 
usual Born formula,’ for a molecule which 
dissociates into a pair of singly charged ions of 
radius b, and b_ we should, in the absence of 
electrical saturation have 


Da=—(— 
2K\b, 
For a pair of small ions in water at room tem- 
perature (2) has a value in the neighborhood 
of 0.1 electron volt. To give the low degrees of 
dissociation of the familiar weak acids, values of 
D greater than 0.2 electron volts are required ; 
for most species Don is greater than D,,. (If the 
neutral atid molecule is surrounded by a certain 
field, because it is polar, the value of D., will 
more strictly be the difference between the final 
and initial energies of the field.) 
If we consider next processes of the type 


(H;0)+ 
(H;0)+ 


(2) 


(3) 


we see that they are different from the dissocia- 
tion of a neutral molecule, because initially an 
electrostatic field exists round the positive ion, 
and the positive charge is merely transferred 
elsewhere. If the total energy of the field round 
the (H;O)* ion is, for a particular species, the 
same as the total energy of the field initially 
surrounding the positive ion, the value of D,, 


7 Born, Zeits. f. Physik. 1, 45 (1920) ; compare Hammett, 
reference 2, p. 102. 
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for this process will be zero. It will be convenient 
to name this an isoelectric transference of a 
proton. For dissociations similar to (3), even 
when they are not approximately isoelectric, the 
value of D.; will at any rate be much smaller 
than in those dissociations where we create a 
positive and a negative charge. 

Intermediate between these two types lie 
dissociations like that of the glycine ion: 


NH 3*CH + . 
NH 2COO-+H 307. (4) 


As in (3) we begin with a positive ion, and end 
with a molecule which as a whole is electrically 
neutral, but which here bears opposite charges 
on its two ends; their electrostatic field may 
possess a considerable energy, though not, of 
course, as much as a pair of separate ions. The 
same applies to the other mode of dissociation of 
glycine: 


H.0+NH:CH:COO-— 
OH-+NH;*CH:COO-. (5) 


In this case the electrostatic work of transferring 
the proton should be diminished by the negative 
charge already on the glycine ion, which will 
attract the proton. In (4) above the electrostatic 
work is diminished by the repulsion of the posi- 
tive charge on the glycine ion; we expect then a 
smaller D,, than that for a neutral ROOH mole- 
cule, as already stated. On the other hand, in the 
dissociation of a negative ion, for example 


we expect D., to be greater than for a neutral 
molecule. 

In the next section we shall examine the rate of 
variation of D with temperature for these four 
types of dissociation. We shall expect no correla- 
tion between this rate and the value of the total 
D. For a large or small value of D.; may be ac- 
companied by a large or small Dion; and it is 
only D,, that changes. Over the range of tempera- 
ture from 0°C to 60°C for which measurements 
are available, the dielectric constant of water 
falls from 88.31 to 66.62;% the value of 1/K 
rises from 0.01133 to 0.01502, that is, by 33 
percent. If then the electrostatic energy were 
strictly proportional to 1/K, the value of D.) 


* Akerlof, J. Am. Chem. Soc. 58, 1842 (1936). 
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would also increase by 33 percent. Though this 
will not be so, we must expect that dD/dT will 
be large when D,, is large, and small when D,, is 
small. 


§3 


In expressing the dissociative equilibrium at 
temperature 7, to ensure that the partition 
functions of the various particles are measured 
from the same zero of energy, the factor e~? *” 
must be inserted. The relation between the 
activities of the positive and negative ions, a. 
and a_, and the activity of the neutral molecule 
a, is then 


a,a_/a=K=Ce?l'r, . (7 


where the statistical factor C depends on the 
solvent only, and is the same for all solutes. 
(In the second dissociation of an acid an addi- 
tional statistical factor } enters.) For substances 
in aqueous solution at room temperature it is 
well known that for some species the dissociation 
constant K increases with rise of temperature, 
and for others falls. The description given here 
accounts for the characteristic behavior of each 
type. For (7) takes the form 


—In (8 


For a dissociation which is approximately iso- 
electric the last term will be absent, and § 
is found to take the form 


—log K=A/T+B, 9 


where A and B are constants.’ For neutra 
molecules, on the other hand, this will not be so 
For five species of neutral molecule Harned anc 
his co-workers find that the value of K passe: 
through a maximum value at a certain tempere- 
ture @ characteristic of each substance, and the 
variation of log K can in each case be representec 
approximately by a parabola 


with the same value of p for all.!° It will appea’ 
that this approximate regularity of behavior © 
due to the change in the dielectric constant of the 
water in which all the substances are dissolve¢ 

* Pedersen, Dansk. Vidensk. Selsk. 14 (1937): als 
measurements for the (NH,)* ion over the range irom ° 
to 45°C will be published shortly by Wynne Jones. 


10 Harned and Embree, J. Am. Chem. Soc. 56, 1 
(1934). 
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Consider any two substances in aqueous solu- 
tion at temperature 7. Let their dissociation 
constants be K and K,, respectively; and at 
temperature 7” let the values be K’ and K,’. 
Then from (7) we shall have the equations 


D'—D=kT K—kT" In K’ 

—kT \In C+kT’' In C’, (11) 
D/—D,=kT K,—kT’ In K,’ 

—kT |In C+kT" In C’. 


On subtracting the second from the first, we 
obtain 


(D'—D) —(D,’—-D,) 
=kT |In (K/K,)—k7” In (K’/K,’). (12) 


If now the dissociation to which K, refers is 
nearly isoelectric, while the other is not, the 
value of (D,’—D,) will be small compared with 
(D’—D). The variation of K, will be given by an 
expression of the form (9) above. If, taking the 
observed values of K for several other species of 
molecules, we substitute them in (12), we can ob- 
tain values of the change in D for each of them 
relative to a single reference substance (or 
relative to each other). This has been done for 
different types in Table I." Similarly the rela- 
tive values of dD/dT at 25°C have been ob- 
tained from the expression 


(d/dT)(D—D,) =(d/dT)(kT In (K,/K) (13) 


using the empirical power series given in the 
original papers. Apart from the figures in the 
last column, the values have been calculated 
taking log K,=A/T+B, with B equal to 0.614, 
the value found recently by Pedersen® for the 


Tama 

D) D 

(d/dT)(D —D,) at 25°C 

Temp. VOLTS €-VOLTS PER 1°C 

(H:PO,)- 20°-50°C | 0.0414 in 30°} 13.1107! | 
Acetic acid 0°-60° | 0.0662 in 60°} 10.8 13.1 
Glycine Kp 10°-45° | 0.0205 in 35°} 5.5 7.8 
Glycine K4 10°-40° | 0.0126 in 30°} 4.2 65 
Anilinium (0.00) (0.0) 2.3 


"I have not included values for the (HSO,)~ ion, de- 
rived from the measurements of Hamer (J. Am. Chem. Soc. 
56, 860 (1934)), because in Chicago I learned from Pro- 
essor T. F. Young that a calorimetric measurement 
of the heat of dissociation had yielded a value twice as 
great as that given by Hamer; a still higher value has 


recently bee i i 
2368 1037) published by Pitzer (J. Am. Chem. Soc. 59, 
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TABLE IT. 
D IN ELECTRON 
Temp VOLTS dD/dT at 25°C 
(H2PO,)~ 50°C 0.534 15.3 X1074 
0° 0.486 
Acetic acid 60° 0.433 13.0 
o° 0.353 
Glycine Kg 45° 0.369 7.7 
10° 0.341 
Chloracetic acid 40° 0.289 10.5 
0° 0.247 
Glycine Ka 40° 0.252 6.4 
10° 0.233 
H3sPO4 50° 0.257 10.2 
0.3° 0.206 
Water * 60° 1.090 15.1 
o° 0.998 
HCl estimate 25° —0.31 


* Harned and Hamer, J. Am. Chem. Soc. 55, 2194 (1933). 
+ Ebert, Naturwiss. 13, 393 (1925); Wynne Jones, J. Chem. Soc. 
1064 (1930). 


anilinium ion. In the last column of the table 
are given for comparison values calculated 
taking B equal to 1.772, the value found by 
Pedersen for the o-chlor-anilinium ion.” 

It will be seen that for the different types the 
values lie in the order predicted above. The most 
rapid variation with temperature is found for 
the negative ion H,PO,-. Next comes the neutral 
molecule of acetic acid; and the values for the 
other similar acids fall here. Then the glycine 
values lie between acetic acid and the anilinium 
ion, as they should. 


§4 

As these results derived from (7) are satis- 
factory, we may discuss (7) further. When disso- 
ciation of a molecule or ion has taken place, the 
result is that the proton is now vibrating in the 
interior of a water molecule instead of in the 
interior Of the original molecule or ion. In a 
solution we need to know the ratio of the number 
of available levels for the proton in the solvent 
molecules to the number of levels in the solute; 
this ratio determines the statistical equilibrium. 
As we are concerned only with the lowest vi- 
brational level in each case, there is one level 
for the proton in each water molecule. The ac- 
tivities a, and a_ are usually expressed in moles 
per liter; since one liter of water contains 55.5 
mole of water, the value to be assigned to the 
factor C in (7) according to the above simple 
scheme is 55.5 for the first dissociation of any 
acid. For the second dissociation of a dibasic 


Pedersen, Dansk. Vidensk. Selsk. 15 (1937). 
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acid, the value will be one-quarter of this, that 
is 13.9. As it seems worth while at least to explore 
the simplest assumptions, values of D calcu- 
lated in this way from experimental values of K 
are given in Table II. If 55.5 is held to represent 
merely an upper limit to the number of available 
levels in the water molecules, each value of D 
in the table must be regarded as an upper limit. 

It will be noticed that the values of dD/dT are 
almost identical with the corresponding values in 
the last column of Table I. The reason for this is 
that Pedersen’s empirical expression for the 
o-chlor-anilinium ion, used above, may be 
written K=59.2e-?/*7, Or, in other words, use 
of the value C=55.5 is equivalent to expressing 
results relative to a reference K, given by 
log K,=A/T+log 55.5. 

Included in the table are values for the 
self-dissociation of water: 2HxO—H;0++OH-, 
which has the equilibrium 


(14) 


where a’ has been taken equal to 55.5*=3080. 
The very small dissociation of water must be due 
to the presence of a large D,.n; there is no reason 
why D,, should be very different, since the disso- 
ciation produces a (H;O)* ion and a negative ion, 
as from any neutral molecule. As we are taking 
the magnitude of dD/dT as a guide to the value 
of D.1, we expect the value of dD/dT for water 
to be similar to that for acetic acid, or a little 
larger, since the OH~ ion has a small radius and 
the energy of its electrostatic field may be some- 
what larger than for a (ROO)~ ion. It is satis- 
factory to see from Table II that this is so. 

We have hitherto spoken of the dissociation 
energy D as if it were necessarily positive. There 
is, however, no difficulty in imagining a neutral 
acid molecule in which the energy level of the 
proton is higher than the vacant level in the 


TABLE III. 
—log K at 25° 3) dD/dT at 25° 

(1)* HsPO. 2.123 —18.0°C | 80 X10-4 
(2) Glycine K, 2.350 +53.9 4.2 

(3) Chloracetic acid 2.861 — 54 8. 

(4) Formic acid 3.752 +24.7 8.7 

(5) Glycine Kg 4.219 93.0 5.5 

(6) Acetic acid 4.757 22.6 10.8 

(7)t n-Butyric acid 4.817 8.0 11.7 

(8) (H2PO,)- 7.206 43.1 13.1 


* Nims, J. Am. Chem. Soc. 56, 1110 (1934). 
+ Harned and Sutherland, J. Am. Chem. Soc. 56, 2040 (1934). 


solvent molecule. For such a substance in dilute 
solution nearly all the protons will have been 
transferred to solvent molecules, the small num- 
ber remaining behind, to form acid molecules, 
being in accordance with the Boltzmann law. 
The electrostatic part of D is necessarily positive; 
but, when D,on has a large negative value, the 
total dissociation energy D may be negative. 


x 
/ 
Fic. 2. 


For a particular species of acid molecule the 
value of D may clearly be positive in one solvent 
and negative in another. For illustration a 
rough value of D for HCl in water has been 
included in Table II, based on estimates of K. 


§5 


The variation of D.; with temperature is 
intimately connected with the characteristic 
temperature @ at which K becomes a maximum. 
As Table III shows, there is no correlation at all 
between the degree of dissociation of a substance 
and the value of @ assigned to it. The ideas out- 
lined in this paper may therefore be subjected 
to a severe test by enquiring whether they ac- 
count for these seemingly haphazard values. 

When with rising temperature the value of 
K reaches a maximum, and begins to fall off, this 
is because the potential minimum of curve é¢ in 
Fig. 1 is becoming rapidly deeper. Eq. (8) may 
be put into the form 


—log K/C=A/T+f(T), 


where A =Dyon/k, and f(T)=Dei/kT. In Fig. 2 
the principal intersecting curves are intended 
to represent —A/T and —f(T), respectively, for 
a particular species of molecule, e.g., acetic acid. 
Adding the ordinates of these two curves, oné 
obtains the lower curve with a maximum, which 
is Harned and Embree’s parabola. The diagra™ 
shows clearly how the characteristic temper 
ture 0, at which the maximum occurs, depend: 
upon the relative magnitudes of the electrostatic 


and 

ener 
mole 
stati 
curv 
will 

addi 
curv 
the s 
the 1 
The 
acids 
tion 
abou 
aceti 
the 1 
to t 
disso 
stati 
acid, 
tion 

aD /c 
dD. ¢ 
value 
the | 
hand 


lieve 
excep 
93°C 
In 
ion w 


ilute 
been 
1um- 
ules, 
law. 
tive; 
, the 
tive. 


FORCES BETWEEN IONS IN SOLUTION 505 


and nonelectrostatic parts of the dissociation 
energy. In comparing two different species of 
molecules, if we have nearly the same electro- 
static forces but weaker exchange forces, a 
curve for Dyon/kT like the broken line in Fig. 2 
will be used for the latter. It is obvious that, on 
adding the ordinates, the maximum of the lower 
curve will be shifted upward and to the left; i.e. 
the substance will be more highly dissociated and 
the maximum will occur at a lower temperature. 
The behavior of phosphoric and chloracetic 
acids, relative to acetic acid, affords an illustra- 
tion of this. Their values of K are respectively 
about 500 times and 100 times as large as for 
acetic acid, while @ lies below the freezing point; 
the values —18° and —5.4° have been assigned 
to them, respectively. If the high degree of 
dissociation were due mainly to weak electro- 
static forces, 8 would be higher than for acetic 
acid, as is clear from Fig. 2. But, as the dissocia- 
tion creates a positive and a negative ion, we 
expect a normal D,;, and indeed the value of 
dD/dT shown in Table II indicates a value of 
D.. differing little from the other acids. The low 
value of D is therefore due to a small Dron, and 
the low value of @ is anticipated. On the other 
hand, for both dissociations of glycine we be- 


Fic. 3. 


lieve the electrostatic forces to be small. Their 
exceptionally high values of 6, namely 53° and 
93°C are therefore satisfactory. 

In the further dissociation of the (H2PO,)- 
1on we have strong electrostatic forces, and these 


by themselves would make the value of @ lower. 
The observed value of @ is, however, 43.1°, as 
compared with 22.6° for acetic acid. This is not 
in disagreement with the theory but affords 
additional evidence in its favor; for the value 
of K is some 300 times smaller than for acetic 
acid, showing the presence of a much larger Dyon, 
which in Fig. 2 shifts the maximum of the curve 
downwards and to the right. This behavior may 
be contrasted with that of the bisulphate ion 
(HSO,)-; here, in spite of the large D.i, the 
value of K is nearly 1000 times greater than that 
of acetic acid; we should expect an exceptionally 
low value of 6; and, indeed, Harned assigned the 
value —27.6° to it, a temperature lower than for 
any other substance investigated. 

The value of 6 for any substance thus depends 
upon the relative magnitudes of the two parts of 
the dissociation energy Dron and D,). But it 
would perhaps be more convincing to express 
this in terms of experimental quantities, with- 
out mentioning D,., and D.;. Now for differ- 
ent species we have been taking the value of 
(d/dT)(kT \In (K,/K)) as a relative measure of 
Di; and for different species the value of 
log K is a relative measure of the total D. Sup- 
pose then that for each dissociation we define a 
quantity x by 


x=(—log K)/(dD/dT), (16) 


both terms being taken at room temperature, 
then x is a measure of the ratio of D to Dx, and @ 
should be a smooth function of x (though it 
shows no correlation with either the numerator 
or the denominator taken separately). The suc- 
cess of this maneuver, shown in Fig. 3, is 
satisfactory. In the figure the numbering of the 
points corresponds to the numbering of the sub- 
stances in Table III, where they are arranged in 
the order of log K. This treatment should have 
equal success for substances in nonaqueous 
solution, when sufficient data are available. 

The simplified discussion given in these last 
two paragraphs is merely a step towards a more 
elaborate statistical treatment. 
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Bjerrum's theory of the influence of substituents on dissociation constants has been extended 
and amplified. The molecules and ions entering into the ionization equilibria are treated as 
cavities of low dielectric constant, rather than as structureless regions of the same dielectric 
constant as the solvent. The theory gives better results than the simple Bjerrum formulation, 
especially for the short chain dicarboxylic acids, and in the fact that it permits a satisfactory 
treatment of the influence of dipolar substituents on dissociation constants. 


I 


HE ratio of the first to the second dissocia- 
tion constant of a dibasic organic acid is 
always greater than four, and approaches four as 
the length of the acid is increased. In 1923 
Bjerrum! suggested that the ratio of the dis- 
sociation constants can be accounted for by the 
combination of a statistical factor of four and the 
electrostatic effect of the negative charge of the 
acid-ion on the dissociation of the second hydro- 
gen, and presented an approximate mathe- 
matical discussion. In this paper, following a 
consideration of the Bjerrum treatment, an 
extension of the electrostatic theory of dissocia- 
tion constants will be developed. 

A thermodynamic comparison of the dissocia- 
tion constants K, and Ke of two organic acids 
HA, and HAz in the same solvent leads to the 
following relation 


RT log K,/K:=F°(HA 1) 
(1) 


where F°(HA,), F°(A2-), F°(HA2), and F°(A1-) 
are the reference values of the chemical po- 
tentials of the respective solute species taken at 
infinite dilution in the solvent. From a molecular 
point of view, it is readily seen that F°(HA;) 
+ F°(A,-) —F°(HA:») —F°(A,) is equal to N Aw’, 
where N is Avogadro’s number and Aw’ is the 
average reversible work expended in the trans- 
fer of a proton from the carboxylate group of an 
ion A, to the carboxylate group of an ion A;, 
situated at an infinite distance from A, in the 
given solvent. A part of Aw’ may be regarded as 


1 Bjerrum, Zeits. f. physik. Chemie 106, 219 (1923); see 
Greenspan, Chem. Rev. 12, 339 (1933). 


purely electrostatic in origin, that is, work done 
in moving the proton in the electric field pro- 
duced by the average charge distributions of the 
two ions. We may denote this part by Aw. An- 
other contribution arises from an entropy change 
associated with differences in molecular sym- 
metry, and equal to log [o(HA;)o(A2) 

and o(A2-) are the symmetry numbers of the 
respective particles. The symmetry number gives 
rise, in the ratio of the dissociation constants of 
a symmetrical dibasic acid, to the well-known 
statistical factor of four. Finally, and by no 
means to be overlooked, is the contribution to 
Aw’ associated with the intrinsic structure of the 
carboxyl group and carboxylate ion. However. 
if the two acids, HA, and HAzg, are, let us say. 
of the form RyCH2COOH and Rx,CH2COOH, and 
the possibility of structural resonance is absent. 
it is reasonable to suppose that the relatively 
distant substituent groups R; and R, will have 
little influence on the structure of the carboxy! 
group and carboxylate ion, except perhaps for 
an electrostatic polarization which may be taken 
into account in Aw, the electrostatic part of Aw’. 
Under these circumstances, the structural con- 
tributions to the work expended in detaching the 
proton from HA, will be canceled by a contribu- 
tion of equal magnitude gained by attaching the 
proton to A;. The net structural contribution 
to Aw’ is then zero, and we may write Eq. |!) 
in the form 


RT log NAw=2.303RTAPK 


Aw= W(HA,)+ W(As-) —W(HA2)—W(Ar) (2 
(HAs), 
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where W(ITA,), W(A2-), W(HA2) and W(A;) 
are the electrostatic free energies of the respec- 
tive molecules in the given solvent. The Bjerrum 
theory, in its most general form, is summarized 
by Eq. (2). Although Bjerrum developed his 
theory on the basis of a somewhat different 
argument, his fundamental assumption may be 
expressed in our terms, namely that the in- 
trinsic structure of the carboxyl group and car- 
boxylate ion is the same in the two acids. 
Although this assumption can certainly be no 
more: than a close approximation to the facts, 
that it is essentially correct is supported by 
the observation that the dissociation con- 
stants of two aliphatic acids, Ri}CH2COOH and 
R,CH2COOH generally differ but slightly when 
R, and R2 are nonpolar radicals. 

If we are interested in a comparison of the 
first and second dissociation constants of a sym- 
metrical dibasic acid, we are concerned with the 
molecular species 172A, HA~, A= and HA-. 
Since the symmetry number ratios, o(1/2A)/ 
o(IIA~) and o(A~)/e(JIA~) each have the value 
two, the statistical factor becomes four, and 
Eq. (2) becomes 


RT log Ki/4K2= NAw=2.303RTAPK (3) 
—2W(HA-). 


Bjerrum’s approximation is obtained if Aw is set 
equal to the electrical work e?/DR’ required to 
bring the two charges, —e, of the acid ions, HA- 
from infinity to their distance of separation, R’, 
in the di-ion A=, in a homogeneous medium 
having the solvent dielectric constant, D. This 
leads to results of fair accuracy? in the case of 
elongated straight chain dibasic acids, but with 
short molecules the values of R’ obtained are, 
in some cases, only a fifth of that expected on the 
basis of structural considerations. 

In the next section we shall describe a more 
refined calculation of Aw, which, although still 
somewhat schematic, avoids the treatment of 
the molecules and ions, in which the charges are 
imbedded, as structureless regions of dielectric 
constant D. Instead, we regard the molecules as 
forming cavities of low dielectric constant, Dj, in 
the solvent. Of course the introduction of an 
internal molecular dielectric constant D;, differ- 
ing from unity is somewhat artificial. Neverthe- 


* Gane and Ingold, J. Chem. Soc. 2153 (1931). 


less, its physical significance is clear, for it 
provides a crude means of allowing for the 
polarization of a molecule produced by the 
average charge distributions of polar or ionic 
substituent groups. A reasonable estimate of D; 
for aliphatic compounds is about 2.00, a value 
in the neighborhood of the dielectric constant 
of the liquid paraffin hydrocarbons. 

For certain purposes it is convenient to write 
Eq. (3) in the simple Bjerrum form, with the use 
of an effective dielectric constant, Dz, defined 
in terms of Aw, where the distance, R, is the 
separation of the protons in 


kT log K,/4K2=e?/DrR (4) 
Dr=e/RAw, 


rather than that of the negative charges in the 
di-ion A~. This mode of expression was intro- 
duced by Eucken* from empirical considerations. 

We shall also be concerned with a comparison 
of the dissociation constants of pairs of mono- 
carboxylic acids similar in structure except for 
the presence of additional polar groups in one 
of them. As an example of such a pair, we may 
mention acetic and chloracetic acids. The 
statistical factor, ¢, is evidently unity in these 
cases. If we adopt the notation of Eucken, using 
an effective dielectric constant, Dz, the ap- 
propriate form of Eq. (2) is 


kT log Ki/K2=eM cos {/DgR?=2.303kRTAPK (5) 
cos 


where M is the dipole moment of the polar 
group in the acid HA,, R is the distance from the 
proton of the point dipole, and ¢ its angle of 
inclination to the line joining it to the proton. 
We shall also make use of a more complicated 
form of Eq. (1) in which the charge distribution 
is not treated merely as a point dipole. Substitu- 
tion of D, the dielectric constant of the solvent, 
for Dz leads, as Eucken showed, to calculated 
values of R as much as six times smaller than 
might be expected. 

Perhaps the most serious approximation in our 
calculations of Aw which follow is the neglect 
of possible deviations of the local dielectric con- 
stant of the solvent in the neighborhood of a 
solute molecule from its macroscopic value, D. 
Such deviations might arise from electrical 


3 Eucken, Zeits. f. angew. Chemie 45, 203 (1932), 
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saturation, electrostriction, and other departures 
of the local density of the solvent from its bulk 
value. Ingold* has attempted to estimate cor- 
rections to the simple Bjerrum formula arising 
from electrical saturation and electrostriction. 
However, his calculations are based upon the 
Clausius-Mosotti equation, which is known to 
be invalid for polar liquids such as water. 
Further, even were the general method accurate 
the effect he calculated would be reduced to a 
fraction of the indicated magnitude if the solvent 
were excluded from the volume actually occupied 
by the molecule. It seems, then, entirely possible 
that he overestimates the importance of these 
effects. While the neglect of electrical saturation 
and electrostriction may lead to some inaccuracy, 
certainly the present treatment must be under- 
taken before finer adjustments are considered. 


II 


In order to calculate the electrostatic free 
energy of a spherical molecule of radius }, with an 
arbitrary charge distribution and immersed in a 
solvent of dielectric constant D, we make use of a 
formula developed by one of us in another con- 
nection.® A slightly different zero of energy is 
used here, since we have to consider changes of 
state involving the addition and removal of 
charges in the molecule. 


* (n+ 1)D-+nD,] 


3,1), 


k=1 l=1 


where the set of point charges e,:--e; charac- 
terizes the average distribution of charge in the 
molecule, 7; is the distance between the charges 
e,, and e:, r, and r,; the distances of e, and e; from 
the center of the sphere, and 3; the angle be- 
tween r, and r;. The functions P, (cos #1) are 
the ordinary Legendre functions, and D; is the 
internal dielectric constant of the molecule. If we 


‘Ingold, J. Chem. Soc. 2179 (1931). 
5 Kirkwood, J. Chem. Phys. 2, 351 (1934). 


consider the set of molecules HA;, HA», A, 
and A; of equal radii b, in which A, contains 
the charges e;---e, in excess of those contained 
in A», a simple calculation shows that 


1 1 (n+1)AG,, 


AG,= > dx), 


k=1 


where e¢ is the protonic charge, R; is the distance 
between e; and the proton in HA,,r and 7; are the 
lengths of the vectors joining the center of the 
sphere with the proton and the charge e,, re- 
spectively, and 3, the angle between r and 7;. We 
note that a similar formulation is possible when 
the excess of charge of A1~ over A» is continuously 
distributed with density p. If the ions are those 
involved in the first and second dissociation of a 
symmetrical dibasic acid, the charge distribution 
of the acid-ion HA~, which corresponds to Ar, 
contains a proton in excess of that of A=, which 
corresponds to Az~. If we denote by x the ratio 
(r/b)?, where r is the distance of each proton 
in H.A from the center of the spherical cavity of 
radius b, Eq. (7) becomes 


1 1 
Aw=—{— 
b (2x(1—cos #))! 
(n+1)x"P,(cos 3) 


1 1 
’ (8) 
+(; D; n+1+nD;/D 


where # is the angle between the lines joining the 
two protons of H2A to the center of the molecule. 
Unless we refer our calculation to the distance 
between the protons rather than the distance 
between negative charges in the ions, we neglect 
the effect of the carboxyl dipole. Although the 
sum in Eq. (8) may be expressed as a definite 
integral, an adequate approximation is obtained 
by expanding in powers of D;/D, and neglecting 
all terms in powers of D;/D higher than the 
first. Then, making use of the properties of the 
Legendre functions, the series involved may be 
summed. When the resulting expression is sub- 
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stituted in Eq. (3), we obtain 


é f 1 f 2 
logio K1/ 
2.303RRTLD D; 


2(2x(1—cos 8))! (2x(1—cos #))} 


(1—2x cos x 
(1—2x cos 3+x?)!+x—cos 


1= 


X log. 
1—cos 3 


(2x(1—cos #))} 
(1—2x cos 8+x?)! 
R=b(2x(1—cos 8))!, (9) 
where R is the distance between the protons in 
H.A. lf we make use of the notation of Eq. (4), 


in which an effective dielectric constant appears, 
we have 


1/De=fi/D+fo/Di. (10) 


In Fig. 1 the values of Dg are shown for all 
values.of ¥x from zero to one, and for three 
values of cos #, namely, —1, —1/2 and 0, 
corresponding to the cases in which the protons 
are on a diameter, placed at an angle, computed 


from the center, of 120°, and at right angles. D 


was taken as 78 for water at 25°; D; was assigned 
the value 2.00. The effective dielectric constant 
rises, with increase of Vx, from D; when the 
charges are in the center of the sphere and 
\ x=0 to values above the dielectric constant of 
water when the charges are on the edge of the 
sphere and ¥x=1. The fact that the effective 
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Fic. 1. Plot of Dg fora charged substituent against ¥ x. 
Ia curve I cos #=—1, in curve II cos 8=—1/2 and in 
curve III cos #=0. D; is 2.00 throughout. 
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Fic. 2. Curve A, plot of Dg against ¥x for a charged 
substituent with D; =2.50. Curve B, same except D;=2.00. 
Curve C, plot of Dg against ¥x for a dipolar substituent. 
In all cases cos 3 = —1. 


dielectric constant is greater, at the edge of the 
sphere, than the dielectric constant of the 
medium may at first glance be surprising. This 
need cause no concern, for Dz is only a dielectric 
constant by convention, being so defined that 
our equations take on the simple Bjerrum form. 
The result is in no way dependent upon the 
approximation made when Eq. (8) was ex- 
panded in powers of D;/D. In part IV of this 
paper, the sum in Eq. (8) is expressed as a def- 
inite integral, and the exact value of Dg for 
charges at the edge of a sphere and on a diameter 
obtained. 

In the computation of Dz, no arbitrary pa- 
rameters enter except the internal molecular 
dielectric constant, D;. However, it is assigned a 
reasonable value, common to all aliphatic mole- 
cules, and is not readjusted to meet the demands 
of individual cases. Further, the value of Dz 
depends primarily upon the existence of a low 
internal dielectric constant, and only secondarily 
upon the exact value assigned to it. In Fig. 2, 
curve A was computed for charges on a diameter, 
with D;=2.5, while curve B was computed for 
charges on a diameter, with D;=2.0. The differ- 
ence, while real, is small. 

Besides the effect of substituents bearing a 
total net charge, we wish to determine the elec- 
trostatic free energy difference, Aw, when HA, 
contains a polar group in excess of HA». We may 
regard the dipole as two equal and opposite 
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charges, so placed and of such a magnitude that 
the product of charge by the distance of separa- 
tion is equal to the dipole moment. Eq. (+8) 
gives the value of Aw in this case. We can also 
regard the charge distribution as a point dipole. 
If we restrict our attention to those cases in 
which the dipole and proton both lie equidistant 
from the center on a diameter, then 


1 1 4x} 8x? 
(1+)? 14x 


D; 


4 
ing (11) 
x} 


where all the symbols have their usual signifi- 
cance. In Fig. 2, curve C isa plot of Dz, obtained 
from Eq. (11), against ¥ x. 


III 


From the known dissociation constants of 
organic acids and Eqs. (4) and (5) an experi- 
mental value of Dz can be obtained for the cases 
in which HA, differs from HA, either by a 
charged or by a dipolar substituent. With 
charged substituents, Dg varies from 15 to 
about the dielectric constant of water, with 
dipolar substituents from about 3 to about 10; in 
both cases Dz is greater for greater values of r/b. 
For values of r/b between 0.5 and 0.9 which 
correspond to actual molecules, Figs. 1 and 2 
show that Dz, obtained by our theoretical 
treatment, varies in a similar manner. We note 
that the effective dielectric constant is much 
smaller for a dipolar than for a charged substit- 
uent at the same value of r/b.*° In this section we 
wish to present the computations for the appli- 
cation of our equations to a few simple cases. 

The exact location of the proton in the 
carboxyl group is still open to question. Available 
evidence’? seems to favor a double potential 
energy minimum for the proton adjacent to each 
of the oxygen atoms, at approximately the 
position of the proton in the normal O—H bond. 
To avoid undue complications in the calculations, 
we shall adopt the extended chain model, and 
conventionally place the proton on the extension 

6 Smallwood [J. Am. Chem. Soc. 54, 3048 (1932) ] has 
shown by means of a Born cycle that Dz will be unity if 
certain cancellations in the free energies of solvation, from 


the hypothetical gaseous state, of the acids and ions occur. 
7 See Price and Evans, Proc. Roy. Soc. A162, 110 (1937). 
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of the bond between the carboxyl carbon atom 
and its nearest carbon neighbor. Some con- 
vention, of which this is the most convenient is 
necessary solely for the purpose of estimating 3, 
the angle between the lines joining the protons to 
the center. 

The selection of a structural distance between 
the protons of a dicarboxylic acid for comparison 
with the value computed from the dissociation 
constant ratio is further complicated by present 
ignorance relative to internal rotations in the 
molecule. We obtain a probable upper limit for 
the structural distance by using the extended 
chain model with the proton conventionally 
located at a distance of 1.45A from the respective 
carboxyl carbon atoms on the extensions of the 
terminal carbon to carbon bonds. The distance 
of 1.45A is actually the greatest distance of the 
projection of the proton on the extension of the 
terminal carbon-carbon bond, when attached to 
either oxygen at the covalent distance 0.94A and 
with an angle of 90° between the C—O and O—H 
bonds. For the purpose of comparison we have 
also calculated the root mean square separation 
of the protons on the basis of free rotation. In 
this computation, we have employed the follow- 
ing formula for the mean square distance R,’ 
between the terminal atoms in a chain containing 
n bonds® 


n k-1 


n 
R,2= L2+2> II Nm, 


k=1 k=1s=1 m=s 


where 1m, m—1 is the cosine of the supplement of 
the angle between the m and m—1 bonds, and 
l, is the length of the k’th bond. We are aware of 
the fact that the free rotation value does not 
provide an absolute lower bound for the average 
structural distance. However, it is not unreason- 
able to suppose that the average structural 
separation of the protons lies between the 
extended chain and the free rotation value. In 
Table I we have presented the data for a few 
acids which are at least approximately spherical. 
It is seen that the distance computed from the 

8A special case of the formula developed by Eyring, 
Phys. Rev. 39, 746 (1932). It is worth while pointing out 
that for the above formula to be valid, free rotation is not 
actually necessary. The only requirement is that the mean 
values of cos gand sin ¢ vanish for adjacent bonds, ¢ being 
the angle specifying the projection of the i+1’th bond on 


a plane perpendicular to the 7’th bond. This may occur 
even if the rotation is highly inhibited. 
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TABLE I. Separation of protons in dicarboxylic acids.* 
Distance in A. 


FREE 
Rota- 
TION 


3.50 
4.12 


4.66 


Tuis 
PAPER 


3.85 
4.15 


4.80 


AcIp MAXIMUM 


4.44 
4.87 


6.66 


Oxalic acid 

Dimethyl malonic acid 

Tetramethyl succinic 
acid 


1.34 
0.96 


* The values in column 4 are given to the nearest 0.05A. Since the 
computed values of R in Table I refer to a single configuration of the 
molecule, they are not strictly suitable for comparison with mean 
distances calculated on the basis of free or inhibited internal rotation. 
Nevertheless, we have stated the maximum and free rotation distances 
in order to exhibit a possible range of variation of the interprotonic 
distance in the acids. With internal rotation, we should properly com- 
pute an average value of ApK over all internal configurations for 
comparison with the experimental value. Such a calculation is very 
tedious and does not seem worth while on the basis of free rotation, 
since in actual molecules the rotation is doubtless inhibited. We are 
unable to make the calculation for inhibited rotation, until future 
information about the hindering torques becomes available. 


dissociation constant ratio by means of our 
equation do actually fall between the suggested 
limits, while the values computed with the simple 
Bjerrum formula lie considerably below the free 
rotation value and are almost certainly too small 
on any basis. We remark that the Bjerrum 
formula is by no means as unsatisfactory for the 
longer chain acids. We have selected for dis- 
cussion a group of short chain acids the shapes of 
which roughly approximate spheres, since it is in 
these cases that the Bjerrum theory needs 
refinement. In a later paper, we shall develop a 
formula, on the basis of an ellipsoidal model, 
which passes asymptotically into the simple 
Bjerrum formula as the chain length increases. A 
knowledge of the torques hindering internal 


rotation would permit a more exact application 


of our formulas involving a calculation of the 
average value of ApK over all internal configu- 
rations of the acids and ions. 

For these computations it was necessary to 
estimate the radii of the spheres. In the absence 
of specific information, we have fallen back on 
Traube’s rule. The partial molar volume of 
dimethyl malonic acid in water, then, is about 
99.1 cc, and the diameter of the molecule, con- 
sidered as a sphere, is 6.77A. Further, since we 
have decided, conventionally, to place the 
protons on the extensions of the terminal carbon 
to carbon bonds, we may assign the angle 
cos = — 0.333 to the lines joining the protons to 
the center. From Eq. (4) we obtain 


e 
Vv xD E> ’ 
2.303bkTApK (2(1—cos #))! 
r/b=NV x. 
* Traube, Saml. chem. chem-tech.. Vortr. 4, 255 (1899). 
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Substituting the value of 2.29 for ApK for 
dimethyl malonic acid found by Gane and 
Ingold? in this expression, we find ¥ xDz equal 
to 19.1. Values of ¥ xDz computed with the aid of 
Eq. (10) are given in Table II for varying Vv x and 
cos 3. Using this table, the value of ¥ x can be 
obtained in any particular case by linear 
interpolation. The computations for the other 
acids were carried out in a similar fashion, again 
using the data of Gane and Ingold. The angle 
between the protons for oxalic and tetramethyl 
succinic acids was chosen as 180°. 

As an example of the effect of a simple dipolar 
substituent we present the computations for the 
ratio of the dissociation constant of acetic to 
chloracetic acids. Since the molecule does not 
approach a true sphere, we arbitrarily placed the 
center at the midpoint of the line joining the 
proton with the midpoint of the carbon-chlorine 
bond. We employed the value 1.56 e.s.u. 
for the dipole moment of the C—Cl bond;!° the 
values of the dissociation constants of acetic 
and chloracetic acids were taken from Landholt- 
Bornstein." The data for chloracetic acid are 
given in Table III. Other calculations have shown 
that, in chloracetic acid, the acidity is relatively 
little effected by multipole moments of higher 
order than the second. 

Finally we may mention that, in common with 
the previous electrostatic theories of dissociation, 
the temperature coefficients experimentally de- 
termined” do not agree with those calculated, on 
the basis of a constant value of R. Although a 
real variation of the average value of R with 


TABLE II. Values offcDe(x). 


SESS 


10Smyth, Am. Chem. 
Structure, p. 209 (1937). 
1 Landholt-Boérnstein, Physikalische-Chemische Tabellen, 


Eg IIIc. 
12 Jones and Soper, J. Chem. Soc. 133 (1936). 
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Vx cosd=—1 cosd=—} cos 3=0 

il 0.00 0.00 0.00 0.00 
0.10 0.25 0.24 0.23 

e .20 4 0.60 0.55 

n .30 1.16 1.02 

40 2.11 1.76 

50 3.82 3.02 

l. .60 7.15 5.42 
.70 14.2 10.6 

80 30.9 23.3 

90 68.9 56.4 

t 95 96.3 83.2 

- 1.00 116.3 103.5 

n 
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TABLE III. and dipole in chloracetic 
ac 


BJERRUM 


0.55 


Minimum THIS PAPER 


1.84 2.90 


MAXIMUM FREE ROTATION 


3.39 2.99 


temperature is not unlikely, it may be that a more 
refined theory is necessary before accurate esti- 
mates of temperature coefficients can be made. 
In any case, the temperature coefficients are not 
large. 


IV 
In this section we shall present some of the 


mathematical details omitted earlier in the paper. 
Eq. (7) may be written in the following form. 


Aw=e> 


k=1 


H,= sian, cos (12) 


(n+1)x"P,(a) 
n=0 1+n(1+w) 


where x; is equal to rr;/b? and w is the ratio, 
D;/D. From the properties of the Legendre 
functions, we may write 


S(x, a) = 


© 1 
So(t, a) = P,(a) = 3 
(t, a) (a) (13) 


It is readily verified that S(x, a) may be ex- 
pressed in the form 


S(x, a) f 
0. 


Simplification of Eq. (14) leads to 


(1+w) 


So(t!*, coat] (14) 


1 
S(x, a) = 


{Su a) a) 
w 1+w 


z 
(+e) 


Si(x, a) Ate) 


G. KIRKWOOD AND F. H. WESTHEIMER 


menoenenage in powers of w with the neglect of 
terms in w gives 


S(x, a) = So(x, a) +o[.S1°(x, a) — So(x, a) 
(1 (16) 


1 
S)°(x, a) =— log 
x 


Evaluation of S;(x, a) by graphical quadrature 
shows that the expansion (16) is adequate in all 
cases of interest although Dg calculated from (10) 
is 2 percent high at x = 1. Substitution in Eq. (12) 
gives 


1 
Db 


——log 


1 — 2x; cos 3, +2, —Ccos 


1—cos 3; 


Eq. (9) for the symmetrical dibasic acid follows 
at once from Eqs. (12) and (17), only one term 
occurring in the sum of Eq. (12). For a dipolar 
distribution consisting of two point charges +e, 
and —e, Eq. (12) takes the form 


Aw=e_H,—H_], 


where H, and H_ are to be computed from Eq. 
(17) for the respective charges. We note that H; 
may be regarded as a function H(r,, r) where t 
and r are the vectors extending from the center 
of HA, to the charge k and the carboxyl proton, 
respectively. Thus for a point dipole, we may 
write 


(18) 


Aw=eM-V,/1(n, r), (19) 


where M is the vector moment and JV, the gradi- 
ent operator associated with the position r; of the 
dipole. The general expression is rather compli- 
cated and will not be given. However, for a 
molecule in which the proton and dipole are 
located on a diameter, Eqs. (17) and (19) yield 
Eq. (11) of Section IT. 
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The Electrostatic Influence of Substituents on the eaneuemees Constants of 
Organic Acids. II. 
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(Received June 4, 1938) 


The theory of the influence of substituents on dissociation constants presented in an earlier 
article of the same title is extended to acids ellipsoidal in shape. Bjerrum’s formula is obtained 
in the limiting case of unit eccentricity of the elliptical section of a prolate ellipsoid of revolution. 
The influence of dipolar as well as ionic substituents can satisfactorily be treated. 


I 


T has been proposed by Bjerrum! that the 

influence of a polar substituent on the dis- 
sociation constant of an organic acid is primarily 
electrostatic in origin. On the basis of this 
hypothesis, the ratio of the dissociation con- 
stants of two similarly constituted acids may be 
calculated from electrostatic work necessary 
to transfer a proton from the first acid to the 
ion of the second. Because of the fruitfulness of 
Bjerrum’s idea, we considered it desirable to 
introduce certain refinements into the calcula- 
tion of the work of transfer. In a previous 
article, we showed? that the electrostatic work 
is dependent not only upon the charges present 
and the distance between them, but also upon 
the shape of the molecule and upon the positions 
of the charges. We have already presented the 
equations for a spherical model, and have ap- 
plied them to several aliphatic acids which 
approximate a sphere in shape. In the present 
paper, similar considerations lead to the equa- 
tions for an ellipisod of revolution, a model well 
fitted to treat the longer aliphatic dibasic acids. 

We are concerned with four ellipsoids, repre- 
senting the molecular species HA, Ai~, 
and A,-. For the sake of simplicity, we have 
assumed that the protons of the carboxyl group 
and the dipoles in the molecules lie on the line 
connecting the foci of the ellipsoid, and, in the 
actual computations, have further assumed that 
the charges are actually at the foci, although 
this latter limitation is not inherent in the equa- 
tions developed. In the present discussion, it is 
convenient to use confocal elliptical coordinates. 


: Bjerrum, Zeits. f. physik. Chemie 106, 219 ne. 
(1938 


en and Westheimer, J. Chem. Phys. 6, 506 


If 7; and re are the distances of a point from the 
two foci, and R the distance between the foci, 
then and ¢ is the 
angle by which the point has been rotated past a 
fixed plane through the major axis. Since we 
have chosen an ellipsoid of revolution with all 
charges on the major axis, the model is axially 
symmetrical and the angle ¢ does not enter the 
formulation. 

Following the same reasoning advanced for 
the spherical model, we approximate the total 
free energy increment in the reaction 7A,+A2~ 
@HA.+A; by the sum of the electrostatic 
contribution, Aw, and the symmetry contribu- 
tion RT log o, where a is the statistical factor, 
equal to four for symmetrical dibasic acids. 


kT log Ki/0K2=2.303 RTAPK =Aw 
AW=W(HA))+W(Az) 
—W(HA2)—W(Ar), 


where W(HA,), W(A2-), W(HA2) and W(A1) 
are the electrostatic free energies of the respective 
molecular species. In order to determine the 
electrostatic work, Aw, it is first necessary to 
solve Laplace’s equation in confocal elliptical 
coordinates. Then, from the electrostatic po- 
tential of the charges, the electrostatic work 


(1) 


TABLE I. Effective dielectric constant for varying do. 
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can be computed. The actual mathematical 
treatment is presented in Part III of this paper. 
In that section it will be shown that, if HA, 
differs from HA, by a protonic charge, then 


e 
2.303kTR 
n=01+(D;/D)C, 
Qn (Xo) 
U,= 2 +d 
(2n ) 


/Pn(Xo) 


Further, if HA, differs from HA: by a point 
dipole of moment M inclined at an angle ¢ to the 
major axis of the ellipsoid, then 


eM cos¢ 
2.303kT R? 


wi 


Qn (Xo) 
P n(Xo) 


In these formulas, P, represents a Legendre 
function of the first kind, D is the dielectric 


ApK = 


Y,=n(n+1)(2n+1) 


(—1)". 


re 


Fic. 1. Curve A is the plot of Dg for a charged substit- 
uent against Ao, curve B is the plot for a dipolar substit- 
uent. The slope of the curves at Ao=1 has not been 
determined. 
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3 Eucken, Zeits. f. angew. Chemie 45, 203 (1932). 
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constant of the solvent, D; the internal molecular 
dielectric constant? and \=Xpo is the equation of 
the ellipsoidal cavity in the solvent created by 
the molecule. 

Q, represents a Legendre function of the 
second kind ; Qy and Q; are given by the equations 


A+1 
Qo(d) =3 log. 


log, ——— 


and the rest can be obtained from the first two 
by means of the recursion formula, 


(M+ 1) Qn+i(d) = (20+ 1)Qn(A) 


As in the previous article, an effective dielec- 
tric constant, Dg, is defined in terms of the 
Bjerrum equation.': * For the case in which //A, 
has a charge in excess of HA2, 


ApK =e2/2.303DeRkT (4) 


and in the case in which HA, contains a dipole in 
excess of HA» 


ApK=eM cos ¢/2.303DeRkT. 


In Table I the values of Dz accurate to 1 percent 
for molecules in water at 25° computed from 
Eqs. (2), (3), (4) and (5) for both dipolar and 
charged substituents are recorded for typical 
values of Xo. Di, as previously, has been set 
equal to 2.00. 

In Fig. 1, Dz has been plotted against Xo. At 
large values of Xo the charges are brought close 
together in the medium of low dielectric constant, 
and Deg approaches the internal molecular 
dielectric constant D;. On the other hand, as \o 
decreases toward 1, the ellipsoid thins out into 
a rod, and at A» equal to 1 disappears altogether, 
leaving the charges immersed only in the solvent. 
As in the spherical case, the values of the effec- 
tive dielectric constant depends primarily upon 
the existence of a low internal dielectric constant, 
only secondarily upon its exact value. 

With the ellipsoid as with the sphere, values of 
De greater than the dielectric constant of the 
solvent are obtained. With the ellipsoid, how- 
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ever, since the effective dielectric constant must 
approach D as the ellipsoid vanishes, the values 
of Dz go through a maximum. 

Here again our theoretical treatment is in 
agreement with the experimental fact that, for 
the same value of Xo the effective dielectric 
constant is much greater in the case in which the 
substituent is charged than in the case in which 
the substituent contains a dipole. 


II 


The equations presented in Part I will now be 
applied to some actual molecules which ap- 
proximate an ellipsoid of revolution. We have 
compared the first and second dissociation con- 
stants of three dicarboxylic acids, succinic, 
glutaric and suberic, and have compared the 
dissociation constants of chloracetic acid and of 
8-chlorporpionic acid with the dissociation 
constants of the corresponding unsubstituted 
acids. Chloracetic acid, intermediate in shape 
between a sphere and an ellipsoid, has been 
treated by both methods. In order to fix the 
boundries of the molecules, we have again 
estimated the partial molar volume in aqueous 
solution from Traube’s rule, and have made 
use of the following equation relating Xo, R, 
and the molecular volume, V, of an ellipsoid of 
revolution : 


Ao? — Ao = (6/7) (V/R'). (6) 


This equation must be solved simultaneously 
with Eqs. (4) and (5) and with Fig. 1 to obtain 
the value of R in particular cases. We note that 


(7) 


Dz (charge) e? 
—-2.303RTAPK \6V 


and 


De (dipole) eM cos¢ ( 7 ) 


We have prepared a table of Dg (charge)/ 
and Dg (dipole) /(Ao?— Xo)! for various 
values of Xo and can, then, estimate the value of 
\o for any particular case by linear interpolation. 
From this value of \» and Eq. (6), R is easily 
obtained. These values are presented in Table II. 
The values of R, computed for the acids selected, 


‘Traube, Saml. chem. chem-tech., Vortr. 4, 255 (1899). 


ACIDS 


TABLE II. 


Dg (CHARGE) /(Ao?—Xo)! (DIPOLE) /(A0? —do)3 


> 


PAPER* 


MAXI- 


AcID MUM BJERRUM 


3.99 
6.53 
10.9 
0.55 
1.06 


6.66 
7.39 
11.46 
3.39 
4.61 


Succinic acid 

Glutaric acid 

Suberic acid 

Chloracetic acid 

8-Chlorpropionic 
acid 


* The values of R in this column are given to the nearest 0.05A. For 
the chloracids, ¢ was computed as the angle between the carbon- 
chlorine bond and the line joining the midpoint of this bond with the 
point of location of the proton on the extension of the terminal carbon- 
carbon bond into the carboxyl group. The values of cos ¢ are 0.59 and 
0.94 for chloracetic and 8-chlorpropionic, respectively. The carbon- 
chlorine bond moment was taken as 1.56 Debye units. 


with the aid of Table II, are presented in Table 
Ill. 

Since the computed values of R in Table III 
refer to a single configuration of the molecule, 
they are not strictly suitable for comparison with 
the mean distances calculated on the basis of 
free or inhibited rotations. Here, as in the pre- 
vious paper,” we have stated the maximum and 
the free rotation distances in order to exhibit a 
possible range of variation of the interprotonic 
distance. The values of ApK were taken from 
Landholt-Bérnstein.® 

It is at once seen that both the Bjerrum for- 
mula and the present treatment give reasonable 
values for the long aliphatic dibasic acids. The 
advantage of the present treatment over the 
simple Bjerrum approximation is most clearly 
exemplified in the cases of the chloro-acids. 

The approximations and possible sources of 
error involved in this formulation have been 
discussed in the previous paper.? They include 
neglect of electrical saturation and electro- 
striction, the use of an internal molecular 
dielectric constant, the use of Traube’s rule, and 


5 Landholt-Bornstein, Physikalische-Chemische Tabellen, 
Erg. IIIc. 
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the difficulty in obtaining a structural value of R 
with which to compare our value, due to present 
ignorance concerning rotation in polyatomic 
molecules in aqueous solution. We are further 
limited by the fact that we have placed the 
charges or dipoles at the foci of the ellipsoid. 


Ill 


In the third part of this paper, we wish to 
derive the equations presented in Part I. 
Laplace’s equation, in confocal elliptical co- 
ordinates for an axially symmetric poteutial 


y, is 
0 0 
On? Or Ou? Ou 


This equation must be satisfied both inside and 
outside the ellipsoid, \o, forming the boundary 
between the molecule of dielectric constant D; 
and the solvent of dielectric constant D. The 
potential ¥; in the interior of the molecule, 
satisfying Eq. (7) may be written as follows 


D;k=1 n=0 


where the B, are constants, ¥; having no singu- 
larities except at the positions of the point 
charges ¢,:--@2, characterizing the average 
charge distribution of the molecule. Since the 
potential must vanish at an infinite distance from 
the molecule, the solution, y., of Eq. (7) outside 
the boundary ellipsoid may be expressed in the 
form 


(9) 


where the P,, and Q, are the Legendre functions 
of the first and second kinds, respectively. The 
boundary conditions namely, continuity of the 
potential and of the normal component of the 
dielectric displacement everywhere on the bound- 
ary ellipsoid, Ao, require the following relations 
to be satisfied. 


Wi(Ao, =e(Ao, #) 


WESTHEIMER AND J. 
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introduce into Eq. (8) the Neumann expansion 
for each of the terms 1/|r—7;| 


ek 


2 
(20+ 1)BnPn(u)Qn(r) 


k=1 lr—r;.| R n=0 


LeP n (ux) 


(11) 


valid when each charge e; is situated at a point 
wu, on the ellipse, \=1, which is the line joining 
the foci. R is the distance between the foci. 
Substitution of (8) and (9) into Eq. (10) and 
use of the property of orthogonality of the 
P,(u) in the interval, —1=y~=+1, yields the 
following set of linear equations for the determi- 
nation of the A, and B, 


(2/R)(2n+1)(Bn/Di)Qn(Xo) 
+B,P,(Xo) =A,Q,(Xo) (12) 
n' (Xo) 
+ (2/R) (2n+ 1)BnQn’ (Xo) =DA nQn’ (Xo) . 
Solution gives the following expression for 
the B,, 


D D; 


x Qn(Ao) Qn’ (Ao) 
Qn! (Xo) — (Di/D) Pn’ (Ao) Qn (Ao) 


where P,,/(Ao) and Q,’(Xo) are the first deriva- 
tives of the Legendre functions. The derivatives 
may be eliminated with the aid of the recursion 
formula, 


(A2—1)P,,’(A) =n[AP, (A) — 
Qn (Xo) 


P,,(Xo) (14) 


Pr_i(do)/Pr 
— Qn—1(Xo)/Qn(Ao) 

The electrostatic work of charging the molecule 

is obtained by substitution of Eqs. (8) and (14) 

into the general formula 


(13) 


W= Mk)s (15) 


k=1 


Before applying the boundary condition we 
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all charges being located on the line \=1. We 
are interested in calculating AW for a set of 
molecules Ay and 


—W(HA2)—W(A,-). (16) 


Use of Eqs. (8), (14) and (15) under the assump- 
tion that all species have the same size and 
shape yields 

s & 


sW=— — 
7% 


Xo — On—1(Ao) /Qn( Ao) 


a, = 
k=1 


where the sums extend over the set of charges 
é\:--é2 which the acid HA, contains in excess of 
those of HA; and 7, is the distance of charge 
e, in HA» from the proton, e, located at the focus 
(1, 1) of the ellipsoid. In the symmetrical dibasic 


acid, the set of charges e;: - -é2 consists of a single 
proton, e, located at the other focus, (1,—1) and 
Eq. (17) becomes, since P,(—1)=(—1)", 


2(2n-+1)(—1)" 


n=0 


P (Xo) DI 


Qn—1(Xo) /Qn(Ao) 
Eqs. (1) and (18) lead to Eq. (2). When HA,» 
contains a dipole moment located at the focus 
(1,—1), in addition to the charges constituting 
HA, we find 


a, =2eM cos ¢ 


li (= P,( (19) 
im 
R 


M cos 4 


where M cos ¢ is the component of the dipole 
moment along the axis of the ellipsoid. By sym- 
metry it is readily seen that its component 
perpendicular to the axis contributes nothing. 
Eqs. (1), (17) and (19) lead to Eq. (3). 


we 
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Photochemical Studies. XXVII. The Effect of Radiation on Mixtures of Nitrogen 
Dioxide and Nitrous Oxide and its Relationship to the Photochemical 
Decomposition of Nitrous Oxide 


F. C. HENRIQUES, Jr., A. B. F. DuNcAN AND W. ALBERT Noyes, JR. 
Metcalf Chemical Laboratory, Brown University, Providence, R. I. 


(Received June 14, 1938) 


Oxygen atoms were produced in NO2—N,O mixtures by the photochemical decomposition 
of NOz alone. No reaction of these atoms with N.O could be detected. Results of calculations 
give a minimum activation energy of 14,000 cal. for the reaction O+ N2O=2 NO and a probable 
maximum activation energy of 24,000 cal. The activation energy of the reaction O+NO, 
= NO-+Oz: is estimated to be 4000 cal. It is concluded that ‘wo primary processes are necessary 
to explain the direct photochemical decomposition of NxO+hvy=N2+O and 
=NO+N. The latter reaction must be followed by N+N:O=NO+N; to explain previously 
determined quantum yields. The energy of activation of this reaction is probably less than 


10,000 cal. 


N a previous article! the number of molecules 
formed upon irradiation of nitrous oxide per 
absorbed quantum, and which were uncondensed 
by liquid air, was determined. It was found that 
the over-all reaction could be represented 
empirically by the equation 


4N,0 =2.8N2+2.5NO+0.702, (1) 


which agrees roughly with results previously 
obtained by Macdonald.” It had also been es- 
tablished definitely that nitric oxide is one of the 
products by Wulf and Melvin* who showed that 
the bands of this substance appear in absorption 
when nitrous oxide is irradiated. Since nitrogen 
dioxide appears after condensation with liquid- 
air oxygen, also, must be formed as shown by 
Macdonald.? 

The value of the quantum yield for nitrous 
oxide decomposition will depend upon the mecha- 
nism assumed for the primary process and for 
the secondary reactions. Molecular oxygen can 
only be formed if at least part of the primary 
process is 


N,O+hy=N2+0. (2) 


Since nitric oxide is formed, if (2) is the sole 
primary process a secondary reaction 


O+N,0=2NO (3) 


1W. A. Noyes, Jr., J. Chem. Phys. 5, 807 (1937). 

2 J. Y. Macdonald, J. Chem. Soc. 1 (1928). 

30. R. Wulf and E. H. Melvin, Phys. Rev. 39, 180 
9 ‘);L. Henry, A.-M., Compt. rendus 200, 656 (1935). 


(1 


must occur. An accompanying reaction 
0+0(+M) =0.(+ M) (4) 


must occur either on the walls or as the result of 
termolecular collisions to explain the known 
production of Ox». 

It was suggested! that an alternative expla- 
nation could be based on two primary processes, 
namely (2) and (5) 


N,O+hvy=NO+N. (5) 
The fate of the nitrogen atoms might be either 
(6) 
or N+N(+M)=N2(+) (7) 


or both. It should be re-emphasized that (5) 
cannot be the sole primary process unless a 
reaction between nitrogen atoms and NO is used 
to explain the appearance of Oz. The latter is 
improbable in view of the small amount of 
reaction which was allowed to occur and the fact 
that the extent of the reaction did not affect in 
any marked manner the value of the yield. 
The present experiments were undertaken 
with a view to obtaining evidence for or against 
reaction (3), thus permitting a decision to be 
made as to whether primary process (2) alone or 
both primary processes occur. It is true that 
Musgrave and Hinshelwood‘ have postulated 


reaction (3) to explain the appearance of NO 


4 Musgrave and C. N. Hinshelwood, Proc. Roy. Soc. 
(London) A135, 23 (1932). 
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during the thermal decomposition of nitrous 
oxide, but their experiments were carried out at 
about 1000°K so that this reaction might have an 
energy of activation too high to permit it to be of 
importance at room temperature. 

A convenient source of oxygen atoms is found 
in the primary dissociation of nitrogen dioxide by 
ultraviolet light, since this substance absorbs in a 
region where nitrous oxide is transparent. Henri® 
reports that nitrogen dioxide possesses two 
regions of predissociation, one near 3700A and 
the other below 2450A leading to the production 
of *P and 'D oxygen atoms, respectively. The 
photochemical decomposition of nitrogen dioxide 
has been investigated by Norrish® and by 
Dickinson and Baxter.” Approximately two mole- 
cules of nitrogen dioxide decompose per quantum 
absorbed, a fact which may be explained by the 
reactions 


NO.+/v=NO-+0, 


(8) 
(9) 


EXPERIMENTAL METHOD 


Nitrogen dioxide was obtained from pure, dry 
lead nitrate and purified by fractional distillation 
with liquid air through anhydrous magnesium 
perchlorate. The middle fraction was vaporized 
into a series of small tubes with capillary ends, 
which were sealed off. The nitrous oxide was a 
commercial product and was purified by bubbling 
slowly through saturated potassium hydroxide. 
It was fractionally distilled several times through 
anhydrous magnesium perchlorate with liquid 
air and the final middle fraction was collected in a 
bulb provided with a stopcock. 

In order to observe possible separate effects 
due to *P and 'D oxygen atoms two light sources 
and reaction vessels were used. The normal *P 
atoms were produced from nitrogen dioxide by 
using a Pyrex tube 100 cm long with plane 
Pyrex windows sealed on and as a light source a 
General Electric type H3 mercury lamp. The 
Pyrex windows did not transmit radiation below 
3000A. A fused silica vessel 20 cm long with 
sealed on plane windows was used for the 'D 

iY. Henri, Nature 125, 202 (1930). 

R. G. W. Norrish, J. Chem. Soc. 761 (1927); 1158, 
1611 (1929), 


. G. Dickinson and W. P. Baxter, J. Am. Chem. Soc. 
50, 744 (1928). 
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atoms. This was illuminated with a cadmium 
spark operated at 2 kva, 25,000 volts. Radiation 
between 2150 and 2350A was isolated by filters. 
Two cm of 0.05 m acetic acid containing 1 
percent acetone absorbed from 2350 to 2900A 
and below 2150A. Five centimeters of chlorine 
gas at atmospheric pressure absorbed from 2800 
to 4300A, above which the nitrogen dioxide is not 
dissociated.* 7 The acetic acid-acetone solution 
was renewed every three hours. The reaction 
cells were attached as shown in Fig. 1. 
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Pressures were measured with a glass Bourdon 
gauge which operated a lever carrying a small 
mirror. Images of a graduated scale placed at a 
distance of fifteen feet were read with a telescope. 
Readings were reproducible to 0.3 cm on the 
scale. Temperature variation over a period of 
time decreased this to perhaps 0.5 cm. Three 
gauges, designated A, B and C, were used during 
the course of the work. Calibrations showed that 
a 1 cm scale reading was equivalent to 0.00460, 
0.00719 and 0.00201 cm of mercury, respectively, 
for the three gauges. 

In Fig. 1, S indicates stopcocks lubricated with 
low pressure grease, T indicates mercury cut-offs 
and A and B constrictions for sealing off. 
Projecting through S2 is a capillary tube sealed 
at the upper end whose use will be apparent. 

To explain the procedure used to obtain the 
results it is necessary to refer to reactions (3), (8) 
and (9) as well as the sum of (8) and (9) 


2NO,=2NO0+0;, 
and the sum of (3) and (8) 
NO.+N,0=3NO. (5’) 


After illumination there will be a pressure 
increase AP,+AP; due to reactions (4’) and (5’). 
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It is known that the termolecular reaction 
2NO+0,=2NOz2 goes to completion at liquid- 
air temperatures. This was verified in our experi- 
ments. Hence the increase in pressure after 
illumination is AP,+AP,; and the residual pres- 
sure after freezing down a number of times is AP3. 

Bourdon gauge readings are called R1, R2, etc. 
Some of these were made merely to check 
constancy of conditions and do not appear in 
later calculations. No readings were taken until 
the temperature difference between the gauge 
and reaction vessel was approximately 0.2-0.3° 
which seemed to be the value most easy to 
maintain. 

The entire system was pumped and baked out 
until the pressure was less than 10-5 mm and the 
reading R1 made. S4 and 73 were closed and 
N,O introduced through S5 to a pressure of 
10-20 cm and the pressure read on 73 after 
which S5 was closed. The N;O was then frozen 
out by placing liquid air on F and the system 
pumped through S4. T3 and S4 were closed, the 
N.O evaporated at —80°C and the pressure was 
read on 73. The reading R2 was made, S1 was 
closed and the apparatus sealed off at A and the 
reading R3 made. 

NO, was introduced by breaking the capillary 
end of a tube with the magnetic hammer, H, and 
R4 made. Side tubes attached to the reaction 
vessel and gauge jacket were simultaneously 
cooled with liquid air (to avoid breaking the 
gauge) and R5 taken. The liquid air was removed 
and one hour allowed for diffusion, after which S1 
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was closed and R6 made. The N20 — NO» mixture 
was next illuminated for a period of four to nine 
hours, after which R7 was made. 

After illumination the gases were continually 
distilled between two wide side tubes attached to 
the reaction vessel. At about two hour intervals 
the liquid air was removed and the gases evapo- 
rated. This was continued until readings R§ 
(with liquid air) and R9 (without liquid air) 
became constant. S2 was then turned, breaking 
the capillary, and S1 and 71 opened to equalize 
the pressure. R10 was made and in every case 
agreed with R1, R2 and R5 within experimental 
error. In experiments with 'D atoms, the system 
was sealed off at B before R6 was made. 

For experiments at higher temperatures, the 
reaction vessel was covered first with wire gauze, 
then asbestos and wrapped with resistance wire. 
The vessel was heated between readings R6 
and R7. 

RESULTS 


The necessary pressures, P(NOz), AP; (final 
pressure increase while frozen down), APs, AP; 
are obtained from gauge readings as follows: 
(a) P(NO:2) = R4—R3; (b) APs +AP;=R7— Ro: 
(c) AP;=R9—R6; (d) AP;=R8—RS5. During 
distillation of the gases by liquid air after 
illumination, the pressure decreased steadily for 
two hours and then remained constant for three 
to five hours. This reading is called AP°. There- 
after the pressure decreased again and finally 
reached AP;. Since the volume of the fused silica 
cell was much less than that of the Pyrex cell, the 


TABLE I. Summary of experimental results. (Pressures are in centimeters of Hg). 
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time to bring AP,+AP, to AP; was only three to 
five hours in the experiments with 'D atoms. The 
relative rapidity of the reverse reaction 2NO+O, 
=2NO; prevented the observation of a AP° here. 
In the experiments with *P atoms, AP® may 
possibly indicate a solid complex, resulting froma 
system containing NxO, NOz, NO and O: which 
has an appreciable vapor pressure at liquid-air 
temperatures. 

The results are recorded in Table I. Experi- 
ments 1A—-5 were made with *P atoms and 
6A-9B with 'D atoms. 

(1) In this run no distillations were made. The 
AP; obtained was probably AP°. 

(2) Although the reaction cell was flamed 
before the NO was introduced, the pressure 
increase was probably due to gas coming off the 
walls. The facts that in 7A there is no difference 
between APs+AP; and AP; and that the 
photochemical decomposition of NO leads to 
NO and Oy verify this belief. 

(3) AP; may be due to that fact that some 
mercury was present. This was avoided in 
subsequent runs by protecting the reaction 
system with a dry-ice-ether mixture during 
evacuation. 

At the nitrogen dioxide pressures used, less 
then four percent of this substance was present as 
N,O; at room temperature and less than one 
percent at the higher temperatures. 


DISCUSSION 


Following the primary process (8), there are 
two competing reactions for oxygen atoms 


(rate =kg(O)(NOs2)), (9) 


0+N,0=2NO (or Ne2+Oz2) 
(rate=k3(O)(N2O)). (3) 


Within experimental error all of the oxygen 
atoms reacted with NOs. Although this does not 
exclude the possibility of reaction (3), 9 must be 
very much greater than k;. It is possible, from the 
experimental data, to fix a lower limit to the 
activation energy of (3). From the previous 
resultst on the photochemical decomposition of 
nitrous oxide (assuming reactions (2), (3) and 
(4)) it is possible to calculate an approximate 
upper limit to the activation energy of (3). 


Assumption of a steady state gives 
d(NO:)/d(N:0) 
= (2ko(NOz) (10) 
=(K(NO:2)/(N20) +1). 
This is of the form dy/dx—Ky/x=1. 
Integration gives 
(NO2)(N2O)-* (11) 


In run 5, which will give the highest value of 
the lower limit, P(N.O) (initial) is 18.5 cm and 
P(NO.,) (initial) is 0.216 cm. AP, is assumed to be 
0.0036 cm. The final P(N.O) is 18.5—AP; 
=18.496 cm, the final P(NO,.)=0.216—2AP; 
—AP,=0.040. (AP;+AP;=0.090 cm, T = 265°C.) 

Therefore 


(0.216) (18.5)-* — (0.040) (18.496)-* 
(18.5) !-K — (18.496) !-« 
(1—K) 


which reduces to 


0.216 —0.040 x (1.000195)* 
= (18.5 — 18.4964 (1.000195)*) /(1—K). (13) 


Remembering that (1+A)* =1+KaA if A is 
small, one finds 


K 22.2X10! and ky ks 21.1104. 


Now k=fk® exp (—E/RT) for a bimolecu- 
lar reaction where f is a steric factor and 
(o is the mean of 
the molecular diameters and the M's are mo- 
lecular weights). Certainly the mean molecular 
diameters and masses are almost the same and 
therefore k3°=k,°. The steric factors should be 
about 0.01 and there is no reason for believing 
that they are different for the two reactions. 

Spealman and Rodebush® estimate that ky 
= 10-5 at 40°C. From their results we calcu- 
late the energy of activation of (9) to be 4300 cal. 
and hence that of (3) to be 14,000 cal. 


8 Schumacher (J. Am. Chem. Soc. 52, 2377 (1930)) il- 
luminated mixtures of nitrogen dioxide and hydrogen and 
assumed from his results that there were two competing 
reactions (a) O+H.2=OH+H; (b) O+NO.=NO+0,. 
In order to account for the fact that all of the oxygen 
atoms seemed to react with the nitrogen dioxide, he states 
that (a) must be at least 10* times as fast as (b). It is 
interesting to compare this result with our value of kg/ks. 

9M. L. Spealman, and W. H. Rodebush, J. Am. Chem. 
Soc. 57, 1474 (1935). 
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We may now calculate the maximum energy of 
activation of (3) from the previous results! using 
reactions (2), (3), (4) and (15) 


O+ walls = 302 (rate=;(O)/2(N2O). (15) 


Assuming a steady state and integrating we 
obtain 


(Oz) (16) 
(NO) = 2k3(O)(N20)¢. (17) 


Assuming the usual formula for k3, we obtain 
from (16) and (17) 


4(O2)(N20) /(NO)? 
= (ks exp (2E3/RT))/t(fs)?(ks°)? 
7 (Ris exp (E3/RT))/fak3°(NO)(N2O). (18) 


It is possible to evaluate ki;/k3;° by Smolu- 
chowsky’s equation Z = 37X?/4L? where Z is the 
number of collisions suffered by the oxygen 
atoms with the nitrous oxide molecules before 
they reach the wall, X is the mean displacement 
of the atoms in reaching the wall and L is the 
mean free path of the oxygen atoms. 

With a nitrous oxide pressure of 10 cm at 
25°C! and assuming diameters of 1.6A and 4.6A 
for oxygen atoms and nitrous oxide molecules 
respectively, 10-* cm is found for L. X is ap- 
proximately 1.25cm. Hence Z=4 X 108. Therefore 


ki;/k3°= (N2O)?/4X 108. (19) 


Reaction (4) must certainly be a very efficient 
process requiring little or no activation energy. 
It will be shown below that (4) is of no im- 
portance even if the steric factor were as large as 
unity ; hence f; exp (— E,/RT) will be assumed to 
be 0.1. Therefore 


X10-%. (20) 


Introducing (19) and (20) into (18) it is found 
that 
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4(O2)(N20)/(NO)?= (exp (2E3/RT))/8 
10%+ (exp (Es/RT))/4X 10°. (21) 


The ratio of oxygen to nitric oxide is 0.28! and 
with a time of 3600 sec. and a final pressure of 
NO of 3.8X10-4 cm, one obtains 


1=(exp (2E;/RT))/8.4X10"" 
+(exp (22) 


The first term on the right (arising from triple 
collisions) obviously can be neglected. Thus 
E;= 9200 cal. 

A comparison of the energy of 14,000 cal. for 
(3) calculated from the data in this article with 
the 9200 cal. obtained from the previous data! 
indicates that practically all of the oxygen 
atoms produced by photochemical decomposition 
of nitrous oxide will form oxygen molecules by 
collisions with the walls. 

Since nitric oxide is known to be formed,'~ 
primary process (2) must be accompanied by 
(5), which will be followed by (6) and (23) 


N+walls=3Noe, (23) 


since the homogeneous termglecular combination 
of nitrogen atoms may be neglected for the same 
reasons as (4). 

Since the number of molecules of Nz and NO 
together formed per quantum absorbed is 1.9;,' 
(6) must occur if (3) is excluded. For the speed of 
(6) to be comparable with that of (23) it can be 
shown by the above methods that E, must be less 
than 10,000 cal. 

From the results of Musgrave and Hinshel- 
wood,‘ the activation energy of (3) must be less 
than 24,000 cal. if this reaction is to compete 
successfully with (15). The 10,000 cal. difference 
between this figure and the minimum found in 
the present article might well explain why 
reaction (3) takes place to a small extent at 
1000°K. 


SEPT 


of tl 

near 

speci 

enery 

pair. 

para 

vith 

that 

to tk 

| the | 

free 

Se 

treat 

latti 

| actic 

sitio 

that 

quar 

gene 

tran 

with 

thec 

pare 

| neat 

q 

and 

3\ 

a 

Soc. 


triple 
Thus 


il. for 
with 
data! 
cygen 
sition 
es by 


d by 


(23) 


ation 
same 


| NO 
1.9;,! 
ed of 
in be 
> less 


shel- 
> less 
pete 
ence 
id in 
why 
t at 


SEPTEMBER, 1938 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 6 


On the Interaction of Atoms in Alloys 
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A generalization of Bethe’s theory of superlattices so as to take into account the interactions 
between atoms which are not nearest neighbors in the lattice is considered. This process will 
not furnish any great improvement in agreement between theory and experiment because it 
predicts, in disagreement with experiment, that many properties of binary alloys are sym- 
metrical between the two constituents. It is further shown to be in disagreement with certain 
fundamental features of the interactions of atoms in alloys. 


ETHE! has given a treatment of the order- 
disorder transformation in binary alloys of 
50 atomic percent composition based upon ideas 
of the interaction between atoms which are 
nearest neighbors in the lattice. If the two 
species of atoms are A and B, then the interaction 
energy between a pair of nearest neighbors is 
Vas, Ven, or Vaz according to the nature of the 
pair. It is assumed that these three energy 
parameters are strictly constants and do not vary 
with the state of order ; and it is further assumed 
that the configurational free energy, i.e., that due 
to the arrangement of the A and B atoms upon 
the lattice sites, can be separated from all other 
free energies and treated separately. 

Several writers?:* have generalized Bethe’s 
treatment to arbitrary compositions while re- 
stricting themselves to the case in which the 
lattice type (but not necessarily the spacing) 
remains the same by assuming that the inter- 
action parameters are independent of compo- 
sition. The results of these investigations show 
that the theories are along the right lines but the 
quantitative agreement with experiment is 
generally poor. In particular, the theory predicts 
that the phase diagrams for the order-disorder 
transformation should be symmetrical about 
fifty atomic percent, a prediction in disagreement 
with experiment. A natural extension of the 
theory consists of introducing three more energy 
parameters describing interactions between next 
nearest neighbors;* indeed even more remote 


wi. A. Bethe, Proc. Roy. Soc. 150A, 552 (1935). 
C. E. Easthope, Proc. Camb. Phil. Soc. 33, 502 (1937) 
and 34, 68 (1938). 
: W. Shockley, J. Chem. Phys. 6, 130 (1938). 
This case has been treated by T. S. Chang, Proc. Roy. 
- London 161A, 546 (1937). 


interactions may be considered. We shall use the 
phrase ‘“‘nth neighbor assumption’ to describe 
the case in which interactions are included be- 
tween all atoms which neighbors of the mth or 
less degree. One intuitively feels that the contri- 
bution of any one atom to the energy of the alloy 
will be dependent upon the arrangement of the 
atoms in its immediate vicinity, i.e., those 
within the neighboring unit cells. For this reason 
we might expect that an mth neighbor assumption 
with a reasonably small value of » would mimic 
quite well the actual situation. It is the purpose of 
this note to point out that such expectations are 
not well founded and to mention some other lines 
of attack. 

In the first place, we shall prove that the 
assumption of interactions between neighbors of 
any order does not interfere with the symmetry 
of the phase diagram about fifty atomic percent. 
Let us denote by Visa, Visas, and Viaz, the 
interactions between nearest neighbors in the 
lattice and by Z; the number of nearest neighbors 
of each atom (these are the quantities used in 
Bethe’s theory) and by Vesa, Voss, and Z2 
corresponding quantities for next nearest neigh- 
bors, etc. In this proof, which is formally identical 
with that given by the writer for nearest neighbor 
interactions,’ it is assumed that the lattice is of a 
simple type, such as simple cubic, face-centered 
cubic, or body-centered cubic so that all lattice 
sites are equivalent as regards their distribution 
of neighbors (for other lattices, such as the 
fluorspar structure, interactions between nearest 
neighbors may give asymmetry). Let us denote 
by Qnaa, Qnes, Qnaz the numbers of pairs of 
atoms which are mth degree neighbors and are 
formed of AA, BB and AB, respectively. Then 
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the contribution of these mth degree interactions 
to the energy E of alloy is: 


En= VnaaQnaat 


If the total number of atoms in the alloy is N 
and if fractions F, and Fz, of them are A and B, 
then 

2A 


If we compare the energy of the alloy with that 
of two crystals of the pure metals, we find that 
the energy due to the mth degree interactions is 


E,'= nQnaB; 
Vin=3(Vaaat Vass) — 


This energy term is obviously symmetrical in A 
and B and so are the terms due to the interactions 
of all other orders. 

Another reason for being skeptical of any of 
the nth neighbor assumptions is that they repre- 
sent emasculated central force treatments of the 
alloy. Let us assume (1) that the atoms interact 
in pairs, the potential energies being Vaa(r), 
Vea(r) and Vaz(r), where r is the distance 
between atoms and (2) that the lattice of sites isa 
rigid unalterable framework. Then the distance 
r, between nearest neighbors is a fixed quantity 
and the same is true of 72, the distance between 
next nearest neighbors, etc. These assumptions 
correspond to the th neighbor assumption with 
Vnaa = Vaa(rn). The fact that this correspondence 
arises directly from (2) constitutes a severe 
criticism of the mth neighbor assumption: for (2) 
is a quite unreasonable assumption denying the 
variations of lattice constant with composition 
which are recognized empirically by Vegard’s law 
and which would otherwise follow from a differ- 
ence in equilibrium radii of the three V’s. 

This reasoning suggests that even the nearest 
neighbor assumption should be modified by 
allowing the interatomic distances to vary, 
entailing corresponding changes in the V’s; that 
is, we should use central force potential energies 
between the nearest neighbors. Hume-Rothery® 
has discussed the formation of a superlattice in 
Fe;Al on this basis: this alloy is characterized by 

having no pairs of Al atoms as nearest neighbors 


where 


5W. Hume-Rothery and H. M. Powell, Zeits. f. 
Krist. 91, 23 (1935). (See also reference 8, p. 40.) 
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and a minimum number as next nearest neigh- 
bors. Although this situation can be readily 
explained by assuming next-nearest-neighbor 
interactions, it may equally well be explained, as 
Hume-Rothery has pointed out, by considering 
the strains produced by the interactions of 
nearest neighbors of different atomic diameters; 
these strains being least when most evenly 
distributed, that is, when the Al atoms are as far 
apart as possible. A. H. Wilson® has carried out 
an analytical investigation of the nearest neighbor 
central force theory for a case corresponding to 
superstructure formation in CuAu. He predicts a 
latent heat and a change of lattice constant with 
order. 

The next-nearest-neighbor assumption without 
allowance for change in atomic spacing has been 
used by Chang,‘ who finds that as compared with 
the nearest-neighbor assumption it gives im- 
proved agreement with experiment by leading to 
higher specific heats near the critical temperature 
of order; and by Bitter,” who shows that it can 
lead to intermetallic compounds. However, the 
predictions of the nearest neighbor assumption 
are greatly altered if change in atomic spacing is 


‘permitted as has been shown both by the 


calculations of Wilson® and by considerations of 
the entropy of lattice vibrations ;* furthermore 
suitable assumptions about the interatomic forces 
lead to the expectation of intermetallic com- 
pounds. From this we conclude that the intro- 
duction of next-nearest-neighbor interactions is 
not required until the nearest-neighbor central- 
force picture has been further investigated. 
Although improvement may be justifiably 
expected to follow the introduction of potential 
energies between nearest-neighbors, it is known 
that the central force picture is inadequate for 
metals. This conclusion is drawn from the failure 
of the Cauchy relations for the elastic constants, 
and the reason for this failure is well understood 
from the work of Fuchs® using the modern theory 
of metals. The work of Mott!® on B-brass has 
shown that for the superlattice formation in this 


6 A. H. Wilson, Proc. Camb. Phil. Soc. 34, 81 (1938). 


7F. Bitter, J. Chem. Phys. 6, 339 (1938). 

8 For a discussion see F. C. Nix and W. Shockley, Rev. 
Mod. Phys. 10, 1 (1938). 

9K. Fuchs, Proc. Roy. Soc. 157, 622 (1936). 

10N. F. Mott, Proc. Phys. Soc. 49, 258 (1937). 
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alloy the central force picture is untenable. He 
finds that the atoms do not interact directly in 
pairs and that the major part of the energy of the 
superlattice arises from electrostatic interactions 
between the Cu and Zn atoms, which are 
unequally charged. This difference in charge 
depends upon the manner of occupation of the 
neighboring sites of each of the two atoms; it is 
greatest when they are occupied as in a super- 


lattice and least when each atom is surrounded 
by neighbors of its own species, the atom tending 
in the latter case to be neutral as in the pure 
metal. Thus the interaction between any two 
atoms depends very considerably upon the other 
atoms which surround them and this is a com- 
plete refutation of the central force picture. 
Mott has suggested an interpretation of the high 
specific heat of 8-brass on this basis. 
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The Vibrational Spectrum and Thermodynamic Properties of Nickel Carbonyl 
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The infra-red absorption spectrum of nickel carbonyl has been measured in the gas phase 
from 4 to 23u, and in the liquid phase from 1 to 8u. Normal coordinate treatments of both the 
square and tetrahedral models of the AB,C, molecule have been applied to the infra-red and 
Raman data; the results give further evidence of the tetrahedral structure of nickel carbonyl. 
An analysis of the vibrational spectrum has been made, based on the tetrahedral model, and 
has been applied to the statistical calculation of the entropy, heat capacity, and free energy 


function. 


INTRODUCTION 


HE theory of directed valence predicts a 

tetrahedral structure for a nickel carbonyl 
molecule consisting of four carbonyl groups 
singly bonded to a neutral nickel atom.' Elec- 
tron diffraction studies? appear to confirm this 
structure. However, the observed nickel-carbon 
distance indicates strong resonance with a 
structure involving nickel-carbon double bonds, 
which may or may not tend to alter the tetra- 
hedral configuration. A square configuration 
has frequently been proposed, in analogy with 
the known square configuration of quadri- 
covalent compounds of bivalent nickel.2 A 
qualitative analysis of the Raman spectrum of 
nickel carbonyl was interpreted as favoring the 


* National Research Fellow in Chemistry. 
t Present address: Brown University, Providence, R. I. 
; L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931). 
. O. Brockway and P. C. Cross, J. Chem. Phys. 3, 
828 (1935). 


square structure ;? however, it has been shown 
that the selection rules were violated in this 
analysis.* 

We have therefore undertaken a quantitative 
analysis of the vibrational spectrum, including 
infra-red data, considering both the square and 
tetrahedral models. It was our hope to obtain a 
new and independent determination of the 
structure; indications favoring the tetrahedral 
configuration were found, but the square con- 


2100 1300 
Frequency 


Fic. 1. Absorption spectrum of nickel carbonyl vapor, 
rocksalt region (20 cm path length, 300 mm pressure). 


. Duncan and J. W. Murray, J. Chem. Phys. 
2, 636 (1934). 
4E. B. Wilson, Jr., J. Chem. Phys. 3, 59 (1935). 
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TABLE I. Observed and calculated infra-red active frequencies. indic. 
unce! 
OBSERVED CALCULATED OBSERVED CALCULATED ance 
INTENSITY ASSIGNMENT (cm™!) INTENSITY (cm™!) ASSIGNMENT 
ne 
90 ws (8) 920 50 
160 (4) 960 (8) pea 
180 2ws (8) 990 4 990 wi 
300 we (v) 1000 wate (6) 
370 wstwe (v) 1060 1 1060 wit2w3+w7 
390 wstwe (vr) 1080 
475 20 470 witws; (vr) 1140 0 1130 Bc 
490 w3tws (5) 1150 w3t+2w; 
510 wstws (5) 1340 1 1340 ws tions 
545 60 540 w; (8) 1590 (liquid) 1590 W3— 4 the 
550 (rv) 1610 0 1620 3w7 
610 (6 iqui uce 
630 (8) || 1682 20 1670 1 
680 10 680 witwe (v) 1740 (liquid) 1750 Ws — ws were 
720 wetws 1960 Wg— Ws 
760 70 760 watwe (rv) 1980 r= 
790 10 800 witwys (vr) 2050 100 2050 ws 
840 50 840 wetw7 (rv) 2120 
840 (8) 2150 20 2140 d 
(5) 2230 4 2230 an 
TABLE II. Assignment of observed Raman active frequencies. vapor from 4 to 14y. The absorption of the vapor T 
in the region from 12 to 23 was also studied, sket 
INTENSITY ASSIGNMENT using the spectrometer at Harvard University.‘ sele 
82 100 p 70 pa The absorption curves are given in Figs. 1 to 3. indi 
90 ws The frequencies of the maxima of the bands are F 
= given in Table I ; the Raman data of Duncan and forc 
601 4 600 2we Murray? are summarized in Table II. calc 
After our measurements had been completed, moc 
872 2 880 ws two a report of a similar study by Bailey and Gordon fun 
appeared.’ The agreement between their ob- 
2043 30 p 2040 we servations and ours is very satisfactory, although 
2132 5 2130 wotw 
3223 3220 there are some differences. Several of the w eak 
bands which we found are not given by Bailey 
Data from reference 3; D=depolarized, P=polarized, p=partly @nd Gordon, and their relative intensities for Her 
a certain bands do not agree with ours. They re- ith 
solved the band at 680 into three com- in 
figuration cannot be conclusively eliminated on : ; 
‘ . : ponents; in our measurements we find some ith 
the basis of this analysis. sad 
Using the tetrahedral model, a satisfactory T 
assignment of fundamental frequencies was 
made, and thermodynamic properties calculated. 
EXPERIMENTAL 
The infra-red absorption in the rocksalt re- 800 700 300 400 
Frequency (cm7') F 
gion was studied using the spectrometer at the : “ 
University of Wisconsin.’ The absorption of the _F16- 2. Absorption spectrum of nickel carbonyl vapor, KB: 
ie ¥ region (20 cm path length, 200 mm pressure). — 
liquid was studied from 1 to that of the —————— 
cacmmecnonesnimes 6 H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. Wi 
5 P. C. Cross and F. Daniels, J. Chem. Phys. 1, 48 (1933). 6, 197 (1938). These measurements were made by B. L. CG. 
These measurements were made by P. C. Cross, who wishes Crawford, Jr., who wishes to thank Professor Wilson and ¥ 
to express his gratitude to Professor Daniels and to the Dr. Gershinowitz for the use of this instrument and {or ae 
Wisconsin Alumni Research Corporation for the grant their continued interest in this work. on 
which made possible this study. 7 Bailey and Gordon, J. Chem. Phys. 6, 225 (1938). 8 
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indication of such structure, but of such an 
uncertain nature that we can attach no signifi- 
cance to it. Certainly the possible side-com- 
ponents are not comparable in intensity to the 
peak at 680 


NORMAL COORDINATE TREATMENT 


Both the square and tetrahedral configura- 
tions of the molecule AB,C, were treated, using 
the methods of group theory* in finding the 
normal coordinates and selection rules. The re- 
duced representations for the internal vibrations 
were found to be, for the square model (D,,), 


r= 2A 
+ 


and for the tetrahedral model (T 4), 
r= 2A i+2£4+47T24+ Ti. 


The normal coordinates for the two models are 
sketched approximately in Figs. 4 and 5; the 
selection rules for the fundamentals are also 
indicated. 

Frequency expressions, based upon valence- 
force type potential energy expressions, were 
calculated for both models. For the tetrahedral 
model, the following five-constant potential 
function was used : 


+H, + (AB;)? 
+hR*> (Aa; i) (Aa 


Here R; is the distance from the nickel to the 
ith carbon atom, 7; the carbon-oxygen distance 
in the ith CO group, a;; the angle between the 
ith and jth nickel-carbon bonds (a;;°=109°28’) 
and 8; the angle between the ith nickel-carbon 


Percentage Transmission 


5000 4500 2000 «1500: 


3000 2 
Frequency (cm-") 


Fic. 3. Absorption of liquid nickel carbonyl (2 mm path 
length). 


*E. Wigner, Géttinger Nachrichten (1930), p. 133; E. B. 
Wilson, Jr., Phys. Rev. 45, 706 (1934); J. E. Vomatiel ans 
G. M. Murphy, Rev. Mod. Phys. 8, 317 (1936). We are 
grateful to Professor Wilson for making available to us 
Some preliminary work on the AB,C, molecule, and for his 
Suggestions and criticisms. 


and carbon-oxygen bonds (6;°=180°). K, and 
K; are the force constants for stretching of the 
nickel-carbon and carbon-oxygen bonds, re- 
spectively; 7, and H, are force constants for 
distortion of the valence angles ; h is the constant 
for interactions between opposite valence angles 
of the nickel atom, the last summation being 
taken over the three pairs of opposite angles a;;. 


G4, Wo, Ay,R W3,E,R 


WER R-IR R-IR 


Ta,R-IR Wg, 7;, In 


Fic. 4. Approximate modes of vibration of the tetrahedral 
AB,C, type molecule, Tz. 


We, Ey AR 


Wig EyIR IR 


Fic. 5. Approximate modes of vibration of the square 
AB,C, type molecule, Dy, 
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This potential function leads to the frequency expressions given below, in which \;=47°cw,. 
Fundamentals marked R are Raman active, those marked JR are infra-red active, and those 
marked Jn, inactive in both Raman and infra-red. 


2A1;w1, we; R(p=3): 


2E; ws, wa; R: 


W5, WE, W7, Ws; R-IR: 


—3MSr 
+(Ki+4W") 
+NQH, 


2(K,—211;") 


—3Mxr 
+2(K;+H,’’) 
+(2K2+H2) 


(K,—2H,") 
+(K.—NH2) 


—(Ke+2NPH2) 


2W(Ki—2H,") 


—(K2—NH2) —(2K2+H2) 


w,; In: 


=H,+(h/2) 
Hy’ =H,—(h/2) 
f=(R+r)/R 


where 


N=1+(4MSr/MoR) 


Frequency expressions were also calculated 
for the square model, using a similar potential 
function (seven constants). Since this model was 
rejected, it does not seem worth while to give 
these expressions in detail. Fig. 6, which shows 
the relation between the frequencies of the 
tetrahedral and square models, plotted qualita- 
tively on an energy scale, will serve to show 
approximately the relative magnitudes of the 
frequencies. It is to be emphasized that this 
figure is a qualitative one ; the actual separations 
of the levels in the square model will depend on 
the relative magnitudes of the bending force 
constants for distortion in the plane of the mole- 
cule and perpendicular to this plane. 


S=Mo/(Mo+4M+4m) 


and Mo, M and m refer to the masses of the A, B, and C atoms, respectively. 


—(Ki+4H,") 
—~QUH, 


—(K,—2H,"”) —(2K2+H2) 


— U(K2—H2) 


— W(K:+4H,") 
+ U(K2+2PH) 


2(K2—PH}) 


U(K2—H2) 


\=H.(M+mf?)/mM 


P=1-—Wg 
Q=2gSUW 
U=(M+m)/m 
W=4m/M> 


VIBRATIONAL ANALYSIS 


Consideration of the square model 


While we have not found conclusive evidence 
from this analysis against the square configura- 
tion, there are certain features of the spectrum 
which are difficult to reconcile with this model. 
The strong Raman line at 82 cm~ appears to be. 
as Duncan and Murray’ describe it, ‘‘very broad, 
unsymmetrical, and of unusually great intensity 
which suggests that it is complex.’”’ This would 
require, on the basis of the square model, two § 
type frequencies lying near 80 cm-. This is 
supported by the appearance of two series of 
bands; one occurring in the infra-red at 920. 
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Fic. 6. Relation between the fundamental frequencies of 
the tetrahedral and square models for nickel carbonyl. 
The frequencies are located on an energy scale; the mag- 
nitudes are noted, in cm™, for the tetrahedral model; the 
locations for the square model are qualitative only. 


990, and 1060 cm~'!, and most satisfactorily 
assigned to combinations involving a g type 
frequency of 70 cm; and the other occurring in 
both the infra-red and the Raman spectra, at 
2040, 2130, and 2220 cm-! (Raman frequencies), 
and similarly suggesting a g type frequency of 
90 

There can be, however, but one g type fre- 
quency of this magnitude. The discussion of the 
Raman data given below in connection with the 
tetrahedral model applies also to the square 
configuration; thus the polarized lines at 382 
and 2043 cm! must be assigned to w; and we, 
and the low lying frequency of about 80 cm~ to 
ws. Several choices may be made as to the other 
assignments of the g type frequencies; we have 
been unable, however, to find a set of assign- 
ments which would give reasonable force con- 
stants and yield another g type frequency near 
80 cm-. (See Fig. 6.) 

These indications, together with the electron- 
diffraction evidence, and the satisfactory assign- 
ment obtained with the tetrahedral model, have 
led us to abandon further consideration of the 
Square configuration. 


Consideration of the tetrahedral model 


The totally symmetric radial frequencies w: 
and we should appear in the Raman effect as 
strong, polarized lines; they must therefore be 
assigned to 382 and 2043 cm™. We may regard 
the broad 82 cm™ line as composed of two com- 
ponents at 70 and 90 cm“, the frequencies indi- 
cated by the series mentioned above. The fre- 
quency expressions then favor the assignment 
w3=70, w,=90 cm. We may then assign w, to 
the next strongest Raman line, at 463 cm~'; or we 
may regard this line as a combination of w;+w; 
and assign w, to 601 cm“. From either assign- 
ment we may calculate K;, Ke, H,’, and He, and 
hence should be able to calculate the other five 
frequencies, adjusting the interaction constant 
h to fit the infra-red spectrum. If we take 
w4= 601, then we must adjust h to make w;= 760, 
since the frequency expressions are such that w7 
must lie above w, (see Fig. 6), and no infra-red 
band between 600 and 760 cm is sufficiently 
intense to encourage the assignment of a funda- 
mental within that region. This would lead to 
a value of about 220 cm™ for ws and about 
580 cm! for ws; these values give no possible 
fundamental or binary combination for the 
intense band at 545 cm™', and hence must be 
rejected. If we take w,=463 cm, on the other 
hand, we are led to the following force constants: 


Ki= 2.52 (105) dynes/cm. 
15.89 
0.089 
0.377 
h=—0.010. 


Table III gives the values of the fundamental 
frequencies thus calculated, and also the assigned 
values, which are used in the calculation of com- 


TABLE III. Fundamental frequencies. 


CALCULATED ASSIGNED 


ws 


Frequencies in parentheses were used to determine force constants. 
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bination frequencies and are to be regarded as 
the ‘‘true’’ fundamental frequencies of the 
molecule. The agreement is satisfactory when the 
approximate nature of the potential function is 
considered. 


TABLE IV. Thermodynamic values. 


Temp. (°K) S* 


97.0 
99.1 
102.7 
107.7 
112.2 
116.4 


—(F°—Eo°)/T Cy 


73.4 32.8 
74.8 33.5 
77.4 34.5 
80.8 35.8 
84.1 36.8 
87.1 37.7 


The calculated values for all infra-red active 
fundamentals, first overtones, and binary com- 
binations which lie below 2250 cm™ are given in 
Table I, compared with the observed spectrum. 
(The table also contains calculated values and 
assignments for certain higher combinations and 
difference tones, which are needed to completely 
fit the observations.) If we divide the nine 
fundamental frequencies approximately into 
those of “‘valence’’ type (w1, we, ws, ws) and those 
of “deformation” type (ws, ws, 5, #7, w9) we 
find an interesting intensity pattern exhibited in 
Table I. The combinations which involve at 
least one valence type component (marked » in 
Table I) are more strongly infra-red active than 
the pure deformation type combinations (marked 
5), which indeed are so weak as not to appear. 
The single real exception to this generalization 
is the fact that no band is found at 720 cm™, the 
calculated value of ws+wy. The possible sig- 
nificance of this intensity pattern is an interest- 
ing subject; however, in view of our present 
inadequate knowledge concerning vibrational 
intensities, we prefer not to attempt any in- 
terpretation. 

It will be noted that the difference tones arising 
from the two low lying fundamentals, with cal- 
culated frequencies at 1960 and 1980 cm~, are 
not found. This is rather surprising in view of the 
appearance of other difference tones. It should 
be pointed out, however, that this is not to be 
regarded as disagreeing with the particular as- 
signment of frequencies given here; any assign- 
ment must include a Raman active frequency 


JR. 


CROSS 


F. 


near 80 cm“, which should give a difference 
frequency at about 1970 cm“. 

The assignments and calculated frequencies, 
which account for the Raman spectrum, are 
given in Table II, together with the observed 
frequencies.® 

The value found for the C—O bond stretching 
force constant, 15.89X10° dyne/cm, lies be- 
tween the value® associated with the C—O 
double bond, 12.3105, and that associated with 
the C—O triple bond, 18.6 X 10°; it is very close 
to the value found in carbon dioxide, 15.24 x 10°. 
This is in agreement with the resonance hybrid 
character of the bond suggested by Brockway 
and Cross? on the basis of the interatomic 
distances. 

Both the C—O and the Ni—C bond force 
constants, when used in Badger’s relation’ be- 
tween force constants and interatomic distance, 
give fairly good agreement: Ni—C distance, 
calculated 1.85A, observed 1.82A; C—O dis- 
tance, calculated 1.17A, observed 1.15A. 


THERMODYNAMICS 


By standard methods of statistical mechanics, 
the virtual entropy, free energy function, and 
heat capacity of nickel carbonyl (perfect gas at 
1 atmos. pressure) have been calculated for the 
temperatures 298.1, 316.1 (boiling point) 350, 
400, 450, and 500°K. The moments of inertia 
were calculated from the interatomic dis- 
tances? Ni—C=1.82A, C—O=1.15A;; the prod- 
uct (ABC)=5.03 (10-") g* cm®. The thermo- 
dynamic values are given in Table IV. 

In view of the fact that several of the funda- 
mentals, including the low lying frequencies, 
were not directly observed, but assigned on the 
basis of combination assignments and frequency 
calculations, these thermodynamic values can- 
not be taken as highly accurate. After consider- 
ing the effects of various frequency changes, 
however, we believe that the error is not greater 
than 2 percent in the entropy and free energy 
function, and somewhat less in the heat capacity: 


%H. W. Thompson and J. W. Linnett, J. Chem. Soc. 
1384 (1937). 
10 R. M. Badger, J. Chem. Phys. 3, 710 (1935). 
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Vibration Spectra and Molecular Structure 
IV. The Infra-Red Absorption Spectra of the Double and Single Molecules of Formic Acid 


Lyman G. BONNER* AND ROBERT HOFSTADTER 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 


(Received May 21, 1938) 


In the present study the absorption spectrum of formic acid vapor between 1 and 15.5y has 
been investigated under a variety of temperature conditions. The spectrum at room tem- 
perature gives essentially the absorption of the double molecules alone, while that at 140°C 
gives the absorption of the single molecules. The results are interpreted in terms of the hydrogen 
bonding existing in the double molecule. In addition an attempt to prepare and investigate 


pure formic acid-d is reported. 


S is well known, the vapor of formic acid 
exists, at the vapor pressure of the liquid, 
largely as double molecules at room temperature, 
but is readily dissociated. A temperature of 
150°C, at this same pressure, suffices to almost 
completely dissociate the material to single mole- 
cules. According to electron diffraction measure- 
ments of Pauling and Brockway! the double 
molecules are associated through the OH to givea 
OHO 
structure HC CH in which the O—O dis- 
OHO 
tance is 2.67A. The purpose of the present work 
has been, in part, to study the structure and 
constants of the formic acid molecule both in the 
single and double states, and in part to investi- 
gate the nature of the hydrogen bond. Formic 
acid is one of the simplest and most symmetrical 
molecules in which this type of bonding exists 
and seemed, therefore, to be a logical starting 
point. 

Depending on the nature and state of aggre- 
gation of the substance under investigation, 
O-H containing molecules seem to fall into 
three classes: those exhibiting a sharp infra-red 
absorption near 3700 cm-!, those with a more 
diffuse band near 3400 cm-!, and those which 
seem to exhibit no absorption characteristic of 
the O—H bond. The first type is exemplified by 
the vapors and dilute solutions of alcohols, acids, 
etc. To the second belong liquid alcohols and 


*Now at Duke University, Durham, N. C. This work 
was performed while the author was the holder of a 
National Research Council Fellowship in Chemistry. 

19. on and Brockway, Proc. Nat. Acad. Sci. 20, 336 


higher acids and in the third class fall chelate 
compounds of the type discussed by Wulf et al.” 
and the lighter fatty acids. These three classes 
and their probable differences have been dis- 
cussed recently by Buswell, Deitz and Rodebush,* 
by Errera and Mollet‘ and by Badger and 
Bauer.’ Undoubtedly the normal undisturbed 
oscillation frequency of the O—H bond is 3700 
cm and the existence of the other two classes 
is a result of varying types or degrees of hydrogen 
bonding. Absence of the characteristic O—H 
absorption (3400 cm in liquid samples) has 
been used as a criterion for the existence of a 
chelate ring or other structure involving the OH, 
but the OH frequency in these extreme cases is a 
matter of some concern since it certainly cannot 
disappear altogether. Buswell, Deitz and Rode- 
bush and Badger and Bauer have suggested that 
the frequency shift is just sufficient to super- 
impose the OH on the CH frequency, near 3000 
cm. It is the belief of the present authors that 
this contention is supported by our observations 
on the spectrum of formic acid vapor. 


MATERIALS AND APPARATUS 


To obtain pure anhydrous formic acid, 
Kahlbaum’s product was treated by the method 
of vacuum distillation and sublimation prescribed 
by Coolidge® until the vapor pressure reached the 


2 Hilbert, Wulf, Hendricks and Liddel, J. Am. Chem. Soc. 
58, 548 (1936). 

3 Buswell, Deitz and Rodebush, J. Chem. Phys. 5, 84 
(1937); J. Chem. Phys. 5, 501 (1937). 

4Errera and Mollet, Nature 138, 882 (1936). 

5 Badger and Bauer, J. Chem. Phys. 5, 839 (1937). 

6 Coolidge, J. Am. Chem. Soc. 50, 2166 (1928). 
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value given by this author. Coolidge mentions 
the possibility of dissociation to CO and H,0 
and since this question was important in the 
present work, a sample of this vapor in a one 
liter flask was maintained at 150°C for 36 hours. 
At the end of this time there was no perceptible 
pressure increase to indicate a dissociation of this 
nature. 


TABLE I. Positions of dimer and monomer absorption bands. 


DIMER MONOMER 
v(cm™!) v(em~) 
3.29 3080 2.80 3570 
3.95 2530 3.40 2940 
4.20 2380 4.30 2325 
1905 
5.75 1740 $3 1740 
7.40 1350 
8.30 1205 
9.15 1093 9.15 1093 
10.9 917 
15.0 667 15.2 658 


For our absorption cell, a 15 cm length of glass 
tubing was sealed at the ends with rocksalt 
plates and 5 cm lengths of tubing, also closed 
with rocksalt, were placed at each end to provide 
a dead air space. The problem of providing glass 
to rocksalt seals, vacuum tight at 150°C was 
considerable since none of the usual waxes or 
resins can withstand this temperature. This 
difficulty was finally surmounted by placing 
small asbestos washers between the glass and the 
salt plate, clamping the whole assembly tightly 
in a brass frame, and painting the joints with 
glyptal. Although the glyptal softened at the 
highest temperatures used in these experiments, 
the asbestos usually kept it from sucking in. 
Two identical cells were made and mounted side 
by side and the furnace wrapped around both. 
One cell was evacuated to serve as a dummy and 
the other provided with a side tube projecting 
beyond the furnace to hold liquid formic acid at 
room temperature, thereby determining the 
pressure of vapor in the cell. In one set of runs 
this side tube was immersed in ice to reduce the 
pressure in the cell. Vapor pressures of formic 
acid are, according to Coolidge, 44 mm at 25°C 
and 11.2 mm at 0°. The temperature of the cells 
was determined by means of a thermocouple 
inside the furnace. 
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The spectrometer used was the rocksalt prism 
instrument described by Barnes, Brattain and 
Seitz.’? Percent transmission measurements were 
made throughout the region 1-15.54 at-0.1p 
intervals. The accuracy of the band positions 
listed in Table I is therefore limited to about this 
value. 


RESULTS 


Figures 1 and 2 show series of runs with the 
vapor at various temperatures between 25° and 
140°C, and with the side tube in air, and in an 
ice bath, respectively. In the first case the side 
tube was warmed by the furnace, increasing the 
vapor pressure to the figures shown on the curves. 
In Table I are listed the wave-lengths and 
frequencies of all bands considered significant in 
the monomer and dimer spectra. As may be seen, 
there is a marked transition from the complex 
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Fic. 1. Absorption bands in formic acid vapor; side tube 
in alr. 
dimer spectrum at room temperature to the 
comparatively simple monomer spectrum at 
higher temperatures. From the point of view of 
O-—H bond considerations the region 2.5-3.5u is 
of particular interest. In the low temperature, 
completely associated spectrum there appears 
one very broad absorption with center at ap- 
proximately 3.3u. As the temperature is increased 
this band decreases in intensity, changes i 


7 Barnes, Brattain and Seitz, Phys. Rev. 48, 582 (1935). 


PERCENT TRANSMISSION 


|| 
shar 
atiot 
weal 
all | 
It is 
| 
equ 
pos 
con 
5.4 
con 
sent 
inte 
exis 


prism 
n and 
were 
O.1y 
sitions 
it this 


h the 
° and 
in an 
e side 
ig the 
urves. 
and 
ant in 
seen, 
mplex 


ABSORPTION SPECTRA OF FORMIC ACID 533 


shape, and ultimately, at the highest tempera- 
ture, corresponding to nearly complete dissoci- 
ation, resolves into two sharp and comparatively 
weak bands at 2.8 and 3.35y, respectively. The 
C—-H band, at 3.35u in the monomer, is almost 
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90 
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Fic. 2. Absorption bands in formic acid vapor; side tube 
in ice bath. 

certainly very little affected by the association. 
It is, therefore, a reasonable assumption that the 
greatly enhanced intensity in this region in the 
dimer spectrum is due to a shifting of the O—H 
absorption down from 2.8u. It is to be noted that 
the monomer O—H, at 3570 cm", falls, as might 
be expected, in the region customarily assigned to 
O-—H in vapors and dilute solutions. 

Using the method of rigid groups and the 
equations presented in a previous paper® it is 
possible to calculate C—O and C=O force 
constants for the monomer. We obtain for these 
5.4 and 11.8105 dynes/cm, respectively. A 
considerable amount of evidence already pre- 
sented*® indicates that the normal value of the 
C-O force constant is not far from 4.5. We 
interpret the above value as indicating the 
existence of single-double bond resonance in the 


* Bonner, J. Chem. Phys. 5, 293 (1937). 


monomer, in accord with the prediction of 
Pauling and Sherman.’ This interpretation is 
further borne out by an apparent shortening of 
the C—O bond in formic acid, recently found by 
Bauer and Badger.'® The amount of resonance 
energy must, however, be less than that pre- 
dicted by Pauling and Sherman. 


DEUTERO-ForRMIC ACID 


Since the intensity arguments presented above 
are perhaps not sufficient proof of the coincidence 
of C—H and O—H —O frequencies in this type of 
molecule, an attempt was made to prepare the 
monodeutero derivative HCOOD. In this mole- 
cule, of course, the C—H frequency should 
remain unchanged and the O—D—O band 
should differ by a factor of approximately root 
two from its position in the normal molecule. The 
reduction to normal intensity of the C—H band 
and the appearance at approximately 4.25y of an 
absorption band would confirm the above- 
mentioned coincidence. 
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Fic. 3. Absorption bands in heavy formic acid vapor. 
The high temperature run shows that this material is 
actually an isotopic mixture of all four possible acids. 


® Pauling and Sherman, J. Chem. Phys. 1, 606 (1933). 
10 Bauer and Badger, J. Chem. Phys. 5, 852 (1937). 
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An attempt was made to prepare HCOOD by 
passing D.S (prepared from D.O and Al,S;*) 
through a column of powdered lead formate at 
100°C. The product was purified by distillation 
and the spectrum obtained over a temperature 
range as before. These spectra are reproduced in 
Fig. 3. It is to be noted that the highest tempera- 
ture curve shows the presence of about equal 
amounts of C—H, C—D, O—H and O—D in the 
monomer. From this we conclude that adsorption 
and exchange has taken place on the surface of 


*We are indebted to Dr. John Turkevitch, of the 
Princeton Chemistry Department for supplying us with 
the pure Al.S3. 


the hot lead formate powder and that the dimer 
curves represent the absorption of a statistical 
mixture of the four possible formic acids. That 
the exchange took place in the preparation and 
not in the cell was demonstrated, we believe, by a 
room temperature run in which the formic acid 
was placed ina cell sealed with glyceryl phthallate 
resin. The curve was identical with that obtained 
using a cell sealed with asbestos and glyptal. 

As a result of further experiments by one of us 
(R. H.) it has been found possible to prepare the 
pure HCOOD by a different method and to 
observe its spectrum. These results will be 
published very shortly. 
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The infra-red absorption spectra of CH;COOH and CHsCOOD between 1 and 15.5u have 
been obtained with a rocksalt prism spectrometer. The association of CH;COOD has been 
studied and a value of the heat of dissociation obtained. Positions of the hydrogen and deu- 
terium bonded O—H—O and O—D-—O frequencies are given. Making certain assumptions, 
a value for the shorter O—H and O—D distances in the dimer rings is found to be 1.07A. 


INTRODUCTION 


T room temperature the vapor of acetic 

acid consists chiefly of double molecules 
whose structure may be assumed to be similar to 
that of formic acid. The structure frequently 
assumed is the following: 


. ..) Preliminary accounts of this work were given at the 
Washington meeting of the American Physical Society, 
April 1938, and in a Letter to the Editor, J. Chem. Phys. 

6, 110 (1938). 


in which the two single molecules are linked by 
two hydrogen bonds. Several years ago Gillette 
and Daniels? studied the effect of association on 
the infra-red absorption spectrum of acetic 
acid vapor. By raising the temperature they 
were able to obtain the spectrum of the single 
molecules since from recent vapor density 
measurements’ it is known that dissociation into 
single molecules is almost complete at about 
150°C. Shifts and changes of intensity of certain 
bands were reported. However, no changes are 
mentioned in the high frequency region where 
vibrations involving hydrogen atoms are known 


2Gillette and Daniels, J. Am. Chem. Soc. 58, 1139 


(1935). 
3F,. H. MacDougall, J. Am. Chem. Soc. 58, 2585 (1936). 
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to occur. Since the dissociation involves the 
breaking of two hydrogen bonds, one would ex- 
pect some changes here if there is an infra-red 
active band associated with the hydrogen bond. 
In order to obtain more information concerning 
this point it was thought desirable to study 
d-substituted acetic acid (CH;COOD). Such a 
substitution would alter the frequency of a 
band associated with the hydrogen bond, shift- 
ing it to about 4.3y if this hydrogen vibration 
originally lay in the region of C—H valence 
vibrations.’ The suggestion that the hydrogen 
bond vibration appears fused with the C—H 
valence vibration was made by Buswell, Deitz 
and Rodebush.} The close proximity of these two 
bands has also been mentioned by Badger and 
Bauer.* Whereas the “‘free’” O—H valence vibra- 
tion (fundamental at 2.754) has been observed, 
only weakly if at all, in the spectrum of the 
double molecule, Badger and Bauer report the 
second harmonic band to be characteristic of the 
single molecule of acetic acid. Errera and 
Mollet? have found corresponding results in the 
liquids of some fatty acids. 


EXPERIMENTAL 


Heavy acetic acid was prepared in vacuum by 
reacting Kahlbaum acetyl chloride with 99.6 
percent pure heavy water. This method was 
previously employed by Engler.’ The DCI as 
well as a small excess of acetyl chloride were 
chiefly removed by pumping. The CH;COOD 
was further purified by three distillations. In 
each case all but a middle fraction was rejected. 
The vapor pressure of the final sample at low 
pressures (4 to 15 mm), agreed with that of 
Kahlbaum glacial CH;COOH. In order to 
determine the ratio of deuterium to hydrogen in 
the acid position, the heavy acetic acid was 
reacted with magnesium and the evolved D.— He 
mixture was analyzed in a mass spectrograph.® 


‘They report that the band at 2985 cm™, which they 
assign to the CH bond, is not affected by increase of 
temperature. 

* Buswell, Deitz and Rodebush, J. Chem. Phys. 5, 84, 
501 (1937). 

* Badger and Bauer, J. Chem. Phys. 5, 839 (1937). 

; Errera and Mollet, J. de phys. et rad. 6, 2811 (1935). 

: W. Engler, Zeits. f. physik. Chemie B32, 471 (1936). 

We wish to thank Drs. Hipple and Delfosse for kindly 
making these analyses. 
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This analysis yielded a D/H ratio of 95+ about 
0.1 percent. 

The cells employed were constructed from 
Pyrex cylinders, 22 cm long, having ground glass 
flanges at either end to which cleaved rocksalt 
plates were sealed by means of a thin coating of 
glyptal No. 1201. This technique differs from 
that of Ebers and Nielsen'® only in the use of 
cleaved instead of polished rocksalt windows. 
The glyptal was usually air dried over night and 
then baked out for several hours. A second coat 
of glyptal was then applied and treated in a 
similar manner. Both the cell for the substance 
and a similarly constructed vacuum cell were 
placed in a furnace provided with rocksalt 
windows. A side tube on the sample cell was led 
outside the furnace through a hole in orfe of the 
end plates. Before distilling the liquid to be 
studied into the side tube, the sample cell was 
baked out at about 140°C and simultaneously 
pumped for several hours. Throughout these 
experiments, the vapor was in equilibrium with 
liquid and its pressure was consequently the 
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Fic. 1. Absorption spectrum of heavy acetic acid. Length 
of cell 22 cm. 


10 Ebers and Nielsen, J. Chem. Phys. 5, 822 (1937). 
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vapor pressure of the liquid. A temperature 
bath surrounding the side tube permitted the 
pressure to be varied. Pressures were measured 
by means of a mercury manometer. A thermo- 
couple was inserted in the furnace to measure 
the temperature which did not vary by more 
than about five degrees over the length of the 
cells. The furnace could be raised and lowered 
on ways so that it was possible to place either 
of the cells in the path of the radiation from the 
source. The technique described prevented the 
condensation of liquid on the windows. The 
rocksalt prism spectrometer employed has been 
previously described." In obtaining the infra-red 
curves the transmission of the sample cell was 
compared to that of the vacuum cell at every 
tenth of a yu, several readings being taken per 
point whenever the transmitted energy was low. 
Certain band centers were located by taking 
readings about every fiftieth of a yw. All the 
spectra shown, except the exchange spectrum, 
have been obtained more than once, in particular 
at different pressures. Room temperature runs 
in cells sealed with Picein instead of glyptal 
gave the same spectra. Likewise, room tempera- 
ture runs of samples that had previously been 
run at high temperatures showed no changes. 


RESULTS AND DIsCUSSION 


1. By the method described we have obtained 
the spectrum of heavy acetic acid at 24°C and 
at 139°C. These results are shown in Fig. 1. 
From dissociation data the calculated degrees 
of dissociation for these two cases are 9 percent 
and 97 percent, so that these curves very nearly 
represent the dimer and monomer spectra. It is to 
be noted that the band appearing at 4.35y 
(2299 cm-') is very sensitive to changes in 
temperature. This fact alone makes it highly 
probable that this band is to be attributed to the 
deuterium bonded O—D-—O vibration, since 
such a high frequency fundamental can be 
assigned to no other bond. If this frequency is 
increased by approximately a factor v2 which 
is demanded by a hydrogen replacement one 
would expect this band to be shifted to about 
3.1 in light acetic acid, and to be temperature 


11 Barnes, Brattain and Seitz, Phys. Rev. 48, 582 (1935). 
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sensitive as well. On comparing this expected 
behavior with the results reported by Gillette 
and Daniels,’ it was felt desirable to restudy 
light acetic acid. The results of this investigation 
are shown in Fig. 2. The degrees of dissociation 
are 9 percent and 98 percent for the low and high 
temperature runs, respectively. The positions of 
the observed bands are in good agreement with 
those found by Gillette and Daniels with the 
exception of the band at 3.204 (3125 cm™) which 
they report at 3.354 (2985 cm). As the curves 
show we again find a marked temperature de- 
pendence in conformity with the expectation 
that this frequency is to be associated with the 
hydrogen bonded O—H —O vibration.” 

As the temperature is raised the ratio of 
single to double molecules increases and while 
the O—H—O or O—D—O bonded frequency 
diminishes in intensity, the ‘free’? O—H or 
O-—D appears. These frequencies appear at 
2.754 (3640 and 3.774 (2653 re- 
spectively. The more or less normal values of 
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Fic. 2. Absorption spectrum of light acetic acid. Length 
of cell 22 cm. 


12 Buswell and Rodebush, J. Chem. Phys. 5, 296 (1938), 
from their studies on the carboxylic acids, believe that the 
hydrogen bond absorption lies near 3.2. 
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these frequencies suggest that the H or D atoms 
are linked to only one oxygen atom in the single 
molecule. 

In both light and heavy acetic acids the 
monomer and dimer spectra are quite different. 
Apart from the changes already mentioned, 
which are concerned with hydrogen vibrations, 
we find numerous changes in intensity, the dis- 
appearance of certain bands, the appearance of 
others, as well as shifts. 

In Table I we have tabulated our frequency 


TABLE I.* 


PROBABLE 
ASSIGN- 
MENT MONOMER 


3840 + 140 


MonoMeR | DIMER 


3640+70 
3125419 
3030 +45 ? 


3125 +50 


2299411 
2127425 
173646 


1387 +4 
131549 
1258 +8} 


1058+5 


3175470 
2653415 
2299+ 15t 
1770+6 
1377 +3 
1266+8 
1176+6§ 
1042 +10 


95345 


177047 
138145 
126648 
117646 
1053 +10 


985 +10 


173949 
1493 +20} 
141844 


128248 
1176+6 
1058 +5 
1010+5 


93847 
820+ 10 


*The italicized C—O, 5|0—H and C—C assignments have been 
suggested by preliminary work on CD;COOD. 
+ Probably due to a small amount of dimer. 
 f Seren due to a small amount of monomer. 
Possibly due to five percent impurity of CHsCOOH. 


measurements of the more prominent bands 
found in CHs;COOH and CH;COOD with rather 
generous estimates of probable errors. The com- 
plexity of all these spectra make impossible a 
complete assignment of bands at the present 
time. However, certain bands may be identified 
with confidence. These identifications have been 
noted in the table. 

In light acetic acid three bands shift apprecia- 
ably. The C=O frequency which appears at 
1739 cm— in the dimer shifts to 1770 cm— in 
the monomer, a change of 31 cm-. A slightly 
larger shift occurs for the band at 1418 cm- 
in the dimer. The fact that so large a shift is 
observed in the latter case is surprising if this is 
really a°C—H bending vibration. The band at 
1282 cm— shifts by 16 cm-, appearing at 1266 
‘tm in the monomer. Gillette and Daniels 
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Fic. 3. Exchange spectrum. Length of cell 15 cm. 


238 4 


report two bands in the dimer spectrum at 1740 
and 1435 cm which shift by 46 and 37 cm“, 
appearing in the monomer at 1786 and 1398 cm“, 
respectively. They also find a band in the dimer 
at 1296 cm~ which is reported at 1288 cm™ in 
the monomer. Thus we find qualitative agree- 
ment in the directions and approximate magni- 
tudes of the shifts although the band centers are 
in slight disagreement. In the region beyond 9.1 
in the monomer spectrum we find the appearance 
of a new band at 985 cm“. 

It is of interest to note the great intensity of 
the O—D—O band in CH;COOD, which sug- 
gests the highly polar character of the hydrogen 
bond." A smaller band occurs at 2127 cm™ on 
the long wave-length side of the O—D—O band. 
A similar result has been found in the spectrum 
of formic acid.“ The O—H—O band at 3125 
cm is likewise of great intensity and shows an 
asymmetry on its long wave-length side. The 
asymmetry can possibly be correlated with the 
2127 cm~ band. The C=O frequency in heavy 
acetic acid is observed at the same position as in 
light acetic acid and shifts in the same direction 
by the same amount, within the experimental 
error. The 1377 cm band shifts by the small 
amount of 10 cm~. The new band in the mono- 
mer spectrum at 953 cm-! probably corresponds 
to the 985 cm frequency of the CH;COOH 
monomer. Among the differences between the 
dimer spectra of the light and heavy acids one 


18M. L. Huggins, J. Org. Chem. 1, 407 (1936). 
14 R. Hofstadter, Paper No. VI of this series. 
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further change appears especially noteworthy. 
The band at 938 cm~! in CH;COOH appears to 
shift to 667 cm in CH;COOD and therefore 
seems to be due to a vibration directly involving 
a hydrogen atom. 

In one run upon CH;COOD some asbestos 
paper was present in the form of washers between 
the windows and cell. After heating the vapor 
in this cell to 124°C the temperature was re- 
duced to room temperature and a run taken. 
The curve obtained is shown in Fig. 3. This is 
clearly the result of an exchange because the 
superposition of the spectra of light and heavy 
acetic acids shows all the bands observed in this 
spectrum. It is not clear whether the exchange 
occurred catalytically between the methyl and 
acid group or is a result of impurity within the 
asbestos. 

2. The heat of dissociation of molecules held 
together by hydrogen bonds can -be obtained 
from the dependence of the equilibrium constant 
on temperature as determined by means of vapor 
density measurements. According to Fenton and 
Garner” the heat of dissociation of acetic acid is 
13,790 cal. per mole, whereas Nernst and von 
Wartenburg"® find 15,000 cal. per mole. Mac- 
Dougall* made what is probably the most careful 
determination and obtained 16,400 cal. per mole. 
Since there are two O—H-—O bridges in the 
dimer ring this last value corresponds to a disso- 
ciation energy of 8200 cal. per mole for each 
hydrogen bond. Recently Badger and Bauer'’ 
investigated the dependence on temperature of 
the intensity of the second harmonic of* the 
O-—H band in acetic acid. They remark that 
because of the difficulty of making photometric 
determinations of intensity, their individual 
determinations do not claim the accuracy of 
MacDougall’s measurements. However, in the 
region where their accuracy was greatest, the 
best straight line through their points does not 
deviate greatly from MacDougall’s which is 
represented by 


logio K »=11.789—3590/T. 
We have made an attempt to find the heat of 


% Fenton and Garner, J. Chem. Soc. 694 (1930). 

16 Nernst and von Wartenburg, Zeits. f. Elektrochem. u. 
angew. physik. chemie 22, 37 (1916). 

1” Badger and Bauer, J. Chem. Phys. 5, 605 (1937). 
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dissociation using one of the fundamental bonds 
of CH;COOD. Since the hydrogen bonded 
O—H-—O frequency falls in the region of C—H 
valence vibrations *: ®' '8 it was not possible to 
investigate light acetic acid by means of its 
fundamental bridge frequency. The substitution 
of a D for an H atom in the acid position shifts 
this frequency, now the deuterium bonded 
O—D-—O frequency, to a region unoccupied by 
other bands, making it possible to study the 
dependence of association on temperature in 
CH;COOD. 

Setting the spectrometer at the center of the 
2299 cm band, intensity measurements were 
made for various total pressures at given 
temperatures between room temperature and 
145°C. In order to obtain a greater percentage 
of dimer in the vapor than we have worked with, 
it would have been necessary to cool the vapor 
below room temperature. This would have been 
very inconvenient with our experimental ar- 
rangement. However, if one were to cool the 
vapor to 5°C the degree of dissociation would 
only change from 9 percent to 6 percent. Since 
it was not possible to obtain 100 percent dimer, 
we could not directly determine the dependence 
of transmission on the number of grams of dimer 
in the cell. This dependence can be obtained by 
means of the ‘free’? O—D band at 3.77u, 
since at or slightly higher than 145°C the vapor 
is almost completely dissociated. However, the 
resolving power of the instrument employed was 
not great enough to study this band; the 
neighboring 4.354 band affected the results. 
Instead, we assumed MacDougall’s values for 
the degree of dissociation a at 38.5°C and 
calculated the partial pressure of dimer Pp from 
the relation 


in which » is the total pressure. Since our total 
pressures were low, rarely greater than 22 mm, 
we may apply the gas law 


PpV=NpRT, 


to find the number of moles of dimer in the 
volume of the cell, V, at a given temperature 7. 


18 Bonner and Hofstadter, J. Chem. Phys. 6, 531 (1938). 
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Fic. 4. Equilibrium constant for CH;COOD. 


This leads to a knowledge of the dependence of 
the transmission on pp=Np/V since we experi- 
mentally determined its dependence on p at a 
given temperature. Employing this calibration 
curve it was possible to find a from the relation 


a=(p—ppRT)/(p+ppRT). 


At each temperature the dissociation constant 
K,, can be evaluated from the equation 


K,=40%p/(1—a2). 


K, is defined as py?/pp in which py and pp are 
the partial pressures of monomer and dimer, 
respectively. Fig. 4 gives the plot of logio K, 
vs. 1/T. The heights of the vertical lines give the 
spread in logig K,, and contain on the average 
about five points. From the slope of the curve 
one finds the heat of dissociation to be 15,900 
cal. per mole or 7,950 cal. per mole per bond. 
This is to be compared with MacDougall’s 
value 8200 calories. Our experimental error 
estimated at roughly +700 cal. is too large to 
distinguish between the heats of dissociation of 
light and heavy acetic acids. 

The use of a deuterium substitution for in- 
vestigating dissociation processes spectroscopi- 
cally may be of some value in examining other 
aliphatic or aromatic compounds. 

3. The electron diffraction experiments of 
Pauling and Brockway'® have shown that in 
formic acid the O—O distance is 2.67A. Since 
this work did not locate the hydrogen atoms in 
the ring some calculations have been made in an 


gah ling and Brockway, Proc. Nat. Acad. Sci. 20, 336 
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attempt to do so. For instance, a quantum-me- 
chanical calculation by Gillette and Sherman”? 
indicated that the dimer ring of formic acid is 
stable if the hydrogen atoms lie on the O—O line 
of centers and at the normal O—H distances. 
The system was shown to be less stable if the 
hydrogen atoms lie midway between the oxygens. 
With the H atoms at the normal O—H distances, 
values of the bridge length and heat of dissocia- 
tion were calculated and found to be in poor 
agreement with experiment. Barnes, Bonner and 
Condon” assumed that the hydrogen lies 
somewhere along the perpendicular bisector of 
the line joining the oxygen atoms considered 
stationary and calculated frequencies for various 
configurations. They obtained an upper limit 
of 2390 cm for the frequency in the case in 
which the hydrogen atom is situated midway 
between the oxygen atoms. This does not agree 
with the experimentally observed bridge fre- 
quency of about 3100 cm, which was not known 
at that time. 

We have made a similar calculation with the 
acetic acid dimer ring considered to have a struc- 
ture identical with that of formic acid. We 
have assumed that the O—-H—O and O—D—O 
groups in the dimer rings of light and heavy 
acetic acids can be treated approximately as 
asymmetric linear triatomic molecules. We 
have further assumed that the force constants 
involved vary inversely with the cube of the 
equilibrium distance according to Badger’s 
formula” 


1.86 X 10° dynes 
(re—dis)* 


’ 
cm 


in which d;; is a constant depending on the 
atoms involved. Calculations were made for 
various O—H and O—D distances. If the shorter 
O-—H or O—D distance is taken as 1.07+0.01A 
the experimental O—H—O and O—D-—O fre- 
quencies, 3125 and 2299 cm~, respectively, are 
checked to within better than 1 percent. This 
distance is to be compared with the normal O—H 
distance of 0.97A. In view of the assumptions 


ass _— and Sherman, J. Am. Chem. Soc. 58, 1135 
a1 Barnes, Bonner and Condon, J. Chem. Phys. 4, 772 


(1936). 
2 R. M. Badger, J. Chem. Phys. 2, 138 (1934). 
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made we do not feel that too much confidence 
can be placed in this value for the distance. 
However, if we use the distances 1.07 and 1.60A, 
the stretching force constants appear to be 4.68 
and 0.91 X 10° dynes/cm. 
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The infra-red absorption spectra between 1 and 15.5u of the monomer and dimer of HCOOD 


have been found with a rocksalt prism spectrometer. The position of the vibration characteristic 
of the ‘‘deuterium bond”’ in HCOOD has been determined and permits an approximate calcula- 
tion of the shorter OD distance in the dimer ring. As in CH;COOD the distance is found 


to be 1.07A. 


INTRODUCTION 


Gomes time ago Dr. L. G. Bonner and the 

author? were interested in investigating light 
formic acid vapor (HCOOH) in the infra-red 
with the hope of finding particulars concerning 
the hydrogen bond and molecular structure of the 
acid. Electron diffraction studies of formic acid® 
indicate that the structure of the double molecule 
is that shown below. 


HC? CH 
\ 


1 These results were reported in part in the Washington 
Meeting of the American Physical Society, April 1938. 

* Phys. Rev. 51, 1017 (1937) and Paper No. IV in this 
series. 

3 Pauling and Brockway, Proc. Nat. Acad. Sci. 20, 336 
(1934). 


Since the vapor of the acid consists of an equi- 
librium mixture of double and single molecules‘ 
it was hoped that by raising the temperature and 
thus promoting the dissociation, the absorption 
band or bands associated with the hydrogen bond 
might be recognized and followed as the transi- 
tion between the double and single molecules 
progressed. In agreement with these ideas, pro- 
found changes were observed in the spectrum on 
passing from room temperature to higher temper- 
atures. However, these results alone were not 
sufficient to fix unambiguously either the position 
or relative intensity of the hydrogen bond band 
in formic acid, for a very great decrease in 
intensity of absorption was observed at about 
3.25u where unfortunately the CH commonly 
appears. One could not say whether this change 
was due to a changing CH intensity or, indeed, to 
the hydrogen bond, although the latter possi- 


4A. S. Coolidge, J. Am. Chem. Soc. 50, 2166 (1928). 
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bility seemed more plausible.®: * It was therefore 
decided to investigate the deuterium substituted 
acid (HCOOD). In this case the ‘deuterium 
bond’’ would have a frequency smaller than that 
of the hydrogen bond by a factor of approxi- 
mately v2. If the band at 3.25 were really to be 
connected with the hydrogen bond, then a new 
one ought now to appear at roughly 4.3y. 

To obtain HCOOD the reaction D.S+PbFo: 
=PbS+2DFo was used. Unfortunately only a 
statistical mixture of the light and heavy acids 
was produced. Apparently internal exchange 
occurred on the surface of the PbFo2 column, so 
that the four forms HCOOH, HCOOD, DCOOH, 
DCOOD were finally present in equal amounts. 
CD and OD, of course, behaved exactly as CH 
and OH ; consequently, no more could be learned 
than from a study of HCOOH. Moreover there 
was a great complexity and richness of bands in 
the spectrum. 

Recently it has been possible to prepare 
moderately pure HCOOD, and it is with this 
compound that the present paper is concerned. 


EXPERIMENTAL 

Preparation 

The first attempt to eliminate the mixed forms 
was to reduce the temperature of the PbFo: 
column from 100°C previously used to room 
temperature. This gave no appreciable improve- 
ment so that an entirely different method of 
preparation was used. At the suggestion of 
Professor Taylor and Dr. Turkevich, a simple 
exchange was tried between HCOOH and 99.5 
percent D,O, with subsequent drying of the acid- 
water mixture with anhydrous CuSO,.’? The 
HCOOD samples obtained by this process were 
still not entirely pure, probably because of 
exchange with light water introduced in the 
drying process. However, the mixtures obtained 
by this method contained only the two mono- 
meric forms HCOOD and HCOOH as shown by 
the absence of any CD absorption in the monomer 
spectrum. The ratios of the amounts of heavy to 
light forms, denoted by D/H, of several samples 


Buswell, Deit 
501 (1937), itz and Rodebush, J. Chem. Phys. 5, 84, 


* Badger and Bauer, J. Chem. Phys. 5, 839 (1937). 


"That this exchange goes was proved by Miinzberg and 
Oberst, Zeits. f. physik. Chemie B31, 18 (1935). a 
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were respectively 3.5, 8, 12.4. The first and last of 
these ratios were determined by analysis in a 
mass spectrograph of the He—Dz mixtures ob- 
tained by reacting Mg strips with the acid 
samples. The middle ratio was found by 
comparing the depths of the OD and OH bands 
found in a high pressure run of the monomer 
spectrum. Good agreement between these two 
methods was obtained for the ratios D/H =1.0, 
3.5. This investigation confirmed the earlier 
conclusions regarding the amounts of the four 
forms present in the mixed acid.? 

The HCOOH was supplied by Kahlbaum and 
was distilled until the vapor pressure agreed with 
the determinations of Coolidge. 


Analysis 

All spectra were obtained with the spectrometer 
described by Barnes, Brattain and Seitz. The 
cells, furnace, etc. were all similar to those used 
in an investigation of acetic acid.'® After the high 
temperature runs were taken, the region 1—-5y 


VV 


PERCENT TRANSMISSION 


WAVE-LENGTH IN & 


Fic. 1. Absorption spectrum of HCOOD. D/H=8.0. 
Length of cell =22 cm. 


8] wish to thank Drs. Hipple and Delfosse for kindly 
doing these analyses. 

® Barnes, Brattain and Seitz, Phys. Rev. 48, 582 (1935). 

10 Herman and Hofstadter, Paper No. V in this series. 
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was re-run at room temperature to see whether 
any exchange had taken place in the cell. No 
significant changes were observed other than 
those which were due to differences in running 
conditions. Even after several days there were no 
important changes. This indicates that spon- 
taneous exchange within or between HCOOD 
molecules is slow, if it exists at all. This is in 
agreement with the work of Miinzberg and 
Oberst mentioned above. 


RESULTS AND DISCUSSION 


Figures 1 and 2 show the spectra obtained for 
two samples with the D/H ratios 8 and 12.4. 
Even though the samples were run at different 
pressures one may see from the curves that 
further purification would probably lead to no 
serious changes." Indeed, the sample for which 
the D/H ratio was 3.5 presented a very similar 
set of spectra. The degree of dissociation, a, was 
determined from the dissociation data of Coolidge. 
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Fic. 2. Absorption spectrum of HCOOD. D/H=12.4. 


Length of cell =22 cm. 


" The ratios of dimers containing one H and one D to 
those containing two D’s are only half the above values, 


namely 4 and 6.2. 
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Since the degree of dissociation at room tempera- 
ture and the pressures shown is in the neighbor- 
hood of 18 percent we may consider the lower 
curves to represent the dimer. The upper curves 
are certainly spectra of the monomer for a is very 
close to 1. Table I lists the frequencies found by 
averaging the band centers for the corresponding 
curves in the two figures. 


TasB_e I. HCOOD. 


MONOMER DIMER 
PROBABLE PROBABLE 
em"! ASSIGNMENT ASSIGNMENT 
3003 +60 C-H 3030 +46 C-H 
2667 35 O-D 
2347 35* 2347 33 O-D-O 
2105 30* 2105 22 
1916 16 
1764 16 C=0 1736 15 C=0 
1381 11 
1253 10 
1163 6 C-O? 1181 7f 
1081 157 
1030 5 1030 5 
980 5 985 10t 
953 8 950 8t 
901 77 


* Very small and probably due to small amount of dimer. 
+ Probably due to 7.5 percent impurity of HCOOH. 
t¢ Probably due to small amount of monomer. 


The dimer spectrum of HCOOD shows a 
strong, broad band at about 4.26y (2347 cm“) 
which is overlapped somewhat by a weaker band 
on the long wave-length side at 4.75 (2105 cm™). 
The band at 4.264 must then be that which 
characterizes the deuterium bond in_ heavy 
formic acid. Further proof that this assignment is 
correct is obtained from the dissociation results. 
for this band is absent in the monomer and 
belongs solely to the dimer, where the deuterium 
bonds occur. No other reasonable choice for such 
a high frequency can be made. One sees how small 
the CH bands of the monomer and dimer are. 
The strong absorption at 3.25u found in light 
formic acid may now be explained as being 
chiefly due to the hydrogen bond, in agreement 
with the earlier ideas.2 Correlated with the 
disappearance of the bonded band, we see the 
emergence of a small band at 3.75u (2667 cm”): 
This is to be attributed to a ‘‘free’” OD valence 
vibration. Results similar to those just mentioned 
have been found in light and heavy acetic acids" 
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TABLE II. 


HCOOD 
Calc. Obs. Calc. 


HCOOH 


1744 cm™ | 1764 1741 
1162 1083 
667 572 


and in the liquid of ethyl alcohol where a similar 
association has been assumed to take place.!” 

From the position of the O—D—O band in the 
dimer, the approximate distance between the D 
and the nearest oxygen in the dimer ring may be 
calculated. It is assumed here that the O—D—O 
vibration frequency in the molecule may be 
computed from a linear asymmetric triatomic 
model where the O—O distance is 2.67A, with 
force constants evaluated by application of 
Badger’s rule.“ A similar calculation can be made 
using the O—H —O frequency at 3077 cm. The 
best fit for the shorter O—H, O—D distance is 
1.07A. Using this distance the two calculated 
frequencies check the observed ones to better 
than 2.6 percent. This same calculation for the 
acetic acid dimers furnishes the identical value, 
1.07A.'° 

The monomer spectrum of HCOOD presents 
an interesting detail in the appearance of an 
intense band at 8.6u (1163 cm~'). In comparison 
with one of the fundamentals of HCOOH (1093 
cm), this band is shifted towards higher fre- 
quency. The similar shape and approximate 
position of this band suggest that it is to be 
attributed to the same bond (C—O), since the 
band at 1093 cm in HCOOH has definitely been 
assigned to this vibration. If this assignment for 
the 8.6u band be correct, the shift towards 


® Errera and Mollet, Comptes rendus 204, 259 (1937). 
*R. M. Badger, J. Chem. Phys. 2, 128 (1934). 


greater frequency cannot be explained as a mere 
mass effect, for one would expect a decrease in 
frequency on substitution of a heavier mass in 
the hydroxyl group (see Table II). 

' One may calculate the frequencies of a 
simplified triatomic model of the HCOOD mole- 
cule by the method of rigid groups. This method 
was previously employed in finding the force 
constants for a similar triatomic model of 
HCOOH.’ The force constants yielded in this 
way, 11.8, 5.4, 0.410° dynes/cm, were used to 
compute the new frequencies shown in Table II, 
where for convenience the data on HCOOH are 
also collected. 

In connection with this anomalous behavior it 
is to be noticed that the calculated frequency for 
the C—O bond lies about midway in the fre- 
quency interval between the bands at 8.6u (1163 
cm~') and 10.2u (980 cm~). No adequate expla- 
nation for this behavior has been worked out, 
although there is a possibility that these two 
bands furnish another example of an accidental 
degeneracy. 

No interpretation of the dimer spectrum has 
been made beyond 6.0x. 
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A. The Raman spectrum of dioxane is shown to include eleven lines, some of which have 
been reported by one or both of two previous investigators. B. Dioxane, when irradiated by 
the ultraviolet light from a mercury arc, gives scattered light whose spectrum shows a broad 
continuous band in the visible region. Intensity of the scattered band diminishes rapidly with 
increasing temperature of the dioxane. Shortened time exposures bring out two maxima at 
about 5800A and 4500A. A similar phenomenon has been observed in work on cyclohexane. 


A. RAMAN SPECTRUM OF DIOXANE 


URING an investigation on the conductance 

of solutions of certain organic salts, it 
became desirable to know something as to the 
Raman frequencies of the solvent dioxane, 1-4 
diethylene oxide. Dioxane has been studied previ- 
ously by two independent experimenters, Villars! 
and Wolkenstein and Syrkin,? but their results 
were not completely in agreement, even though 
both used the same region of the mercury 
spectrum for their work. The former reported ten 
frequencies, the latter eleven, with agreement on 
only six values. Since the frequency shifts meas- 
ured in the first work were duplicated almost 
identically in the second, the disagreement 
obviously arose from doubtful assignment of 
shifted lines to proper ‘‘parent’’ mercury 
frequencies. 

The values reported here rest upon measure- 
ments, made in the near ultraviolet, of frequency 
shifts from the pair of mercury lines 3125.6A and 
3131.6A. These two lines, being almost equal in 
intensity, gave modified lines that were paired in 
characteristic fashion, and it was possible to 
relate each shifted line with its unmodified 
“parent’’ without ambiguity. 

Dioxane was prepared by Dasher.* The com- 
mercial product was frozen (freezing point 11.7— 
11.8°C) to remove the major part of the water, 
melted and refluxed with sodium hydroxide to 
remove methylene acetal, and then repeatedly 
refluxed, first with quick-lime and then with 
sodium, until the resultant solvent did not 


1 Villars, J. Am. Chem. Soc. 52, 4612 (1930). 
? Wolkenstein, Syrkin, J. Chem. Phys. 3, 594 (1935). 
3 Paul J. Dasher, thesis, Indiana University (1937). 


tarnish bright sodium. The final product had a 
specific conductivity of 1-2 X10- ohm-, anda 
density of 1.02766 (25°C). 

The spectrograph employed was a Hilger E-3 
quartz spectrograph, having a dispersion in the 
region of the spectrum used (3000-4000A) of 
about twenty angstroms per millimeter. Ex- 
posure time varied from one to eight hours. The 
light source was a quartz mercury arc with a 
luminous column ten centimeters long. Both 
glass and quartz containing tubes were used, and 
their shape and arrangement were according to 
that described by Wood.‘ The temperature of the 
dioxane was controlled by a smooth sheet of 
water that flowed over the tube and was carried 
away by an aluminum reflector below. 

Results are indicated in Table I. Eleven 
frequency shifts were measured, the modified 
TABLE I. Raman frequency shifts for dioxane in wave 


numbers per centimeter. Numbers in parentheses 
indicate relative intensities. 


I II III IV 


Tus DETER- 
MINATION 


288 (0) 290 
493 (1) 490 
837 (5) 837 


1018 (4) 1018 
1113 (3) 1115 
1211 (2) 1213 
1303 (4) 1304 
1450 (5) 1450 
2724 (2) 2722 
2858 (10) 2861 
2976 (10) 2972 


WOLKEN- 


STEIN, SYRKIN AVERAGE 


VILLARS 


291 (0) 
434 (00)? 


519 (00)? 
837 (4) 


427 (3) 
487 (3) 


836 (10) 


1117 (1) 
1214 (1) 


1442 (2) 
2720 (00) 
2864 (3) 
2967 (3) 


1450 (6) 


2862 (10) 
2971 (10) 


«R. W. Wood, Phil. Mag. 6, 729 (1928); Phys. Rev. 36, 
1421 (1930). 
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lines in every Case appearing in pairs having the 
configuration of the originating mercury lines, 
3125.6 and 3131.6A. Measurements were accu- 
rate to +6 wave numbers per centimeter. The six 
values previously reported both by Villars! and 
by Wolkenstein and Syrkin®? were confirmed 
closely; of the remaining five, two agreed with 
Villars’ original report, and three agreed with the 
values given by Wolkenstein and Syrkin. 

In Table I, the first two columns list final 
values taken from the two earlier reports. The 
results of the present determination appear in 
the third column, and shifts verified definitely by 
two independent determinations are given in 
column four. The six frequency shifts confirmed 
in all three cases are italicized. 


B. FLUORESCENCE 


In the course of the determination of the 
Raman spectrum of dioxane, it was observed 
that the spectrum of the scattered light exhibited 
a strong continuous band in the visible region 
when a quartz mercury arc and a quartz con- 
taining tube were used. Reduced time exposures 
brought out distinct bands, with maxima at 
about 5800A and 4500A. The band of shorter 
wave-length was the broader of the two, ex- 
tending well below 4000A. 
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When a glass containing tube of the same 
dimensions was substituted for the quartz tube, 
the continuous background disappeared. A thin 
glass filter produced the same result, an indi- 
cation that the fluorescence was due to the 
irradiation of the dioxane by ultraviolet light 
below 2900A. 

The intensity of the fluorescence was found to 
diminish rapidly with increasing temperature of 
the dioxane. At 12°C (freezing point 11.8°C) 
there resulted a strong continuous spectrum 
extending from 3400A entirely through the 
visible range. With higher temperatures of the 
dioxane, and with all other conditions unchanged, 
the two maxima were evident, and the extent of 
the spectrum was reduced at both ends. At 
temperatures as high as 50°C the long wave- 
length component disappeared, and the re- 
maining band was limited to a region between 
about 4000A and 5000A. These results are shown 
in Fig. 1 which represents three exposures made 
consecutively on the same plate. The only 
factor changed during the series was the temper- 
ature of the dioxane, as indicated. Exposure was 
for three hours. 

Continued irradiation of the dioxane with 
ultraviolet light did not affect the transparency of 
the liquid. Dioxane exposed for a period of three 
days showed the same fluorescence as that 
exhibited by a fresh sample. 

A phenomenon similar to that just described 
was found by Karl Haberl® during an investi- 
gation on the intensity of the Raman lines of 
cyclohexane, CsHi2. He reported that the in- 
tensity of the fluorescence from cyclohexane was 
negligible for temperatures above 20°C, in- 
creasing rapidly with diminishing temperatures 
to 8°C, the liquid freezing at 5°C. He also found 
two maxima, at 4000A and 2900A. 


5 Karl Haberl, Ann. d. Physik 21, 301 (1934). 


| 
hada 
anda 
r E-3 
n the 
A) of | 
Ex 
. The 
itha — 
Both 
1, and 
ng to 
of the 
et of 
urried 
leven 
dified 
ware 
} 
| 7 
Vv. 36, 


SEPTEMBER, 1938 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 6 
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The near infra-red spectra of three liquid crystals have been studied in the isotropic and 
anisotropic phases, and several band shifts have been observed. A detailed study ot the trans- 
mission as a function of temperature was made at various wave-lengths. Pronounced changes 
in transmission were observed at melting and clearing points, the magnitude of the changes 
being greater at short wave-lengths. Rapid variations in transmission occurred in the meso- 
phase. A qualitative explanation of the observed changes is given in terms of Ornstein’s theory 
of scattering. The application of infra-red methods to the study of changes of state is discussed. 


NVESTIGATIONS of the infra-red spectra 

of associated liquids have yielded information 
concerning the nature of the associational groups 
and the ‘‘structure’’ of the liquids themselves. 
The vibrational spectra of the liquids differ from 
those of the corresponding vapors in two prin- 
cipal respects: (1) the frequencies of certain 
atomic groups may be modified, and (2) new 
intermolecular bands may appear. The variation 
in the 3u alcohol bands is illustrative of the first 
effect, while the 4.74, 20u, and 60u bands in the 
spectrum of liquid water are examples of the 
second. It has been suggested! that valuable 
information might be obtained from a study of 
the infra-red spectra of liquid crystals, in which 
the degree of association is known to be exceed- 
ingly high. There are numerous difficulties in- 
volved in making such a study both on account of 
the complexity of the molecules themselves and 
because of the relatively large magnitude of 
the local lattices or swarms. Rawlins and Taylor? 
have reported their inability to detect differences 
between the spectra of the isotropic and aniso- 


TABLE I. Melting and clearing points of compounds. 


MELTING | CLEARING 
Point Point 
Compound I. Iso-amyl p-(p-ethoxy benzalamino) cin- 
namate : 81°C 137°C 
Compound II. Anisal-amino a-methyl-n-propyl cin- 
namic ester 50°C 85°C 
Compound IIL. Iso-amyl p-(p-ethoxy benzalamino) a- 
methyl cinnamate 83°C 90°C 


(Hereafter the compounds will be referred to by number.) 


1F. 1. G. Rawlins, Trans. Faraday Soc. 29, 993 (1933). 

F. I. G. Rawlins and A. M. Taylor, Infrared Analysis of 
Molecular Spectra (Cambridge University Press, 1929), 
p. 63. 


tropic phases of compounds having liquid 
crystalline properties. 

The compounds chosen for study because of 
their low melting points are listed in Table | 
and were prepared by Mr. J. P. Bain of the local 
chemistry department. The melting and clearing 
points listed are taken from the International 
Critical Tables. For compounds II and III the 
melting points could be checked visually, but 
measurements of the melting point of compound 
I were difficult on account of marked super- 
cooling. Cooling-curve data on compound I gave 
the melting point 77°C, indicating the presence 
of an impurity. The orienting action of container 
walls introduces difficulties in the determination 
of melting and transition temperatures. 

The spectrometer previously described’ has 
been modified by the addition of an optical 
amplifier.* The temperature of the sample being 
studied was controlled by means of a small oil 
bath consisting of concentric brass cylinders. 
The rocksalt cell containing the sample was 
placed in a massive brass support in the center 
of the inner cylinder, cell temperatures being 
determined by a thermocouple brazed into the 
support close to the cell. It was impracticable te 
place a thermocouple in the thin absorbing layer 
itself. The control bath, mounted on horizontal 
ways, could be moved into and out of a vertical 
light beam. 

The spectra of the compounds in the isotropic 
and anisotropic liquid states are shown in Fig. |. 


— Williams and R. Taschek, J. App. Phys. 8, 497 
(1937). 

* The writers wish to express their appreciation to the 
A. A. A. S. for a grant which made possible the construc- 
tion of the amplifier. 
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the temperature of the sample was first raised 
above the clearing point and then lowered by the 
desired amount. Cell thicknesses were of the 
order of several hundredths of a millimeter. Since 


% Transmission 


%o Transmission 


% Transmission 


the temperatures being indicated. In every case 


no marked differences of 
spectra of isotropic and 


from 2y to 7.54 and 156°C 


INFRA-RED STUDY OF LIQUID CRYSTALS 


band positions in the 
anisotropic phases of 
compound I were observed, only a part of the 
spectrum of each phase is shown, 108°C data 
data from 7.5u to 12. 


A / 


:0.00H, 


Wavelength 


Fic. 1. The percentage transmission of compounds I, II, and III in isotropic and aniso- 
tropic liquid phases. The temperatures are indicated. 
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In order to observe bands which could not be 
resolved in thick samples, the 113°C data for an 
extremely thin cell were taken between 6.54 and 
11yu. For compound II no transmission curve for 
the anisotropic phase between 6.54 and 7.7 is 
shown, since the spectra of the two phases appear 
to be similar in this region. Although the spectra 
of the compounds are quite complex, the ab- 
sorption arising from the inner vibrational 
frequencies of certain groups can be recognized. 
The bands near 3.34 can be attributed to C—H 
vibrations, the broadness of the bands possibly 
arising from absorption by both aromatic and 
aliphatic groups. The bands near 6.9u and 7.3u 
in the spectra of most hydrocarbons were not 
observed, their presence possibly being masked 
by intense bands in adjacent spectral regions. 
The intense absorption with maximum at 6.25y 
probably has its origin in the vibrations of the 
benzene ring, although the C=O, C=C, and 
C=N vibrations would be expected to give rise 
to bands in neighboring regions. It is difficult to 
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assign the bands at 8.6u, 9.0u, and to 
particular inner vibrations. 

Although many slight shifts in band positions 
were observed, the shift of the 9.05 band in the 
spectrum of compound III is especially notice- 
able. This band was studied at various tempera- 
tures and the transmission curves are given in 
Fig. 2. In the same figure are shown similar data 
for the 3.34 band of compound I; in order to 
determine whether the orienting property of the 
surfaces had any effect on the observed spectrum, 
the curves on the left were obtained with a thin 
film between two rocksalt plates, while for 
those on the right the compound was spread 
upon a single plate. A study of the transmission 
at numerous temperatures shows the position of 
maximum absorption to be 3.27 in the solid, 
3.374% in anisotropic, and 3.30u in the liquid 
phase. The 6.254 band showed a similar shift of 
0.14 to longer wave-lengths in the mesophase. 
Fig. 2 clearly indicates an increase in transmis- 
sion and a change in band shape at higher tem- 
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Fic. 2. The change in the position of certain bands with temperature variation: 
(a) 3.34 band observed with compound I between rocksalt plates; (b) the same for a 
thin layer spread upon a single plate; (c) the 9.05u band for compound III. 
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Fic. 3. Isochromatic curves for compound I showing percentage transmission as a 
function of temperature. The solid curve shows the data obtained while the tempera- 
ture was being raised; the dotted line represents the data taken during cooling. Note: 
Isochromatics for a normal compound such as stearic acid show marked transmission 


increases only at the melting point. 


peratures. The shift of the maximum apparently 
does not depend upon the method of supporting 
the sample. Although minor changes in the 
shape of the transmission curves are observed, 
no conclusions can be drawn as to the effects of 
the orienting action of the plates. The observed 
wave-length variation of the absorption maxima 
cannot be attributed to changes in the tempera- 
ture of the prism, as the calibration was checked 
at each temperature by means of atmospheric 
bands. 

The most pronounced difference in the spectra 
of the isotropic and anisotropic phases was a 
decrease in background transmission in the latter 
phase between 1p and 6y. Beyond 6y there is 
little difference in intensity. The curves for 
compound III in Fig. 1 are illustrative of this 
decrease in transmission, as they were taken 
with the same sample. In order to study the 
variation in background transmission in more 
detail the data represented in Figs. 3, 4, and 5 
were obtained. These curves indicate the trans- 
mission for given wave-lengths as the temperature 


was varied continuously. Readings were taken at 
two degree intervals from room temperature to 
temperatures well above the clearing point during 
heating and cooling ; the rate of heating was ap- 
proximately one degree per minute. All curves 
show an increase in transmission varying from 
5 percent to 60 percent, the most rapid increase 
taking place over an interval including the clear- 
ing point. The transmission in the solid and 
isotropic liquid phases remains practically con- 
stant, as one might expect. Over the temperature 
intervals in which the compounds are in the 
anisotropic phase one or more “‘arrests’’ occur 
in the heating curves for compounds I and II, 
while sharp maxima and minima appear in those _ 
for compound III. In the cooling curves a 
marked “‘hysteresis”’ in transmission is apparent, 
the percentage transmission remaining constant 
until the temperatures of solidification are ap- 
proached ; at these points transmission decreases 
rapidly and subsequently approaches the values 
for the solid or anisotropic liquid. For example, 
the isochromatics for compound I show that on 
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Fic. 4. Isochromatic curves for compound II showing percentage transmission as a 
function of temperature. Solid line represents heating measurements; the dotted line, 


cooling data. 


cooling the sample invariably solidifies. For 
compound II two sets of 14 data are shown, the 
first curve representing data obtained on initial 
heating; the other curves for this compound 
indicate that the low temperature data corre- 
spond to the transmission value for the aniso- 
tropic state. Other measurements show that 
thin films of compound II remain in a super- 
cooled anisotropic state for several days. A 
similar effect was observed in the case of com- 
pound III, although supercooling was not so 
prolonged. The second set of 1p data for com- 
pound III was taken immediately after the 
first set; a rather interesting change may be ob- 
served in that a rise in temperature apparently 
caused the supercooled anisotropic liquid to 
solidify. In all compounds the fotal change in 
-transmission is greatest at 1u and decreases 
rapidly beyond 6u. Most of the data presented in 
the figures are for regions in which strong char- 
acteristic absorption is not present, but iso- 
chromatic curves taken in the center of absorp- 
tion bands are not noticeably different. (See the 
3.374 curve of Fig. 3.) 

Table II lists the temperature intervals in 
which marked transmission changes and arrests 


occur. Although the previous history of the 
sample and rates of heating. and cooling un- 
doubtedly have an influence on the temperatures 
at which critical changes occur, the variations 
caused by these factors are surprisingly small in 
view of the drastic temperature changes to 
which the samples were subjected. Data ob- 
tained with different samples checked excellently. 
The cooling rate had little effect on the ‘‘hystere- 
sis’’ mentioned above. In determining the exact 
temperatures at which “‘arrests’’ and minima 
occur a slow heating rate and a continuous re- 
corder would be desirable. An interesting effect 
observed with compound III is the growth of the 
high temperature minimum with increasing 
wave-length. It is especially noticeable in the 
isochromatics for compounds I and II that 
the transmissions of the solid and anisotropic 
liquid phases approach the same value at long 
wave-lengths. 

The increase of transmission with temperature 
can be explained qualitatively in terms of 
scattering by the swarms present in the mes0- 
phase, the low transmission of the solid probably 
arising from scattering and reflection by small 
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crystals. Einstein‘ and Smoluchowski> have 
derived an expression for the scattering of light 
by dense media. The intensity of the scattered 
raditation is given in terms of the index of 
refraction of the medium, and hence the ap- 
plicability of this expression to the case of 
doubly refracting liquids is questionable. As- 
suming that the Einstein expression can be ap- 
plied and that an average of mo and 1, is a suit- 
able index, one finds that a decreasing scattering 
volume is the chief factor which could account 
for the observed increases in transmission; there | 
are other factors in the Einstein-Smoluchowski 
equation which indicate a decrease in transmis- 
sion with increasing temperature. Approximate 
calculations of the change in transmission with 
wave-length reveal that the observed changes are 
not in agreement with the values given by the 
Einstein relation. 

A theory more readily applicable to the present 
problem has been developed by Ornstein.® Be- 
cause of birefringence and random distribution 
of swarms, gradients of refractive index occur, 

‘A. Einstein, Ann. d. Physik 33, 1275 (1910). 


5M. V. Smoluchowski, Ann. d. Physik 25, 205 (1908). 
°L. S. Ornstein, Ann. d. Physik 74, 445 (1924). 
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which, according to Ornstein, give rise to a 
scattering coefficient proportional to the square 
of the difference between ordinary and extra- 
ordinary indices. Stumpf? and Van Wyk®* have 
measured these indices in the visible region for a 
number of anisotropic liquids and find that the 
difference between the indices is a decreasing 
function of temperature. A discontinuity of index 
occurs at the clearing point. If one assumes a 
random distribution of swarms in the samples 
used in this study, Ornstein’s scattering theory 
can be applied. If Stumpf’s data can be extra- 
polated to the infra-red, the coefficient decreases 
with increasing temperature and becomes zero 
at the clearing point. 

Although probably due to a combination of 
swarm volume change and birefringence varia- 
tion, the greater part of the increase in trans- 
mission with temperature must be attributed to 
the decrease in birefringence (%»—n,)?. Otherwise 
at wave-lengths as long as 10u most of the in- 
crease should take place as soon as the compound 
melts.® Actually, however, at long wave-lengths 


7F. Stumpf, Ann. d. Physik 37, 351 (1912). 

8 A. Van Wyk, Ann. d. Physik 3, 879 (1929). 

®Measurements of the dielectric and diamagnetic 
anisotropy of liquid crystals show that the volume of the 
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Fic. 5. Isochromatic curves for compound III showing percentage transmission as a func- 


tion of temperature. Solid line indicates heating data; dotted line, cooling data. 
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TABLE I]. Temperatures of characteristic 
transmission changes. 


CooLine 
DeEcREASE 


Seconp 
Rise 


INITIAL 


WAVE-LENGTH ARRESTS 


Compound I 
131-141°C 


133-140 
130-139 
133-142 
134-140 
133-141 
134-142 
135-143 


78-88°C 
100-110 


102-88 

111-95 
(slow decrease) 

89-75 
(slow decrease) 


(from solid) 
lu (from liquid 
crystal) 
2.55 


4.27 
4.76 
6.91 
7.26 


the transmissions of solid and mesophases are 
approximately the same and the increase occurs 
at the clearing point.!° The comparatively small 
total change in transmission can possibly be ex- 
plained in terms of a decreasing (and in the 
case of compound I an almost constant) bire- 
fringence factor (mo—n,)? with increasing wave- 
length. This type of decrease was found in the 
visible region by Stumpf and Van Wyk? and can 
probably be extrapolated to the infra-red." 

Since there is no evidence to indicate that the 
indices mp and n, vary discontinuously with 
swarms is of the order of 10~ cm’, G. Foex and L. Royer, 
Comptes Rendus 180, 1912 (1925). W. Kast, Ann. d. 
Physik 73, 145 (1924). 


10 This is especially noticeable in the isochromatics of 


compound I. 
u LL. R. Ingersoll (Phys. Rev. 44, 399 (1933)) reports a 
decrease of this type in the case of the Kerr constant. 
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temperature in any one phase, the arrests and 
minima in the isochromatics may arise from 
polymorphism in the liquid crystalline phase, as 
observed by Dorn and Lohman" their 
measurements of refractive index, the variation 
of arrests and minima with wave-length depend- 
ing upon the variation of (m)—m,)? with wave- 
length. The observed “‘hysteresis” in transmis- 
sion is probably due to a tendency of the mole- 
cules to maintain the random distribution char- 
acteristic of the isotropic liquid phase. As 
mentioned above, the temperature interval over 
which supercooling persists is apparently inde- 
pendent of the cooling rate. 

Attention is called to the fact that the presence 
of the increased transmission of compound III 
in the temperature interval 63°C to 84°C may 
indicate the existence of an intermediate phase 
not previously detected. On cooling to the 
normal solid phase the compound apparently 
does not pass through this state. The presence of 
the intermediate phase may be due to poly- 
morphism in the solid. 

Although the observed variations in band 
shape and position with change of state may be 
attributed to associational effects, it would be 
necessary to study a large number of substituted 
compounds before any definite conclusions as to 
the associational mechanisms could be drawn. 
The isochromatic data indicate pronounced 
variations in transmission with change of state, 
the liquid crystalline phase being more closely 
related to the solid than to the true liquid. The 
present method provides a sensitive means of 
studying change of state and for anisotropic 
liquids is superior to ordinary cooling-curve 
measurements and visual observations. In com- 
mon with methods involving polarized light, 
the present scheme necessitates the use of thin 
layers. If the interpretation in terms of Orn- 
stein’s theory is correct, isochromatic data 
should enable one to estimate the magnitude of 
birefringence at infra-red wave-lengths. 

The writers wish to express their appreciation 
to Mr. J. P. Bain for preparing the compounds, 
to Mr. H. B. Messec for aid in construction of 
apparatus, and to Mr. R. G. Larrick, who as- 
sisted in taking some of the data. 


2E. Dorn and W. Lohman, Ann. d. Physik 29, 533 
(1909). 
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Using a grating spectrometer having a KBr foreprism 
the spectrum of methyl alcohol vapor was studied in the 
region from 2.5 to 26u. This molecule has bands at 3683, 
2978, 2845, 2054, 1477, 1455, 1340, 1034.18 cm~ and a 
very broad band extending from 860 to beyond 380 cm™, 
These bands are typical perpendicular and parallel bands 
and indicate that the molecule is only slightly asymmetric. 

From the fine structure of the 1034.18 cm™ parallel band 
and from certain assumptions about the structure of the 
molecule the two largest moments of inertia, A and B, 
were found to be 35.18 and 33.83 X10~ g cm?, respectively. 
It was not possible to measure C directly but it was esti- 
mated to be approximately 6.8 X10-* g cm?. 

The low frequency band is probably due to the vibration 
of the hydroxyl hydrogen atom perpendicular to the O—H 


(Received June 20, 1938) 


bond and perpendicular to the figure axis of the molecule. 
The presence of this band indicates that the hydroxyl 
group is not free to rotate, at least in the ground state. 
There is evidence of rotation in the structure of the band 
near the high frequency edge. The spacing increases toward 
high frequencies and finally there is a group of lines 
having a spacing of about 40 cm~ between 600 and 860 
cm~!, Since these lines are very weak and occur at such 
high frequencies it appears that free rotation exists only in 
states of high excitation. A quantitative study of free 
rotation in a symmetric molecule has been made. This 
investigation indicates that in a first approximation free 
rotation should give rise to groups of lines having a spacing 
of about 40 cm™ between each group. 


INTRODUCTION 


STUDY of the infra-red absorption spec- 

trum of methyl alcohol has been made by 
Bartholomé and Sachsse! using a prism spec- 
trometer. Earlier studies have also been made by 
Titeica? and Lecomte.’ In addition, the Raman 
spectrum of methyl alcohol has been well es- 
tablished by the work of a number of investiga- 
tors. These studies were of value in establishing 
the positions of the bands but they were unable 
to give information about their character. 
Moreover, due to the absorption of the rocksalt 
prisms used in the infra-red investigations it 
was not possible to work with wave-lengths 
longer than 16y. 

In the present work the absorption spectrum 
was studied under high dispersion using the 
grating spectrometer designed by Hardy.® With 
this instrument it was possible to resolve lines 
having a separation of less than 1 cm=! over most 


'E. Bartholomé and H. Sachsse, Zeits. f. physik. 
Chemie B30, 40 (1935). 
(1935) Titeica, Comptes rendus Acad. Sci., Paris 196, 391 
a : J, Lecomte, Comptes rendus Acad. Sci., Paris 180, 825 


‘J. O. Halford, L. C. Anderson, G. H. Kissin, J. Chem. 


Phys. 5, 927 (1937). V. S. Mizushima, Y. Morina, G. 
Okamoto, Bull. Chem. Soc. Japan 11, 698 (1936). K. W. F. 
Kolhrausch, Der Smekal-Raman Effekt, p. 309. 

* J. D. Hardy, Phys. Rev. 38, 2162 (1931). 
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of the range. By using a KBr foreprism it was 
possible to extend the spectral region to 26uz. 
This extension of the spectral region and the 
ability to resolve the fine structure of the bands 
made it possible to assign the frequencies more 
accurately and to gain further information about 
the structure of the molecule. 

The form of the methyl alcohol molecule is 
known from sterochemical evidence and from 
infra-red data on other similar molecules. The 
carbon atom lies within a tetrahedron formed by 
hydrogen atoms at three of its corners and by the 
hydroxyl group at the fourth. The- hydroxyl 
hydrogen atom stands out at an angle, somewhat 
greater than a right angle, from the C—O line. 

Such a molecule should possess twelve funda- 
mental frequencies. Since the molecule is very 
nearly symmetrical, however, it turns out that 
only nine of the frequencies are distinct, three 
being double perpendicular frequencies. Of these 
nine, six correspond to the frequencies found for 
the methyl halides. Following Dennison’s nota- 
tion,’ the first six frequencies are assigned to 
those modes of vibration which also occur in the 
methyl halides. v;, vg and vy are assigned to the 
frequencies associated with vibrations within the 
hydroxyl group. 


6D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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EXPERIMENTAL RESULTS 


The positions of the bands, the approximate 
amounts of gas required, the gratings used and 
the slit widths are given in Table I. A 25 cm 


TABLE I. Position of the bands. 


GRATING LINES 


FREQUENCY PR 
PER INCH 


SSURE 
G 


860—650 
630—380 


{ 1200 
360 (2nd order) 


absorption cell with KBr windows was used for 
all the bands except for the region between 860 
and 650 cm-!. Here it was necessary to use a 
meter cell to obtain sufficient absorption. It was 
impossible, however, to measure the pressure of 
the gas in the long cell since it could not be 
evacuated. 

v; and ve are very intense bands which have 
their centers at 2845 and 2978 cm-, respectively, 
Fig. 1. These two frequencies are associated with 
the vibrations of the three hydrogen atoms in the 
methyl group, »; being a parallel vibration and 
ve a double perpendicular vibration. Although 
the amplifier was used at its limit of sensitivity 
and the slits were made as narrow as possible it 
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S 


an 
S 


Percent Absorption 
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was impossible to resolve the fine structure of 
either of these bands. 

vs and v4 are also associated with vibrations 
within the methyl group. v3 is a parallel band 
having its zero branch at 1455 cm“. v4 is another 
double perpendicular band which has its center 
at 1477 cm~, Fig. 2. These two bands fall in a 
region where the atmospheric absorption due to 
water vapor is very intense. For this reason it 
was very difficult to resolve the fine structure of 
the bands except for a few lines near the center 
of the perpendicular band. These lines have the 
frequencies: 1463.2, 1471.7, 1479.9 cm~. It was 
impossible to obtain the fine structure of the 
parallel band even with a finer grating and 
narrower slits. 

These four methyl frequencies found in 
methyl alcohol are in good agreement with those 
found in the methyl halides and are almost 
identical with the frequencies of methy] fluoride.’ 
In the spectra of the methyl halides, however, 
another parallel band is usually observed be- 
tween vy; and v2, which has been interpreted as 
the overtone of v4.8 This frequency was not 
observed in the spectrum of methyl alcohol 
although the region was examined carefully with 
different amounts of gas. The perpendicular 
band is not symmetrical, however, and it may 
be that there is a faint band on the low fre- 
quency side of it. 

The double perpendicular band »¢ consists of 
the vibration between the methyl and hydroxyl 
"7 W. H. Bennett and C. F. Meyer, Phys. Rev. 32, 888 


(1928). 
8 A. Adel and E. F. Barker, J. Chem. Phys. 2, 627 (1934). 
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Fic. 1. The bands »; at 2845 cm=! and v2 at 2978 cm™. 


554 2 

SLIT 
WIDTH 
IN J 
cm"! 
vr 3683 
vo 2978 | 
2845 
2054 5 | 
ve 1477 
vs 1455 
ve 1340 
vs 1034.18 
| 
exp 
Nie 
app 
is, 
int 
ture 
y side 
pen 
ing 
tho 
1351 
ture 
| freq 
T 
and 
ban 
j in tl 
fits 
whe 
the 
que; 
ale 
witl 
2800 2840 3040 cm" ban 


SPECTRUM OF CH;0H 


Percent Absorption 


1220 1260 1300 


1340 


1380 1420 1460 1500 


Fic. 2. The bands »; at 1455 cm™, »4 at 1477 cm™ and v¢ at 1340 cm™. 


groups and lies at 1340 cm“, Fig. 2. One might 
expect the two frequencies in this band to be 
separated, as was the case in the corresponding 
perpendicular band of formaldehyde, studied by 
Nielsen. Such a splitting is not immediately 
apparent from the structure of the band. There 
is, however, considerable atmospheric absorption 
out to 1350 cm~' which introduces uncertainties 
in the positions and intensities of the fine struc- 
ture lines. The band does appear to be con- 
siderably broader than most of the other per- 
pendicular bands of this molecule, thus indicat- 
ing the possibility of a double frequency. Al- 
though the background is fairly clear beyond 
1350 cm~! the band does not show much struc- 
ture except near the center and far out on the low 
frequency edge. The central lines have the fre- 
quencies: 1321.3, 1331.8, 1348.4, 1359.5 cm7. 

The parallel vibration v; between the methyl 
and hydroxyl groups gives rise to a very intense 
band at 1034.18 cm-, Fig. 3. The fine structure 
in this band is well resolved in both branches and 
fits the quadratic equation : 


v=1034.18+1.5975m—0.0088m? 


where m is the number of the line counting from 
the center. The calculated and observed fre- 
quencies are given in Table II. Most of the 
calculated values fit the observed frequencies 
within 0.05 cm-!. 

The overtone of »; occurs as a weak parallel 
band at 2054 em-, Fig. 4. The fine structure 


°H. H. Nielsen, J. Chem. Phys. 5, 818 (1937). 


could be resolved only in the negative branch. 
Here the spacing agrees well with the spacing 
found in »;. 

v; consists of the motion of the hydroxyl 
hydrogen atom along the O—H bond and lies 
at 3683 cm~", Fig. 5. The atmospheric absorption 
due to water vapor is particularly strong here 
and it was impossible to obtain more than just 
the envelope of the band. The central line, how- 
ever, occurs between two water vapor lines and 
was observed several times. This band has the 


TABLE II. Frequencies of the fine structure lines in the 
1034.18 cm™ parallel band. 


CALcu- 


OBSERVED OBSERVED LATED 


1031.94 
1032.57 
1034.18 
1035.77 
1037.34 
1038.89 
1040.43 
1041.94 
1043.45 
1044.93 
1046.40 
1047.85 
1049.28 
1050.69 


1038.79 
1040.34 
1041.92 
1043.41 
1044.94 
1046.40 
1047.75 
1049.23 
1050.64 
1052.12 
1053.55 
1054.86 
1056.18 
1057.55 
1058.82 
1060.02 
1061.23 
1062.50 
1063.80 


1029.36 1029.31 
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-il 1015.62 1015.55 14 1054.93 
—10 1017.40 1017.32 15 1056.16 
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va -5 1025.94 1025.98 20 1062.61 
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appearance of a perpendicular band with pos- 
sibly a weak parallel component. 

In the Raman spectrum of liquid methyl 
alcohol v; was found at 3400 cm instead of at 
3683 cm™'. This discrepancy is not an error in the 
results but is characteristic of the molecule when 
it passes from the liquid to the gaseous state. 
Errera and Mollet’ have plotted this band in the 
infra-red using liquid alcohol and found it at 3y. 
When they diluted the alcohol with CCl, this 
band disappeared and a faint one was detected 
at 2.74. They conclude that the 3u band is due to 
association and that the fundamental is at 2.7. 
The dilution of the alcohol breaks up the asso- 
ciation of the molecules. This agrees well with 
the results of the present investigation since no 
band was found at 3 using alcohol vapor. 

Another hydroxyl frequency v9 consists of the 
motion of the hydrogen atom perpendicular to 
the O—H bond and nearly parallel to the figure 
axis of the molecule. This band was never found 
in the infra-red. The Raman spectrum gives this 
frequency at 1107 cm~. In the molecule CH;0D 
this band also is affected by association." Its 
frequency is about 80 cm™ lower than the 
Raman spectrum predicts. If a similar shift takes 
place in ordinary methy] alcohol this band should 


cm 


an 


Percent Absorption 
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Fic. 4. The band 2»; at 2054 cm7. 


10 J, Errera, Trans. Faraday Soc. London 33, 122 (1937). 
G,. Bosschieter and E. F. Barker, paper immediately 
following in J. Chem, Phys, 


Fic. 3. The parallel bands »; at 1034.18 cm. 
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lie near 1030 cm which is in the midst of the 
very intense 1034.18 cm band, v;. The 1034.18 
cm band does have irregularities in its fine 
structure which may be due to the presence of 
this other band. But the hydroxyl frequency is 
essentially parallel and one should expect to 
find a zero branch somewhere. Although the 
region was studied a number of times with 
different amounts of gas no zero branch could be 
found unless it coincides with the zero branch 
of v;. These two frequencies are of the same 
symmetry class, however, and one would ex- 
pect a resonance interaction to separate them. 
It is not possible to use much gas in this region 
before the absorption becomes complete. The 
hydroxyl frequency may be too weak to be 
detected here. 

Finally, there is a frequency vs in which the 
hydrogen atom moves perpendicular to the O—H 
bond and perpendicular to the axis of the mole- 
cule. If there were no restoring force this motion 
would result in the free rotation of the hydroxyl 
group about the axis of the molecule. The ab- 
sorption band between 16 and 26y indicates, 
however, that the hydroxyl group is not entirely 
free, Fig. 6. At 26 the absorption of the KBr 
foreprism made it impossible to work at longer 
wave-lengths although the band probably con- 
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Fic. 5. The band »; at 3683 cm. 
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tinues for some distance. There are a number of 
water lines in this region. These become more 
numerous and more intense as one proceeds 
toward longer wave-lengths. The curve is 
drawn with broken lines where the atmospheric 
absorption makes the readings uncertain. 

Toward the high frequency edge of the band 
the spacing between the lines becomes greater 
and the lines fall off in intensity. Between 650 
and 860 cm! is a group of very faint lines with a 
regular spacing of about 40 cm! which are 
probably a continuation of vs, Fig. 7. Although 
these lines could be observed with the 25 cm 
absorption cell, their structure was better de- 
fined when the meter cell was used. In the region 
between 700 and 630 cm~ there is intense at- 
mospheric absorption due to COs. The line at 
694.7 cm! was strong enough to stand out 
between the COz lines. It appeared to be fairly 
broad but it was impossible to give its exact 
structure. In addition to the fines shown in 
Fig. 7, there is a region of absorption near 650 
cm. This absorption, also, was broad but no 
part of it was strong enough to stand out between 
the CO, lines. The frequencies of the most im- 
portant lines of vg are listed in Table III. 


DISCUSSION OF RESULTS 


The energy of an asymmetric molecule has the 
form :2 
he? 


E=(n+})hv+— 
82? 


where W is a parameter which makes the energy 
matrix diagonal and is not in general an integer. 
When the values of W are computed for the 
methyl alcohol molecule it turns out that only in 
the lowest states do they differ appreciably from 
squares of integers. This results in a displacement 
and splitting of the lowest levels but the effect is 
negligible for values of W greater than 4. This 
splitting affects the structure of the perpendicular 
bands slightly but, except for a slight broadening 


of the lines, it is unimportant in the parallel 
bands. 


*i. Hi. Nielsen, Phys. Rev. 38, 1432 (1931). 


1 1 1 


Fic. 6. The perpendicular band ys. 


From the fine structure of the 1034.18 cm- 
parallel band an average value of the two largest 
moments of inertia of the molecule could be 
obtained. From further assumptions about the 
internuclear distances, using infra-red data on 
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TABLE III. Lines of the perpendicular band vs. 


382.45 474.65 cm™ 582.60 
383.20 484.60 599.50 
384.65 486.05 611.45 
392.15 494.95 617.55 
401.45 500.80 623.90 
416.10 506.40 (650) 
439.70 515.73 694.7 
444.30 516.42 736.7 
452.20 539.60 778.3 
460.55 540.90 780.0 
462.20 548.65 820.6 
463.38 574.50 860.1 


the methyl halides'* and the water molecule," it 
was possible to compute A and B with a fair 
degree of accuracy. A =35.1810-*° g cm? and 
B=33.83X10-* The C—H and H—H 
distances in the methyl group were assumed to 
be 1.11A and 1.81A, respectively. The O—H dis- 
tance was taken to be 0.9558A as in water. In 
their experiments on electron diffraction Pauling 
and his co-workers give 1.44+0.02A for the 
C—O distance in CH;ONHz:. Results from the 
moments of inertia of methyl! alcohol indicate 
that the C—O distance in methyl] alcohol is 
between 1.42 and 1.43A and the angle between 
the O—H and the C—O bonds lies between 
105° and 115°. These values cannot be known 
exactly until a more exact measurement of the 
least moment of inertia C is made. Although the 
present observations do not yield C directly its 
value can be estimated from the dimensions of 
the molecule. In making computations on the 
free rotation of the two parts of the molecule it is 
also necessary to know the approximate values of 
the moments of inertia of the hydroxyl and 

. M. Johnston and D. M. Dennison, Phys. Rev. 48, 868 
ars ft M. Randall, D. M. Dennison, N. Ginsburg, L. R. 
Weber, Phys. Rev. 52, 160 (1937). 


_ L. O. Brockway, J. Y. Beach, L. Pauling, Am. Chem. 
J. 57, 2693 (1935). 
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respect to a set of fixed axes is best described by 


methyl groups separately about the figure axis of F 
the molecule. For the hydroxyl group C; is 
between 1.3 and 1.4 10-* g cm?; for the methyl 
group C.=5.4X10-" g cm*. The moment off 
inertia C is equal to the sum of C; and Cz and 
lies between 6.7 and 6.8 g cm?. . 

The perpendicular band »s is a single frequency F 
and should be unaffected by the internal angular | 
momentum which changes the spacing in the > 
double perpendicular bands of symmetric mole- 
cules. A simple calculation shows that if the 
hydroxyl group were completely bound one J 
should expect an ordinary perpendicular band > 
having a spacing of approximately 6.5 cm. If, 
however, the hydroxyl group were a free rotator 
one should expect a series of rotation bands. 
with a spacing of about 40 cm~ between each 
band. The low frequency band shows an irregular | 
spacing which increases from about 10 cm™ to J 
30 cm™ in the region from 380 to 600 cm™ | 
and finally has a spacing of about 40 cm~ be- 
tween 600 and 860 cm~'. The increased spacing 
toward higher frequencies probably indicates 
that for higher energies the hydroxy] group is 
approaching rotation levels. 


FREE ROTATION BANDS OF METHYL ALCOHOL 


The problem of the torsion oscillator-rotator 
has already been treated by Nielsen." For sim- 
plicity it will be assumed that the hydroxyl 
group can be approximated by a symmetrical 
top. Then the whole molecule can be represented 
by two symmetrical rotators having a common ( 
axis. It will be convenient to refer to the two 
systems with the subscripts 1 and 2. The posi 
tions of the two parts of the molecule with 


1° H. H. Nielsen, Phys. Rev. 40, 445 (1932). 
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Fic. 7. Rotation lines of the band vs. 
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the Eulerean angles, 0, and 0, g2; where 


fis the angle between the fixed Z axis and the C 


axis of the molecule, y is the angle which de- 
scribes the precession of the molecule, ¢; and ¢2 
are the angles of rotation about the C axis. 
Following Nielsen’s procedure it will be ad- 


vantageous to introduce the quantities g and a 


defined as follows: a=¢e2 
- —g, where C; and C2 are the moments of 


inertia of the two parts of the molecule and 
C:+C2=C, the total moment of inertia about 


_ the figure axis of the molecule. Since the potential 
_ energy depends only upon the angular displace- 
- ment between the two parts of the molecule it 


isa function of @ only. 

Proceeding as Nielsen did the Schrodinger 
equation is separable, one part giving the energy 
of rotation of the molecule as a whole: 


h? 
A C A 


The other part leads to the differential equation : 
@M 
da? h?C 


{ E2— V(a)} M=0. 


Since the hydroxyl hydrogen atom has three 
positions of equilibrium between the hydrogen 
atoms of the methyl group the simplest potential 
function is: V= Vo (1—cos 3a). This leads to the 
Mathieu equation. As was pointed out by Nielsen 
exact solutions for this equation are known only 
for every third state. The intermediate solutions 
can only be estimated. 

It is of interest to investigate the solutions 
for extreme values of the potential. If the barrier 
is indefinitely high there will be simple harmonic 
motion. This will give rise to an ordinary perpen- 
dicular band with a spacing (h/4n?)(1/C—1/A). 
If there is no barrier, however, there will be free 
rotation of the two parts of the molecule. Then 
the energy of the molecule may be written : 


h? | 24+4J—K? 
Sr? A 


The eigenfunction has the form: 
P= F(O)e reer, 


where F(6) is the hypergeometric function. The 
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notation used here is somewhat different from 
that employed by Nielsen. \; and 2 are the 
azimuthal quantum numbers for the two parts 
of the top and may take on all integral values. 
K is the total azimuthal quantum number and 
|K| = 

The selection rules are similar to those obtained 
for ordinary perpendicular bands. If the electric 
moment changes in the system denoted by the 
subscript 1: AJ=0, +1, AK=+1, Ad\i=+1, 
Ad\2=0. If the electric moment changes in the 
other part of the molecule, then Ad\,;=0 and 
Ad\:= +1. In the following discussion it will be 
assumed that the electric moment changes in the 
hydroxyl group, which will be assigned the sub- 
script 1. 

Using these selection rules the zero branches 
on the low frequency side of the center of the 
band have the frequencies: 


Ari KK OJ Ar? A 
where =| K]. 


The fine structure lines associated with the 
Kth zero branch have the frequencies: 


High frequency side of the center 
Ar K J Ar? C. A 
|J|=|K|. 


where 
Low frequency side of the center 

h K-34J7 
 A | 
|J-1|=|K]. 


K-1J7 
Ao iK 


where 


Since J is never less than K certain of the 
fine structure lines will be missing. It is the 
same situation that arises in perpendicular bands 
of symmetric molecules and the band is con- 
structed in the same way. On the low frequency 
side the first zero branch has the first rotation 
line on the low frequency side missing. The 
second zero branch has two lines missing on the 
low frequency side and one on the high frequency 
side. On the high frequency side the situation is 
just the reverse. A typical rotation band showing 
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how the zero branches and fine structure lines 
go together is drawn in Fig. 8. The asymmetry 
of the band will be explained later. 

It is apparent from these equations and the 
diagram that the fine structure lines and the 
zero branches have the same spacing and in most 
cases a number of lines coincide. In computing the 
intensities it is necessary to find the sum of 
the intensities of all the lines falling at a particu- 
lar place. It will be convenient to find the in- 
tensities of the zero branches and fine structure 
lines separately. The intensity formulae for 
perpendicular bands given by Dennison® may be 
adapted for use here. For the zero branches on the 
low frequency side of the center the intensity 
formula is: 


(2J+1)(J+K)(J-—K+1) 
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1A2K— 
Al AWK 


=A[x] 


J=K 


2J(J+1) 
Xexp 
h? A A 


where ’ 
C; C2 


A is a constant which depends upon the popula- 


3 
+10 


| 


Fic. 8. The free rotation band \,=6 
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tion in the initial state and the Boltzmann factor. 
The quantity [x] is equal to one except when 
the quantum number ), is zero or any multiple 
of three, in which case [x] takes on the value 
two. This enhancement of every third line is 
due to the presence of the three protons (like 
particles with spin }(4/27)) in the methyl group 
which doubles the population in every third 
state. 

In the approximation needed here it is con- 
venient to wn the sum into an integral: 


Ai—1A2 K— 
Al 


Xexp (K —1)*]) 
e-o(K2+K) | 


(K*- K) —dal. 


o o(K°+K) U 


The integral which occurs here is the logarithmic 
integral for which tables have been computed. 
A similar expression is obtained for the zero 
branches on the high frequency side: 


M-1eK J A 
Di 


Xexp K*]) 


o o(K°—K) U 


The intensities of the fine structure lines can 
best be treated in two parts. First consider the 
fine structure associated with zero branches on 
the opposite side of the center of the band, and 
then consider the fine structure associated with 
zero branches on the same side of the center. 
In either case it is necessary to sum over all 
rotation lines of that class falling on the Kth 
zero branch. 

Consider all rotation lines on the high fre- 
quency side of the band having zero branches on 
the low frequency side, which fall on the Kth 
zero branch on the high frequency side. The Kth 
zero branch has the frequency : 


Ci A 


where K—1=—(Ai+2). The fine structure line 
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associated with the K’th zero branch on the low 
frequency side has the frequency : 


h 
Ci A 
where K’=)i+)2’ and |J|=|K’|. It is clear 
that this fine structure line coincides with the 
Kth zero branch if K’=J—K+1 and d\.’=J—K 


+1—;. Each fine structure line has the in- 
tensity : 


de! K’—1 
K’ 2 


’ 


(J+K'—1)(J+K’) 
2J 
Xexp (—o[J?—J—(K’—1)? 
1)? +-y2d2"? ]) 
(2J—K)(2J—K-+1) 
2J 
Xexp (—oy2(J+a)’), 


A,-1 do’ K’-1 J-1 
= x 
he’ K’ J A[ 


= A[x 


where 
2K —1—2y2(K+)A,-—1) 
a= 
272 


ye (K +1 


The intensity of all the fine structure lines falling 
on the Kth zero branch may be found by sum- 
ming this expression over J, from J=K to 
infinity. It will be more convenient, however, 
to use the integral form, where 4/3 has been 
written for the average value of [x ]. 

2J 


Xexp (—oy2(J+a)*)dJ. 


4 ea 
I=-Ae** f | 0K -1+20)+ 
K 


The first term may be integrated at once and the 
other two may be written in a simpler form: 


4 
3 


2K —1+2a 
(o72)*(K +a) 


K?-K 7” 
2 +a) X+(o72)!a 


The second integral is the probability integral 
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for which tables have been computed. The third 
may be evaluated numerically. 

Exactly the same expression may be obtained 
for the intensities of the fine structure lines on 
the low frequency side which are associated with 
zero branches on the high frequency side. Here, 
however, 


2K —1—272(K 
a= 


272 


and 
1)? (K —A1)?—a? ]—(K 


In some cases a is negative and greater than K. 
Here the integral does not well approximate the 
sum. It is then more convenient to sum the first 
few terms until J=a. The integral form may be 
used for the remaining terms. 

Now consider all the fine structure lines on the 
low frequency side associated with zero branches 
also on the low frequency side. The intensity of 
one of these lines which is associated with the 
K'th zero branch is: 


(J—K’'+1)(J—K’) 

2J 
Xexp (—o[J?+J—(K’—1)? 


A fine structure line associated with the K’th 


Ai—1A2’ K’-1 J 
Ai 2’ K’ 


TABLE IV. Intensities of lines in the rotation band \,=6. 


Low FREQUENCY SIDE HiGH FREQUENCY SIDE 


FINE St. FINE St. 
K Branco ZB HiGu Z B Low 
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FINE St. FINE St. 
K Brancu Z BLow Z B Hicu 


12.60 17.21 
4.50 16.18 
3.03 15.23 
3.89 14.41 
1.20 13.45 

12.82 


11.37 
9.82 
8.44 
6.66 
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zero branch will coincide with the Kth zero 
branch if K’=K—J and \o’=K—J—),. It is 
only high K values, Fig. 8, which receive con- 
tributions from this class of fine structure since 
J=K'+1, or J=(K+1)/2 if K is odd and 
J=(K+2)/2 if K is even. Also, since J=K—K', 
only those values of /=K—1 can give contribu- 
tions to the Kth zero branch. The intensity of all 
these fine structure lines falling on the Kth zero 
branch is: 


K-1 (2J —K)(2J—K-+1) 


J=(K+1)/2 K odd. QI 
J=(K+2)/2 K even. 


Xexp (—ov2(J+a)*) 


2K —1—2y2(K 
where a= 


272 
and 
1)? (K —1)?—a? J— (K—1)?. 


The same sum is obtained for the intensities of 
the fine structure lines on the high frequency side 
of the band which are associated with zero 
branches also on the high frequency side. In this 
case, however, 


2K 
a= 


272 
and ]—K?. 


Since it is only in states of high excitation 
that the hydroxyl group approaches free rotation 
the intensities of some of the lines in the band 
where \,;=6 are computed in Table IV. This 
band is drawn in Fig. 8. The quantum number 
K changes from K—1 to K for the lines on the 
low frequency side of the center of the band and 
from K to K—1 for the lines on the high fre- 
quency side. In each rotation band there are two 
particularly strong zero branch lines having the 
quantum number \,.=0 and —3. In the first 
rotation band where \i=1 these strong lines 
would occur at the center of the band and the 
fine structure lines would be distributed evenly 
on both sides. In each successive band, however, 
these strong lines are displaced farther from 
the center toward lower frequencies by an 
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amount h/(47°A). In high rotation states this 
displacement brings most of the important zero 
branches on the low frequency side of the center. 
Thus the important contributions from the fine 
structure in the high rotation bands is that 
associated with zero branches on the low fre- 
quency side. These fine structure lines always 
begin at a distance h/(87?A) from the center 
of the band on the high frequency side but they 
build up more slowly on the low frequency side. 
For this reason the bands have an unsymmetrical 
appearance and are very intense on the high 
frequency side since so many fine structure lines 
coincide here. This effect is enhanced in each 
successive band as ), increases. 

These bands do not compare very well with 
the structure of the lines in Fig. 7. These are not 
single lines, however. It may be that, when 
interactions between the two parts of the mole- 
cule are considered, the rotation bands will have 
quite a different form. 

The results of the present investigation are 
highly suggestive from a qualitative point of 
view. It has been shown that, if the motion of 
the hydroxyl group were entirely free, one should 
expect a series of rotation bands with a separa- 
tion of about 40 cm-! between each band. This 
is clearly not the case here. The very fact of the 
existence of a low frequency band shows that the 
levels must in some measure possess a vibration- 
like character. This indicates that the hydroxy! 
atom is confined within a threefold potential 
barrier produced by the three hydrogen atoms 
of the methyl group. If this barrier were very 
high the particle would vibrate in one of its 
positions of equilibrium and give rise to a 
typical perpendicular band. This band would 
not be affected by vibration-rotation interaction 
and should have a spacing of about 6.5 cm". 
The spacing in the low frequency band is be- 
tween these two extremes. At the low frequency 
edge of the band the spacing approaches the 
bound case but near the high frequency edge the 
distance between the lines indicates that the 
hydroxyl atom is nearly free. 

While the evidence is not at all conclusive, it 
suggests that the height of the barrier might be 
estimated by correlating it with that part of the 
band where the distance between the lines begins 
to approach the spacing expected on the basis 
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of free rotation. This would appear to corre- 
spond to a quantum number of the order of 
magnitude of 5 in the energy of the free hydroxy] 
rotator. This estimate of the barrier is accord- 
ingly between 400 and 600 cm-'. A study of the 
far infra-red would be helpful in fixing this 
value. 

It is beyond the scope of the present paper to 
discuss the transition between vibrational and 
rotational states. Since the potential function of 
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the alcohol molecule is intermediate between the 
two extremes discussed here, it is clear that both 
the rotational and vibrational states are com- 
plicated by the presence of the three minimum 
potential barrier. 

In conclusion the authors wish to express 
their gratitude to Professor D. M. Dennison for 
his interest in the work and for many helpful 
discussions during the development of the 
theory. 
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The Infra-Red Absorption Spectra of CH;0D and CH,.DOD 


E. F. BARKER AND G. BOSSCHIETER 
Department of Physics, University of Michigan, Ann Arbor, Michigan 


The infra-red absorption spectra of CH;OD and CH2DOD 
between 2.54 and 244 have been examined with a KBr 
prism spectrometer, and with appropriate gratings. The 
observed bands represent all of the fundamental vibrations 
except the one of lowest frequency which is associated with 
torsional vibrations about the C—O bond. Since these 
molecules depart only slightly from axial symmetry, the 
bands, with the exception of three due to the hydroxyl 
radical, correspond in position and appearance to those of 
the methyl halides. The rotational structure for the 10u 
band (vs) of CH;OD has been resolved, and partial resolu- 
tion is obtained in some other bands. The band ys in 


HE Raman spectra of ordinary methyl alco- 
hol in the liquid phase, and of two methyl 
alcohols containing deuterium, have recently 
been described by Halford, Anderson and 
Kissin.! While the particular modes of vibration 
to be associated with each observed frequency 
may be decided almost by inspection, it is only 
through measurements in the infra-red that the 
character of the motions is revealed. Moreover 
the true molecular frequencies must be observed 
for the gaseous state. A comparison of infra-red 
data for the gases with Raman data for the 
liquids yields information of interest particularly 

in connection with the problem of association. 
The molecule of methyl alcohol departs 
slightly from axial symmetry, since the line 
‘ Halford, Anderson and Kissin, J. Chem. Phys. 5, 927 


{t807). See also Mizushima, Morino and Okamoto, Bull. 
hem. Soc. Japan 11, 698 (1936). 


(Received June 23, 1938) 


CH2DOD has two components arising probably from two 
forms of the molecule in which the hydroxyl D atom occu- 
pies different valleys of the threefold potential. The defor- 
mation vibration (v7) is single for CH;0D but has four 
components in CH,DOD, indicating a separation of levels 
which for the former molecule are degenerate. A comparison 
of the frequencies obtained with gaseous and liquid 
samples indicates large displacements of the bands arising 
from the hydroxyl valence and deformation vibrations, the 
former toward greater wave-lengths in the liquid, and the 
latter toward smaller wave-lengths. 


connecting the atoms O and H is nearly per- 
pendicular to that connecting C and O. The 
minimum moment of inertia corresponds to an 
axis which is very nearly parallel to the C—O 
line, and which, for convenience, may be re- 
ferred to as the symmetry axis. The relatively 
large electric moment may be attributed prin- 
cipally to the hydroxyl group. Twelve funda- 
mental vibrations are to be expected, some of 
which occur in pairs, corresponding to two 
dimensional oscillators with nearly equal fre- 
quencies in the plane determined by the atoms 
C, O and H, and normal to it. Six of these are 
associated primarily with the methyl radical, 
three with the C—O bond, and three with the 
hydroxyl group, the first nine being conveniently 
designated as in the methyl halides.” 


2 Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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The infra-red spectrum of CH;OH has been 
reported in the preceding paper.’ The observa- 
tions here presented include certain bands of 
particular interest in the spectra of CH;OD 
and of CH,DOD. We are indebted to Professors 
Halford and Anderson for the preparation of 
these materials; the specimens are in fact the 
same ones upon which the Raman observations 
were made. Professor Randall kindly placed at 
our disposal a recording prism spectrograph for 
exploratory measurements, the absorption cell 
being a metal tube with windows of KBr and an 
optical path length of 80 cm. For measurements 
at high resolution a prism-grating spectrometer 
was used, with a 60° NaCl foreprism, and grat- 
ings ruled 4800 and 2400 lines per inch. With this 
instrument the absorption cells were 25 cm glass 
tubes without stopcocks, the windows being 
attached with Apiezon wax. A side arm held the 
liquid sample in sufficient quantity to saturate 
the chamber at room temperature ; when intense 
absorption bands were under observation the 
vapor pressure was reduced by cooling this 
reservoir. 


THE METHYL VIBRATIONS 


In the region of 3.54 the two methyl! frequen- 
cies »; and v2 occur. These consist primarily of 
motions in which the hydrogen atoms approach 
and recede from one another. In »; the change in 
electric moment is along the methyl axis, and 
the band is of the parallel type. It appears at 
2850 cm-! in CH;OD, and at 2874 cm“ in 
CH.DOD. The displacement toward higher 
frequencies in the latter case in spite of increased 
mass is caused by a slight change in the direction 


3 Borden and Barker, J. Chem. Phys. 6, 553 (1938). 
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Fic. 1. Transmission at 3.44 and 4.8u. a, vi and v2 for CH;0OD, 25 cm cell at —30°C. 


of vibration. The second fundamental vibration 
v2, is perpendicular in character. In CH;OD it 
involves two oscillations normal to one another 
and of almost identical frequencies. The band 
center, at approximately 2965 cm™, can be 
determined only to within three or four wave 
numbers, since almost no structure is apparent. 
Whether or not there is also absorption in this 
region due to the parallel component of 2y,, as 
indicated by Raman scattering, cannot be de- 
termined because of the breadth and intensity 
of the perpendicular band. The observations are 
represented by curve a of Fig. 1. The two corre- 
sponding motions in CH2,DOD are of very 
different frequencies. One is only slightly dis- 
placed, giving rise to a band at about 2944 
cm, but the second appears at about 2180 
cm~'!. These are shown on curves band ¢ of Fig. 1. 
In the latter curve there appears also a fairly 
intense zero branch at 2158.4 cm which may 
be the harmonic 2»;, although in CH;OD the 
latter appears at 2065.4 cm~, as indicated in 
curve d. 

The second pair of methyl frequencies, v3; and 
vs, lie in the region between 6.8 and 8.0u, and are 
difficult to locate because of the strong atmos- 
pheric absorption due to water vapor. The 
frequencies v3 are 1459 cm! for CH;OD and 
1357 for CH2.DOD, corresponding to vibra- 
tions in which the three hydrogen atoms move 
in phase approximately at right angles to the 
carbon bonds. In v4 the motions are similar ex- 
cept that the three hydrogens oscillate with 
phases differing by equal intervals, producing 
bands of the perpendicular type. In the spectrum 
these appear as groups of lines in the neighbor- 
hood of 1480 cm with a second group appearing 
in CH,DOD near 1330 cm“. 


2075 2060 cM” 


b, v, and v2 for CH,DOD, at —30°C. c, v2’ for CH»,DOD, at 20°C. d, 25 for CH;OD, 


at 20°C. 
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Fic. 2. Transmission at 10u, showing v;. a, CH;0D (with some CH;OH). 8, CH.DOD. c, CH;OH. 


THE C—O VIBRATIONS 


In the methyl halides two frequencies are 
associated with the single carbon bond, one (vs) 
giving rise to a parallel band and the other (v6) 
to one of the pependicular type. The absorption 
due to v; is very intense in the methyl alcohols, 
and is indicated in Fig. 2. With CH;0D two 
zero branches were observed, one at 1034 cm™ 
and the other at 1040 cm~. The fact that the 
former coincides exactly with the band of 
CH;OH as indicated by the broken line, sug- 
gests that this sample had been contaminated 
with ordinary hydroxyl hydrogen, the tendency 
to exchange in this position being well known. 
Observations at 2.74 showed that the OH band 
was present also. The second maximum at 1040 
is for CH;OD. It is somewhat broader 
than the one at 1034, indicating that the replace- 
ment of the hydroxyl hydrogen by deuterium 
produces a molecule somewhat farther removed 
from axial symmetry. The increase in frequency 
again indicates a larger restoring force due to the 
change in direction of vibration. 

In the absorption curve for CH2DOD there are 
also two maxima which appear to be rather broad 
zero branches at 1023 and 1049 cm~'. These are 
probably to be explained by the fact that this 
compound exists in two forms depending upon 
the azimuthal position of the hydroxyl D with 
respect to the methyl D. In approximately two- 
thirds of the molecules the hydroxyl D will lie 
between an H and a D of the methyl group, 
the obliquity of the vibration with respect to the 
C—O line with being even greater than in 
CH;OD, and the frequency correspondingly 


Discontinuities indicate changes in prism setting. 


higher. The fact that no maxima appear in this 
spectrum at 1034 and 1040 cm indicates that 
no exchange of hydrogens occurs in the methyl 
group, and also demonstrates the effectiveness 
of the procedure employed for the preparation 
of the material. 

The band obtained with CH;OD shows a dis- 
tinct rotational structure, somewhat falsified, no 
doubt, by lines due to CH;OH. As shown in 
Table I, the frequencies agree reasonably well 
with those computed from the equation 


v= 1040.61.58: —0.008 m?, 


where m is the ordinal number of any line count- 
ing from the center. This indicates a mean 
moment of inertia for axes normal to the sym- 
metry line, of approximately 35.0 10-* g cm’, 
a little more than 1 percent greater than for 
CH,;QH. 

The perpendicular band vg has not been ob- 
served for CH;OD. In CH:DOD two broad 
maxima consisting of groups of lines were located 
at about 1330 and 1270 Whether the 
higher frequency group is a component of v¢, 
which should be double, or belongs to v, is 
not clear. 


THe O—H anp O—D VIBRATIONS 


The characteristic valence oscillation of the 
hydroxyl radical (v;) occurs at approximately 
2.74 for OH, and at 4.04 for OD. The former is 
indicated in Fig. 3 as an absorption maximum 
around 3680 cm~ in the lower curve (C) taken 
with low dispersion, using CH;O0H vapor. The 
determination of this band under high dispersion 
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TABLE I. Observed and computed frequencies for the band vs. 


v COMP. 


1054.2 
1052.8 
1051.3 
1049.8 
1048.4 
1046.8 
1045.3 
1043.8 


1030.8 
1029.1 
1027.4 
1025.8 
1024.2 
1022.2 
1020.4 


m v OBS. 


1053.7 
1052.5 
1051.3 
1050.0 
1048.2 
1047.0 
1045.8 
1044.0 


1030.5 
1028.7 
1027.7 
1025.3 
1023.5 
1021.8 
1020.1 


is a matter of some difficulty because of the 
intense absorption due to atmospheric water 
vapor, which is indicated in curve A. Measure- 
ments taken under the same conditions with 
CH;OD in the cell are recorded in curve B. 
The very considerable absorption is evidence 
of the presence of hydroxyl H. 

The region in which the O—D valence vibra- 
tion occurs is a much more favorable one, since 
no atmospheric absorption interferes. Fig. 4 
shows the band observed with CH:DOD, its 
center having the frequency 2718.8. An almost 


3745 3640 


Fic. 3. Transmission at 2.74. a, Atmospheric water. 
6b, CH;OD (with some CH;OH). c, CH;OH. 
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identical curve was obtained for CH;OD, with 
the center at 2720.4 cm“. The change in electric 
moment associated with this motion has com- 
ponents parallel to and perpendicular to the 
symmetry axis, both of course having the same 
frequency. The perpendicular one is probably 
greater in magnitude, and the band associated 
with it should exhibit a strong central maximum 
consisting of a family of zero branches spaced at 
intervals of approximately 5 cm~'. These zero 
branches correspond to AK=+1, where K 
measures the angular momentum about the 
symmetry axis. There is some indication of such 
structure, but superposed upon it are a set of 
lines with spacing appropriate to the parallel 
band, as in v;. The small number of rotation lines 
observed on each side is evidence that this 
parallel band is superposed upon a strong and 


—az150 


2780 


Fic. 4. Transmission at 3.74 for CH,DOD. Hydroxy! 
valence vibration »;=2718.8 


rather narrow background, probably attributa- 
ble to the perpendicular transitions. 

A second degree of freedom (vs) for the hy- 
droxyl D atom involves motions of low fre- 
quency normal to the D—O-—C plane. In the 
absence of restraints this would become free 
rotation on the single carbon band. While the 
evidence here presented is not conclusive, it 
appears probable that free rotation occurs in only 
a relatively small number of molecules which 
have been excited to upper states of this motion, 
and that the normal state is an oscillation in one 
of the potential minima provided by the methyl 
radical. There is at least the negative argument 
that free rotation should introduce low states of 
considerable population, and in consequence, 
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combination bands which could hardly escape 
detection. For CH;OD the three potential valleys 
are identical, and the lowest vibration state 
should be single while the second is double but 
degenerate. In CH2DOD, however, one of the 
valleys differs from the other two, and this 
degeneracy would be removed, so that two low 
frequency bands should appear. Prism curves 
including the region from 13 to 22y, taken with 
the 80 cm cell, indicated no absorption by 
CH,OD. In CH:DOD two very weak bands 
appeared near 16u; these, however, could not 
be observed with the grating instrument in 
which it was necessary to use the shorter cell. 
The so-called hydroxyl deformation vibration 
is normal to the O—D line, but in the D—O—C 
plane. This introduces a change in electric 
moment with its major component parallel to 
the symmetry axis. The resulting bands have 
strong zero branches, somewhat broadened, and 
degraded toward higher frequencies. They are 
shown in Fig. 5. The CH;OD gives a single band 
at 863 cm~'. In CH,DOD four frequencies are 
observed, the two different types of molecules 
previously mentioned each being responsible for 
a pair of bands. The frequencies are 952.9, 
892.3, 829.5, and 810.1 cm. It might be ex- 
pected that the splitting of the excited state 
would be greater in the more prevalent than in 


the less prevalent molecular type, the latter 
having a plane of symmetry, but it is not yet 
clear how these four bands are to be assigned. 


THE EFFECT OF ASSOCIATION 


Condensation from the vapor to the liquid 
state modifies some of these frequencies very 
considerably. This effect is apparent from a 
comparison of Raman and infra-red data, and 
has also been verified by observations upon the 
absorption of liquids. The methyl bands, and 
v; due to the C—O bond, are not appreciably 
shifted, but the hydroxyl frequencies show a 
large displacement. The O—H valence vibration 
vz, found at 3650 cm~! in CH;OH vapor is 
reduced in frequency to about 3300 cm™ in the 
liquid, as has been shown by Errera and others. 
The corresponding bands for CH;OD and 
CH:2DOD are shifted from 2720 cm~ to about 
2500 cm-'. The deformation bands y» in the 
deuterium alcohols move to higher frequencies in 
the liquid phase, from 863 cm to about 940 
cm for CH;OD.* This effect could not be 
verified for CH;OH, since vy has not been ob- 
served in the vapor. Its position is indicated by 
Raman measurements as 1104 cm~ in the liquid, 
and in the vapor it is probably obscured by the 


4G. Bosschieter, J. Chem. Phys. 5, 992 (1937). 


ths 


Fic. 5. Transmission at 124 showing vs. a, CH;OD. b, CH,DOD. 
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intense 10u band. A considerable intermolecular 
influence upon vs might also be expected in the 
liquid. A bond between hydroxyl groups in 
adjacent molecules’ would result in a further 
separation of the two components of this motion, 
though one of the components might be ex- 
cluded from Raman transitions by considera- 


5 W. H. Zachariasen, J. Chem. Phys. 3, 158 (1935). 
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tions of symmetry. The only line in the Raman 
spectrum which can be assigned to v6 has a 
frequency of approximately 1160 cm~ for both 
CH;OH and CH;0D, as compared with about 
1340 cm~ for the gases. 

This investigation owes partial support to 
the Faculty Research Fund of the University 
of Michigan, which assistance is gratefully 
acknowledged. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


The Free Energy of the Producer-Gas Reaction, and the 
“Zero-Point’’ Entropy of Graphite! 


In 1933, calculations made by Gordon? on the change in 
free energy for the producer-gas reaction indicated the 
possibility of the existence in graphite of a ‘‘zero-point”’ 
entropy of about 0.5 calorie per degree per mole. In an 
extensive study reported in 1934, Terebesi* compared, for 
five different reactions involving graphite, the experi- 
mentally determined equilibrium constants with those 
evaluated from calorimetric data and statistical calcula- 
tions. According to Terebesi, the data on three of these 
reactions pointed strongly toward the nonexistence of a 
“zero-point’’ entropy for graphite, while the data on the 
other two reactions, one of them the producer-gas reaction, 
were inconclusive. 

New calorimetric data on the heats of combustion of 
graphite’: > yield new values, significantly different from 
the older ones, for the heats and free energies of formation 
of carbon dioxide* and carbon monoxide.’ It seems de- 
sirable, therefore, to calculate a new value for the heat and 
free energy of the producer-gas reaction, and to ascertain 
the effect of the new values on the question of the ‘‘zero- 
point’’ entropy of graphite. 

The new values for the heat and free energy of formation 
of carbon dioxide are, per mole:® 


C (c, graphite) +O2 (g) =COs (g), (1) 

ATT 16 = —393,355 +46 NBS international 
joules, or —94,030+11 calories, (2) 

AF 16 = —394,228+58 NBS international 
joules, or —94,239+14 calories. (3) 


For carbon monoxide, the writer’ gives the following 
recalculated value for the heat of combustion, per mole: 


CO (g)+3 Oz (g) =COs (g), (4) 
AF 295 1g = — 282,937 +120 NBS international 
joules, or —67,635+29 calories, (5) 


and for the heat and free energy of formation, per mole: 


C (c, graphite) + Oz (g) =CO (g), (6) 

AH 16 = —110,418+129 NBS international 
joules, or —26,395+21 calories, (7) 

AF x95 16= —137,160+133 NBS international 
joules, or —32,788+32 calories. (8) 
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Combination of the foregoing values, together with 
appropriate values®:7 of H%295.i16s—H%, yields for the 
producer-gas reaction, per mole: 


C (c, graphite) +COz (g) =2CO (g), (9) 

AFP 16 =172,519+244 NBS international 
joules, or 41,240+58 calories. (10) 

= 119,908+248 NBS international 
joules, or 28,663 +59 calories, (11) 

AH, = 165,603 +244 NBS international 
joules, or 39,587 +58 calories. (12) 


Combination of the above value for AH% with the 
appropriate values of (F°—H%)/T, reported by Clayton 
and Giauque,’ Johnston and Davis,!° and Kassel," permits 
calculation of the equilibrium constant for reaction 9 for 
any given temperature: 


logio K = —(1/2.3026R)[AH%/T+A(F°—H%)/T]. (13) 


In this manner, there may be calculated the following 
values of logio K : 1000°K, 0.2689; 1100°K, 1.0772; 1200°K, 
1.7467; 1300°K, 2.3093; 1400°K, 2.7891; 1500°K, 3.2022. 
These values of the equilibrium constant are found to be in 
excellent accord with the experimentally determined 
equilibrium constants within the relatively large limits of 
uncertainty of the latter (see the summaries by Eastman," 
Chipman,” and Terebesi*). It may be concluded, therefore, 
that the data on the producer-gas reaction are in accord 
with the nonexistence of a ‘‘zero-point'’ entropy for 
graphite. 

FREDERICK D. RossINI 

National Bureau of Standards, 


Washington, D. C., 
August 10, 1938. 
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